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PREFACE 


--  The  following  report  is  the  last  of  a  series  of  five 
reports  dealing  with  research  related  to  the  problem  of 
deep-sea  mooring  lines » under  the  stated  contract.  These 
studies  were  initiated  in  January,  1959  with  Dr.  Basil  W. 
Wilson  as  principal  investigator.  A  large  part  of  the 
efforts,  as  summarized  in  the  first  four  reports,  were 
related  to  the  problem  of  the  equilibrium  configuration. 

The  present  study,  which  was  been  completed  under  the  support 
of  the  Texas  A  &  M  Research  Foundation  in  fulfillment  of  the 
original  contract  obligations,/)  deals  with  the  dynamics  of 
deep-sea  mooring  lines.  The  primary  emphasis  is  on 
methodology  in  the  solution  of  this  problem.)  Although 
some  application  is  made,  an  extensive  parametric  study 
is  well  beyond  the  scope  of  the  present  study. 

Several  publications  based  upon  prior  studies  under 
this  contract  have  appeared  in  the  literature  and  are 
cited  in  the  list  of  references. 

The  author  wishes  to  express  his  appreciation  for 
the  patience  of  the  original  sponsor.  Hs  is  particularly 
grateful  for  the  preserverance  of  his  secretary,  Mrs. 

Florence  G.  McCnlly  in  coping  with  the  lengthy  set  of 
equations  endlessly  interwoven  into  the  fabric  of  the 
presentation.  It  is  hoped  that  the  reader  is  equally 
presevering. 


College  Station,  Texas 
July,  1968 


R.  0.  Reid 

Professor  of  Oceanography 
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I.  INTRODUCTION 


1 .  Background  and  Scope 

The  problem  considered  in  this  report  deals  with 
the  motions  of  and  tensions  within  a  quasi-elastic  mooring 
line  which  is  anchored  at  the  sea  floor  in  oceanic  depths 
while  attached  to  a  ship  or  buoy  at  or  near  the  sea 
surface  and  subiect  to  the  influence  of  time  varying 
currents.  This  is  a  natural  extension  of  the  many  pre¬ 
vious  studies  dealing  with  the  equilibrium  configuration 
of  an  anchored  cable  in  the  presence  of  steady,  coplanar 
currents.  The  mooring  motion  problem  has  received  much 
less  attention. 

One  of  the  more  recent  investigations  of  the 
equilibrium  problem  is  that  of  Wilson  (1964,  1965)* 
in  which  the  effect  of  vertical  shear  of  the  current  is 
included  (in  the  second  of  the  cited  papers).  A  com¬ 
prehensive  summary  of  previous  studies  related  to  the 
equilibrium  configuration  problem  is  given  by  Wilson  (1964) 
and  will  not  be  repeated  here.  The  results  of  such 
investigations  of  the  equilibrium  configuration  of 
mooring  lines  serve  as  a  useful  base  on  which  one  might 

*  A  list  of  references  is  given  in  alphabetical 
ord^r  by  author  at  the  end  of  this  report. 
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superimpose  a  dynamic  perturbation  theory  related  to 
departures  of  the  current  from  its  mean  state. 

In  the  last  decade  there  has  been  increased  interest 
in  the  general  mooring  motion  problem,  both  in  connect^  on 
with  the  mooring  of  vessels  in  the  deep  sea  and  with  the 
use  of  surface  and  subsurface  anchored  buoys,  which  s^rve 
as  platforms  for  oceanographic  and/or  meteorological 
measurements.  Recent  studies  which  are  particularly 
germane  to  the  mooring  motion  problem  are  those  of 
Whicker  (1958),  Walton  and  Polacheck  (i960),  Polacheck, 
et  al.  (1963),  Paquette  and  Henderson  (1965),  and  Fofonoff 
(1965,  1966).  Although  not  directly  related  to  the  mooring 
problem,  the  paper  of  Zajac  (1957)  on  the  dynamics  of 
laying  and  recovery  of  submarine  cable  should  not  be  over¬ 
looked. 

An  important  aspect  of  the  problem  of  the  transient 
mooring  line,  which  does  not  enter  directly  into  con¬ 
siderations  of  the  equilibrium  configuration,  is  that  of 
the  elastic  and  anelastic  properties  of  the  mooring  line. 
Although  the  qualitative  features  of  the  stress-strain 
processes  in  stranded  cables  and  ropes  are  well  known, 
the  available  information  is  rather  deficient  in  a 
quantitative  sense.  Moreover  an  adequate  theory  of  the 
time-dependent  properties  of  a  rope  or  stranded  cable 
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is  lacking.  Such  a  theory  should  be  able  to  allow  for 
such  phenomena  as  creep  under  nearly  constant  load, 
time -dependent  relaxation  of  elastic  strain  under  no  load 
and  hysteresis  under  conditions  of  cyclic  loading. 

In  a  recent  paper  by  Wilson  (1967)  a  summary  of  some 
of  the  important  properties  of  nylon  and  coir  ropes  and  of 
stranded  steel  cables  is  presented.  This  includes  in¬ 
formation  in  regard  to  the  weight  per  unit  length  in 
water,  the  dynamic  modulus  of  elasticity  (which  is  a 
function  of  the  mean  load)  and  the  permanent  or  an- 
elastic  strain  (which  depends  upon  the  maximum  tension 
that  the  material  has  previously  experienced).  Quantita¬ 
tive  measurements  in  regard  to  hysteresis  and  creep,  as 
well  as  dynamic  modulus  of  elasticity  for  nylon  rope, 
under  cyclic  loading  at  various  mean  tensions,  have  been 
carried  out  recently  by  Paquette  and  Henderson  (1965). 
However,  the  influence  of  hysteresis  was  not  allowed 
for  in  their  transient  state  analysis  of  a  moored  buoy. 

An  attempt  is  made  in  the  present  study  to  include  the 
hysteresis  phenomenon  by  adopting  a  time-dependent  model 
for  the  physical  properties  of  the  mooring  line,  consistent 
with  the  above  information.  This  is  one  of  the  new  aspects 
of  the  present  theory. 
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The  dynamical  model  adopted  in  this  study  allows  for 
three  degrees  of  freedom  of  displacement  of  each  material 
point  of  the  mooring  line  in  response  to  currents,  which 
are  not  necessarily  in  the  same  direction  at  all  depths, 
and/or  in  response  to  three-dimensional  motions  of  the 
floatation  unit  (ship  or  buoy)  to  which  the  mooring  line 
is  attached.  The  mooring  line  is  regarded  as  completely 
flexible  but  extendable  (both  in  the  elastic  and  anelastic 
sense) .  Hysteresis  and  relaxational  phenomena  related  to 
the  time-dependent  stress-strain  processes  are  admitted 
via  a  generalized.  Maxwell  type,  visco-elastic  model  of 
the  physical  behavior  of  the  mooring  line.  This  model  in¬ 
corporates  a  "memory  effect",  reflecting  the  influence 
of  the  prior  history  of  the  stress-strain  process. 

The  damping  of  longitudinal  transient  oscillations 
of  the  mooring  line  is  related  primarily  to  the  hysteresis 
properties  of  the  line.  On  the  other  hand,  the  damping 
of  transverse  oscillations  is  primarily  related  to  the 
mechanical  "radiation"  of  energy  to  the  fluid  through 
the  action  of  hydrodynamic  form  drag,  which  is  dependent 
upon  the  motion  of  the  mooring  line  relative  to  the  fluid. 
These  two  effects  are  considered  to  be  the  dominant 
mechanisms  for  damping  in  the  present  theory.  Specifically, 
the  tangential  skin  drag  on  the  cable  (Reid  and  Wilson, 

1963)  is  regarded  as  negligible  compared  with  the  normal 
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form  drag.  Moreover,  as  in  most  other  studies  of  cable 
dynamics,  thermal  effects  are  not  considered,  nor  is  the 
influence  of  hydrostatic  pressure  on  the  properties  of 
the  cable  (this  being  unknown  at  the  present  time). 

2.  Order  of  Presentation 

Within  the  framework  of  the  above  restrictions,  a 
"general"  theory  of  the  three-dimensional  mooring  line 
dynamics  is  set  forth.  The  presentation  starts  with  a 
formulation  of  the  equations  of  motion  of  the  mooring 
line  in  vector  form  and  a  discussion  of  the  hydrodynamical 
fluid  forces  on  the  line.  This  is  followed  by  a  dis¬ 
cussion  of  the  time-dependent  stress-strain  properties  of 
ropes  and  stranded  cables  and  presentation  of  an 
anelastic  model.  The  implications  of  this  model  in 
respect  to  the  dynamic  modulus  of  elasticity  and 
hysteresis  are  then  compared  with  the  data  of  Paquette 
and  Henderson  (1965)  and  Wilson  (1967). 

A  general  consideration  of  the  energy  conservation 
relation  for  the  system  is  presented.  The  damping  effects 
associated  with  both  form  drag  and  internal  hysteresis 
are  displayed  explicitly  along  with  the  energy  supply 
terms  due  to  work  on  the  mooring  system  at  its  boundaries. 
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In  the  section  which  follows,  the,  component  equations 
of  motion  of  the  mooring  line  are  given  explicitly  in  a 
natural  coordinate  system  similar  to  that  suggested  by 
Zajac  (1957),  which  is  based  on  the  local  orientation  of 
the  mooring  line.  This  system  has  the  advantage  of 
facilitating  the  evaluation  of  the  normal  drag  force  on 
the  mooring  line.  Kinematical  compatibility  conditions 
are  derived  which  relate  the  component  velocities  of  the 
mooring  line  to  other  dependent  parameters  characterizing 
the  cable  orientation  and  strain.  The  system  of  equations 
are  of  quasi-linear,  hyperbolic  type  and  accordingly 
can  be  recast  in  the  so-called  characteristic  form 
(Freeman,  1951)*  The  latter  form  has  the  advantage 
that  a  stable  numerical  solution  of  the  mooring  line 
configuration  at  any  time  t,  starting  with  an  arbitrary 
initial  three-dimensional  configuration  is  greatly 
facilitated.  The  characteristic  form  of  the  basic 
relations  also  clarifies  the  roles  played  by  the  trans¬ 
verse  and  longitudinal  modes  of  motion  of  the  ii  coring 
line  and  hence  adds  to  our  understanding  of  the  mooring 
line  dynamics.  The  numerical  method  of  solution  via 
the  method  of  characteristics  is  discussed. 


ry 
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For  certain  problems  the  linearized  equations  of 
motion  of  the  mooring  line  are  of  advantage.  These 
equations  are  set  forth  in  a  subsequent  section  and 
some  examples  of  their  use  for  small  perturbations  of 
the  mooring  line  from  a  straight  equilibrium  state  are 
discussed. 

The  use  of  the  perturbation  approximation  is  then 
extended  to  the  more  general  case  where  the  equilibrium 
configuration  of  the  mooring  line  is  a  curved  line  in  the 
vertical  plane.  In  this  case  the  longitudinal  and 
transverse  modes  of  oscillation  of  +he  mooring  line 
relative  to  the  equilibrium  state  are  coupled  due  to  the 
curvature.  The  linearized  perturbation  approximation 
is  useful  in  respect  to  considerations  of  the  stochastic 
aspects  of  the  response  of  the  mooring  line  to  irratic 
motion  of  the  floatation  unit  to  which  it  is  attached. 

An  application  of  the  method  is  made  for  the  case  where 
the  motion  of  the  floatation  unit  is  prescribed  in 
terms  of  a  stochastically  stationary  time  sequence 
with  stipulated  variance  spectrum. 


II.  EQUATION  OP  MOTION  OP  MOORING  LINES 


1.  Assumptions 

The  primary  assumptions  in  the  present  analysis 
concerning  the  dynamics  of  a  mooring  line  are  as 
follows : 

a.  The  mooring  line  is  treated  as  perfectly 
flexible,  i.e.,  it  is  supposed  that  the 
line  cannot  support  shearing  stresses 
normal  to  the  line; 

b.  Torsional  strains  in  the  line  are  not 
considered,  it  being  assumed  that  there 
is  no  effect  of  such  strains  on  the 
longitudinal  or  transverse  displacements 
of  the  line; 

c.  Longitudinal  skin  drag  on  the  line  is 
neglected  compared  with  normal  drag 
associated  with  relative  fluid  motion 
past  the  line; 

d.  It  is  considered  that  the  mooring  line  is 
fixed  at  its  lower  end;  i.e.,  it  does  not 
drag  anchor; 
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e.  It  is  supposed  that  the  longitudinal 
strain  in  the  line  is  related  to  the 
axial  tension  and  its  time  history  only, 
thermal  strains  and  exchange  of  heat  with 
the  environmental  fluid  being  ignored; 

f.  It  is  assumed  that  radial  strain  of  the 
line  is  related  directly  to  the  longitudinal 
strain  and  is  independent  of  the  pressure 
exerted  by  the  surrounding  fluid; 

g.  Finally  the  Coriolis  force  on  the  mooring 
line  is  considered  altogether  negligible 
compared  with  other  forces. 

In  the  general  analysis,  no  restriction  is  imposed 
on  the  nature  of  the  current  acting  on  the  mooring  line. 
Accordingly  allowance  is  made  for  three  degrees  of 
freedom  of  displacement  of  any  material  point  of  the 
line.  In  general  the  mooring  line  is  a  three-dimensional 
curve  at  any  instant.  However,  certain  special  cases 
will  be  considered  in  which  the  equilibrium  configuration 
lies  in  a  vertical  plane. 

Assumptions  (e)  and  (f)  are  imposed  largely  in  view 


of  the  limited  state  of  knowledge  in  regard  to  the  physical 
and  thermal  properties  of  typical  mooring  line  materials. 
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The  radial  strain  is  needed  only  in  respect  to  the 
evaluation  of  the  buoyancy  and  normal  drag  force  on 
the  line,  which  depends  on  its  effective  diameter  among 
other  factors. 

2.  Material  Coordinate  and  the  Vector  Equation  of  Motion 
Consider  a  mooring  line  which  is  securely  anchored 
at  point  0  on  the  sea  bed  (Figure  1).  Let  P  be  any 

A 

material  point  on  the  axis  of  the  mooring  line.*  Let  r 
denote  the  position  vector  of  with  respect  to  anchor 
point  0  as  a  reference.  Clearly  the  mass,  M,  of  line 
between  0  and  P  is  conserved  regardless  of  the 
dynamic  state  of  the  line.  One  could  accordingly  adopt 
M  as  a  logical  material  coordinate  of  the  point  P. 
However,  it  is  more  customary  to  employ  an  arc  length 
to  identify  the  point  P.  In  the  present  analysis  we  will 
let  s  denote  the  arc  length  between  the  two  material 
points  0  and  P  of  the  mooring  line  in  its  original, 
completely  relaxed  state,  prior  to  any  permanent  strain. 

We  also  define  m  as  the  mass  per  unit  length  at  s 
in  the  same  state ,  such  that 


*  A  list  of  symbols  appears  at  the  end  of  this 
report  for  convenience  to  the  reader. 
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s 

M  =  J  m(s)  ds,  (1) 

0 

in  which  m  (for  a  tapered  line)  may  in  general  depend 
on  s.  In  this  sense  s  is  merely  an  alias  of  M. 

Note  that  s  is  not  to  be  confused  with  the  actual 
arc  length,  a,  of  the  line  between  0  and  P  in  its 
general  strained  state.  In  general  a  >  s.  Moreover 
the  mass  per  unit  actual  length  under  strained  conditions 
will  generally  be  less  than  m,  and  is  equal  to 
m  ds/da. 

All  dependent  variables  such  as  r,  a,  etc.  will  be 

regarded  as  functions  of  the  independent  variables  s 

and  t.  To  avoid  confusion  with  a  we  will  hereafter 

refer  to  s  as  the  material  coordinate. 

The  vector  velocity  of  any  material  point  of 

* 

the  line  will  be  denoted  by  V  and  is  related  to 

a 

r(s,t)  by 


V  =  dr/  dt.  (2) 

The  translational  momentum  of  the  line  between  0  and  P 

is  accordingly  given  by 

s 

J  mV  ds. 
o 
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A  A 

Let  T  denote  the  vector  tension  at  P  and  TQ 
the  vector  tension  at  0,  both  taken  in  the  sense  cf 
increasing  s  along  the  line.  These  tensile  loads 
are,  in  viev;  of  assumption  (a),  tangent  to  the  axis  of 

A 

the  line  at  the  points  in  question.  Let  k  denote  the 
vertical  unit  vector  (upwards)  and  g  the  acceleration 
due  to  gravity  such  that  in  the  absence  of  fluid  the 

A 

line  between  0  and  P  has  the  weight  -  M  g  k  or 

s 

f*  A 

-  I  m  g  k  ds. 
o 

A 

Finally  let  F  denote  the  residual  surface  force  acting 
on  the  line  between  0  and  P  due  to  the  action  of 
the  surrounding  fluid. 

Application  of  Newton's  law,  bearing  in  mind 
assumptions  (a)  and  (g),  then  leads  to  the  relation 

s  s 

a/dt  |  m  V  ds  =  T  -  Tq  -  J  m  g  k  ds  +  F.  ( 3) 

A  A 

Now  m  is  independent  of  t  and  Tq  =  T  (o, t), 
so  it  follows  from  (3)  by  differentiating  with  respect 
to  s  that 

A  4^ 

m  av/at  =  aT/  as  -  m  g  k  +  f,  (4) 


* 

in  which  T  has  the  direction  of  the  vector  dr/  os 
*  * 

and  f  =  dF/ds. 

3.  Hydrodynamic  Forces  Due  to  the  Fluid 

In  the  absence  of  any  motion  of  the  fluid  or  of 

a 

the  mooring  line,  the  force  F  reduces  to  the  static 

A 

Archemedian  buoyance  force,  B,  given  by 

a 

B  =  gk  j  pA  do  (5) 

o 

.  which  p  is  the  fluid  density  (mass  per  unit  volume) 
and  A  is  the  cross-sectional  area  of  the  material  of 
the  line  at  s  under  strained  conditions.  The  buc>yancy 
per  unit  material  coordinate  is  accordingly  given  by 

f  =  dR/ds  =  n  cr  k  A  drr/ds  *  (6) 

For  simplicity  it  will  be  supposed  that  p  is  a 
constant.  However  A  and  do/ ds  can  in  general 
•  vary  with  s  or  t. 

In  the  presence  of  motion  of  the  fluid  relative  to 


the  mooring  line,  an  additional  dynamic  force  can  occur. 
This  is  associated  with  anomalies  of  pressure  distribution 
around  the  surface  of  the  immersed  line.  Such  pressure 
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anomalies  can  be  cieated  by  acceleration  of  the  fluid 
and/or  by  acceleration  of  the  line  through  the  fluid. 

In  addition,  at  moderate  to  large  Reynolds  number,  a 
wake  phenomenon  can  exist  on  the  lee  side  of  the  line, 
which  leads  to  a  form  drag  effect  above  and  beyond  the 
influence  of  acceleration.  Also  it  is  well  known  that 
periodic  vortex  shedding,  which  occurs  at  moderate 
Reynolds  number,  can  lead  to  an  alternating  lateral  force 
on  the  line  (normal  to  the  relative  current)  at  a  fre¬ 
quency  deper.dent  upon  the  the  Reynolds  number  (Strouhal 
frequency) . 

A 

The  drag  force,  fd,  per  unit  material  coordinate 
is  presumed  to  be  of  the  form 

fd  =  (V2)P  D  |Un  -  Vn|  (U„  -  V/  (?) 


where  D  is  the  effective  diameter  of  the  line  at 
position  s  and  time  t,  is  the  drag  coefficient 

which  is  a  function  of  the  Reynolds  number  associated 
with  the  relative  speed,  is  the  component  of  V 

normal  to  the  mooring  line  at  s,t  and  U  is  the 

A 

component  of  the  fluid  velocity  U  normal  to  the  mooring 
line  at  s,  t*.  The  above  relation  gives  a  vector  force 


*  Actually  U  is  the  fluid  velocity  that  would 
exist  in  the  absence  of  the  mooring  line. 
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normal  to  the  line  at  s,  t,  the  tangential  skin  drag 
being  neglected.  Values  of  c^  for  normal  drag  on 
stranded  cable,  as  a  function  of  the  Reynolds  number, 
are  summarized  by  Wilson  (1964).  The  tangential  drag 
coefficient  is  of  the  order  of  only  one  per  cent  of  that 
for  normal  drag  for  usual  ranges  of  Reynolds  number 
(see  Wilson,  1964  or  Reid  and  Wilson,  1963).  Accordingly 
the  omission  of  tangential  skin  drag  is  not  considered 

A 

serious.  The  vector  Un  can  be  expressed  in  the  form 


A 


(8) 


where  t  is  a  unit  vector  tangent  to  the  axis  of  the 

A  A 

mooring  line  at  s,  t  and  UT  =  t  •  U.  A  similar 

A  A 

relation  of  course  holds  for  V  .  The  unit  vector  t 

A 

is  related  to  r  by  the  relation 


T 


dr  ydo 

Ssfk  ’ 


(9) 


The  effect  of  vortex  shedding  can  be  simulated  by 

A 

allowing  an  oscillatory  component  of  f^  normal  to  the 

a  A 

direction  of  f,  and  the  local  vector  t.  For  the 
d 

evaluation  of  the  SLrouhai  frequency  at  given  Reynolds 
number,  see  for  example  Schlichting  (i960).  The  amplitude 
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of  the  resulting  transverse,  alternating  force  is 
presumably  of  tihe  order  of  magnitude  of  ten  per  cent 

A 

of  the  local  mean  value  of  fd.  If  we  let  ti  denote  the 
ratio  of  the  amplitude  of  lateral  force  to  the  drag  Jforce 
and  let  ios  be  the  Strouhal  frequency,  then  the  alter¬ 
nating  lateral  thrust  per  unit  length  due  to  vortex 
shedding  can  be  approximated  by 

*  A  A 

f g  =  T]  T  X  fd  COS  (u>3  t  +  6)  ,  (10) 

where  6  is  an  arbitrary  phase  angle. 

To  the  above  drag  force,  which  can  exist  under 
steady  relative  motion,  we  must  allow  for  the  possibility 
of  an  added  inertial  force.  Consider  for  the  time  being 

A 

that  the  fluid  motion  U  (away  from  the  influence  of  the 
mooring  line)  is  steady,  but  that  the  mooring  line  is 
accelerated  in  a  direction  normal  to  the  line  like  a 
rigid  cylinder.  In  this  case  a  portion  t>£  the  fluid 
in  the  neighborhood  of  the  line  must  also  be  accelerated, 
thus  leading  to  sun  effective  retarding  force  on  the  line. 
The  magnitude  of  this  force  per  unit  length  of  line  can 
be  expressed  as  -  m„  a„  where  a  is  the  acceleration 
normal  to  the  line  and  mQ  is  the  effective  added  mass 
of  fluid  being  accelerated  at  the  rate  an.  If  the  line 


.8 


experiences  an  acceleration  tangential  to  the  line  and  if 
viscous  effects  are  ignored  then  there  will  be  essentially 
no  acceleration  of  the  fluid  if  end  effects  are  also 
ignored  (very  long  line).  Thus  we  should  expect  the 
inertial  force  to  be  highly  directional,  like  the  drag 
force. 

A 

For  a  general  acceleration  of  the  line  but  with  U 
steady,  we  propose  the  following  expression  for  the 
inertial  force  per  unit  length  due  to  the  presence  of 
the  fluid 


f *  =  -  m  (dV/dt  -  t  dVT/dt)  ,  (11) 

i  a 


where  VT  =  V  •  t,  this  being  the  tangential  component 
of  the  mooring  line  velocity  at  s,  t.  The  above  relation 

A 

is  not  the  only  possible  choice  for  f%  Other  logical 
choices  are 


-  m  [  dV/8t  -  t  (dV/dt)  *  t  ] 

cL 


and 

-  m  a/at  (V  -  t  VT)  . 

ct 

A 

If  t  were  independent  of  time  then  each  of  the  three 

A 

possible  forms  yield  identical  values  of  f '  . 

i 


Thus 
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the  differences  result  only  from  terms  involving  br/bt, 
which  represents  a  turning  rate  associated  with  the 
local  orientation  of  the  axis  of  the  mooring  line.  Of 

A  A 

the  above  three  possible  forms  for  f’  (for  steady  U), 
only  that  given  by  (11)  leads  to  a  physically  acceptable 
energy  equation  for  the  system  as  we  will  see  in  a  later 
section.  The  choice  of  relation  (11)  is_  justified  pri¬ 
marily  on  this  basis. 

In  the  more  general  case  where  the  velocity  of  the 
fluid,  U  (in  the  absence  of  the  mooring  line)  is  not 

A 

steady,  then  the  acceleration  term  dU/dt  will  contribute 
to  the  inertial  force.  Its  effect,  above  and  beyond 

A 

that  of  bV/bt,  is  analysed  as  follows.  In  the  absence 
of  the  mooring  line,  the  fluid  occupying  the  position 
of  the  mooring  line  would  experience  a  force  pA  ba/bs 

A 

dU/Bt,  per  unit  s,  where  A  da/ds  represents  the 
volume  per  unit  s.  This  force  exists  by  virtue  of  the 
pressure  gradient  required  to  produce  the  acceleration 

A 

dU/dt.  In  the  presence  of  the  mooring  line,  a  disruption 
of  the  flow  occurs  in  its  vicinity  and  this  incrouaces 
an  added  inertial  force  similar  to  that  caused  by 
accelerating  the  mooring  line  through  the  fluid.  The 
added  inertial  force  is  taken  as 
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m  (dU/dt  —  t  &U  /dt)  , 
a  T 

by  analogy  with  (11).  Thus  in  the  absence  of  any 
acceleration  of  the  mooring  line  the  inertial  force 
per  unit  of  s  is  considered  to  be  given  by 

* 

f^  =  p  A  ba/bs  bU/bt 

+  m  (dU/dt  -  t  dUT/  bt)  .  (12) 

a 

The  total  inertial  force  per  unit  of  s  exerted 
by  the  fluid,  in  the  presence  of  acceleration  of  the 
fluid  as  well  as  acceleration  of  the  mooring  line,  is 
accordingly  taken  as  the  sum  of  the  right  hand 
sides  of  (11)  and  (12). 

The  added  mass  coefficient  m  per  unit  s  can 

St 

be  expressed  in  the  form 

m  =  p  A  c  bo/bs  (13) 

St  3. 

where  c  is  a  non-dimensional  added  mass  coefficient. 

St 

The  theoretical  value  of  c&  for  pure  potential  flow 
around  a  circular  cylinder  according  to  Lamb  (19^5) 


is  exactly  unit.  This  value  applies  for  very  small 
Reynolds  number.  For  large  Reynold  number,  the  value 
of  c  is  generally  somewhat  less  than  unity. 

CL 


4.  Alternative  Form  of  the  Equation  of  Motion 

The  total  force  per  unit  material  coordinate  due 
to  the  fluid  is  given  by 


A  A  A 


f  =  f.  -f  f  „  +  f  4-  f  *  +  f" 

b  d  s  i  i 


(14) 


Making  use  of  (6),  (7),  (10),  (11),  (12),  and  (14)  in 
(4)  and  combining  terms  involving  the  acceleration 
of  the  mooring  line  yields 


/\  .  "V -  /N 


(m  -t-  iua)  oV/ot  -  m  t  0VT/dt  =  dx/ds  -  wk  +  K,  (lb) 

£L 


where  w  is  the  net  weight  per  unit  s  given  by 


w  (m  -  md)  g 


(16) 


in  which  m^  is  the  displaced  mass  of  fluid  per  unit  s 


m^  =  p  A  da/bs. 


(17) 
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The  residual  fluid  force  R  is  given  by 

A  ^ 

R  =  Z  cd  ®|f  I Rn  -  vnl  [1  +  +  6)  «]  (Un  -  Vn) 

A 

+  (m  +  m  )  4?  -  m  r  .  (18) 

'  d  a'  dt  a  IFtT 


It  will  be  recalled  that 


m  =  c  in 
a  ad 

(19) 

and 

^  A  A 

Un  =  u  -  T  uT 

A  A  A 

V  =  V  -  T  V 

n  t 

(20) 

A  ^  A  A 

where  UT  =  U  •  t  and  Vt  =  V  *  t.  Moreover, 

A 

and  da/d s  are  related  to  r  as  follows 

A  A 

V,  T 

A  ^ 

Y  =  dr/dt 

(21) 

A 

a  dr  So 

T  -  'sr/  3? 

(22) 

da  _  I  dr. 
3?  "  '3s  I  * 


(23) 


23 


The  vector  tensile  force  can  be  expressed  as 

A  A 

T  =  T  t  (24) 

where  T  is  the  scalar  tension  at  s,  t  and  is  related 
to  the  longitudinal  strain  (da/ds  -  1)  by  un  equation 
of  state  appropriate  to  the  particular  mooring  line 
concerned.  This  will  be  considered  in  some  detail  in  a 
later  section. 

A  A  A  A 

The  dependent  variables  r,  V,  t,  T  and  a  are 
sought  as  functions  of  s  and  t  for  given  initial 

A 

and  end  conditions.  The  velocity  U  is  understood  to 

* 

be  a  prescribed  function  of  r  and  t.  Its  dependence 
on  s  is  not  given  explicitly  but  is  implied  through 
the  dependence  of  r  on  s. 

5.  End  Conditions  -  Vessel  Motion 

At  the  anchor  point  for  each  mooring  line  we  require 

simply 

r  (o,  t)  =  0,  (25) 

A 

which  implies  in  turn  that  V  (o,  t)  also  vanishes. 

We  will  restrict  our  analysis  to  mooring  lines 
without  any  intermediate  attachments  between  the  anchor 
and  the  surface  vessel  to  which  each  is  secured.  Let 
L  denote  the  material  coordinate  at  the  point  of 
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attachment  of  the  given  mooring  line  to  the  vessel  (i.e., 

L  is  the  original  unstrained  total  length  of  the  mooring 
line).  A  possible  simple  end  condition  at  s  =  L  is 

A  * 

to  specify  V  (L,  t)  as  a  prescribed  function  of  time. 

a 

This  together  with  specified  initial  values  of  r  and 

a 

V  versus  s  is  admissible  from  a  mathematical  point  of 
view.  However,  from  a  physical  point  of  view,  it  does 
not  seem  rational  to  specify  the  motion  of  the  moored 
vessel  a  priori.  Clearly  the  action  of  the  mooring  line 
plays  a  role  in  constraining  the  vessel's  motion  and 
indeed  can  introduce-  natural  frequencies  of  oscillation 
in  addition  to  those  which  are  relevant  to  an  unmoored 
vessel.  The  unmoored  vessel,  of  course,  possesses 
distinct  natural  frequencies  for  roll,  pitch,  and  heave 
only.  The  constraint  imposed  by  the  mooring  line  (or 
lines)  can  lead  to  additional  natural  frequencies  in 
surge,  sway,  and  yaw  (O'Brien  and  Muga,  1963). 

Thus  it  would  be  desirable  to  deduce  the  motion 
of  the  moored  vessel  itself  in  terms  of  the  hydrodynamic 
forces  acting  on  it,  together  with  the  constraining  effect 
of  the  mooring  line.  This  implies  that  r  (L,  t)  anu 
T  (L,  t)  for  the  mooring  line  (or  lines)  are  related 
by  an  appropriate  differential  equation  in  which  the 
wave  action  and/or  steady  current  acting  on  the  vessel 
should  enter  as  the  prescribed  exciting  force.  The 
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approach  employed  here  is  a  compromise  in  which  the  motion 
of  the  vessel  in  the  absence  of  the  mooring  line  (the 
potential  motion)  is  prescribed  rather  than  the  exciting 
forces.  This  of  course  requires  the  specification  of  six 
scalar  functions  of  time  corresponding  to  the  six  degrees 
of  freedom  of  motion  for  a  rigid  vessel.  Typical  measure¬ 
ments  of  roll,  pitch,  and  heave  for  large  vessels  underway 
in  a  seaway  are  given,  for  example  by  Gover  (1955)* 
Measurements  of  all  six  motions  for  a  large  moored  vessel 
are  presented  by  O'Brien  and  Muga  (1963). 

Approximate  relations  connecting  the  motion  of  a 
moored  vessel  with  the  prescribed  potential  motion  for 
given  mooring  conditions  are  derived  in'Appendix  A 
based  upon  the  linearized  equations  of  motion  of  the  type 
employed  by  St.  Denis  and  Pierson  (1955),  but  allowing 
for  the  constraint  imposed  by  the  mooring  line  (or  lines). 
Allowance  is  made  for  torque  exerted  by  the  mooring  lines 
as  well  as  the  constraining  force  since  in  general  the 
attachment  of  a  given  mooring  line  is  not  located  at  the 
center  of  the  mass  of  the  vessel.  The  resulting  equations 
are  valid  for  small  translational  and  rotational  displace¬ 
ments  of  the  vessel.  These  relations  are  summarized 
below. 

a 

Let  r^  denote  the  position  vector  of  the  point 
of  attachment  on  the  vessel  of  mooring  line  k.  relative 
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to  the  anchor  point  of  mooring  line  k.  Let  pk  denote 
the  separation  vector  from  the  vessel's  center  of  mass  (C) 
to  the  attachment  point  of  mooring  line  k  (Fig.  2). 

A 

The  vector  p^  has  constant  magnitude  but  varies  in 
direction  depending  on  the  angular  displacement  of  the 

A  A 

vessel.  The  mean  values  of  the  vectors  r^  and  p^ 

A  A 

are  denoted  by  rk°  and  pk°  respectively.  Now  let 

A  A 

X  and  a  denote  the  potential  translational  and 

rotational  displacements*  of  the  vessel,  in  a  given 

seaway,  if  it  were  not  moored  and  not  underway.  These 

two  vectors  are  regarded  as  prescribed  functions  of 

time.  For  small  rotational  displacements,  the  position 
* 

vectors  r^  at  the  attachment  points  can  be  approximated 

by 


A  A  A  A  ^ 

r  =  r  °  +  X  +  a  x  p  0 

k  k  k 

+  R'  +  *'  x  pk°  (26) 

A  A 

where  R 1  and  if  ’  represent  anomalies  of  the  translational 
and  rotational  displacements  of  the  vessel  associated  with 

A  A 

the  constraint  of  the  mooring  lines.  Thus  X  +  R'  is  the 
vector  trans la tional  displacement  of  the  center  of  mass 


*  Referred  to  equilibrium  state. 
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of  the  moored  vessel  relative  to  its  equilibrium  position. 

A  * 

Likewise  a  +  f '  is  the  vector  rotational  displacement 
of  the  moored  vessel  relative  to  its  equilibrium 
orientation. 

Let  fij  ',  Rg  * ,  R3  '  denote  respectively  the  surge, 

a 

sway,  and  heave  components  of  the  vector  R’  and  let 
♦x  «,  ♦a  *,  f3 1  denote  respectively  the  components  of 

A 

♦  ’  in  roll,  pitch,  and  yaw.  These  components  represent 
a  Cartesian  system,  the  reference  axes  of  which  are  taken 
as  the  equilibrium  orientation  of  the  principal  axes  of 
the  vessel.  As  shown  in  Appendix  A,  the  approximate 
equations  of  motion  governing  the  above  quantities  are 
(in  compact  form)  - 

(d2/dta  +  2$.  d/dt  +  a.s)  Rj»  =  (27) 

k 

and 

(da/dt3  +  20  *  d/dt  +  TS  )♦  J  *  =  ’I  J,kjAj  (28) 

°  k 

where  j  =  1,  2,  or  3.  The  sums  on  the  right  are  taken 
over  the  total  number  of  mooring  lines.  The  terms 
T1  represent  the  Cartesian  components*  of  the  vector 

kj 

*In  the  same  system  as  that  for  R  - 1  and  $  .  • . 

<3  J 
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T^ ’  where 


(29) 


Tk  being  the  vector  tensile  force  at  the  point  of  attach- 

A 

ment  on  the  vessel  for  line  k  (Fig.  2)  and  Tk° 

A 

is  the  mean  value  of  T  (or  equilibrium  value).  Like- 

K. 

wise  the  terms  J*  .  are  the  Cartesian  components,  in 

kj 

a 

the  same  system  of  the  vector  J^.'  defined  by 


v 


(30) 


where 


Jk  »  (Pk  +  (*'  +  a)  *  Pk  ]  x  Tk  (31) 

and  Jy°  is  the  mean  value  of  J^.  Physically  -Jk 

represents  the  torque  exerted  by  line  k  on  the  vessel. 

The  terms  M .  are  the  effective  masses  of  the 
J 

vessel  (including  effective  added  mass  of  water)  for 

motion  in  surge,  sway,  and  heave  (j  =  1,  2,  3  respectively). 

The  terms  I .  are  the  effective  moments  of  inertia  about 
J 

the  principal  axes  for  roll,  pitch,  and  yaw  (j  -  1,  2,  3 
respectively) .  The  coefficients  p  and  p  . 1  are 

j  J  , 
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are  damping  rates  (units  of  reciprocal  time)  for  the  six 

components  of  motion.  Finally  o  and  a . 1  are  the 

d  J 

nominal  natural  frequencies  (radians  per  unit  time)  for 
the  unmoored  vessel.  Note  that  ax  =  a2  =  a3  *  =  0 
(for  surge,  sway,  and  yaw). 

a  ^ 

Equations  (2?)  and  (28)  imply  that  R*  and  ij» ' 

* 

are  appropriate  integrals  of  the  dependent  vector 

A 

and  specified  function  a.  Thus  (26)  implies  a  relation 
between  the  dependent  variable  r^  at  the  upper  end  of 

a 

line  k  and  the  T  for  all  lines.  Moreover,  there  are 

Xv 

as  many  relations  (26)  as  there  are  lines,  which  implies 

that  the  system  should  be  determinate .  As  a  special  case, 

* 

if  R 1  and  $ '  are  negligible  compared  with  X  and 
* 

a  respectively,  then  (26)  reduces  to  the  simple  relation 

A  A  A  *  A 

rj.  =  r  0  +  X  +  a  x  p  o  (26a) 

K  k  k 

A  A 

in  which  X  and  a  are  prescribed. 

Whether  one  chooses  to  specify  X  and  a  or 

A  A  A  A 

(X  +  R')  and  (a  +$’)  is  perhaps  optional;  in  either 
case  the  motion  of  the  upper  end  of  the  mooring  lines  (in 
the  case  of  multiple  mooring)  should  be  consistent  with 
(26). 
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III.  STRESS-STRAIN  RELATIONS  FOR 
ROPES  AND  STRANDED  CABLES 


1.  Known  Properties 

Some  examples  of  load  versus  elongation  for  ropes 
and  stranded  cables  are  shown  in  Figs.  3  and  4.  In 
Fig.  3>  which  is  taken  from  Wilson  (1967),  typical 
processes  for  steel  wire  cables  (upper  panels)  and  for 
coir  mooring  rope  (lower  panels)  are  presented.  Fig.  4, 
which.  io  taken  from  Paquette  and  Henderson  (1965), 
represents  a  test  carried  out  for  a  one-half  inch 
diameter  nylon  rope  of  55- inch  original  length.  In  the 
latter  test  the  tension  was  cycled  over  ranges  of  +  180 
lbs.  relative  to  ten  selected  mean  tensions.  Such 
tests  demonstrate  certain  features  common  to  both 
stranded  cables  and  ropes  which  are  summarized  below: 

a.  The  stress-strain  relation  is  definitely 
process  dependent;  i.e.,  the  relation  between  stress 
and  stain  at  time  t  is  dependent  upon  the  prior 
stress-strain  conditions; 

b.  The  material  suffers  a  permanent  anelastic 
strain  which  is  dependent  upon  the  maximum  load 
previously  experienced; 
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c.  For  cyclic  loading  and  unloading  a  hysteresis 
loop  exists  which  implies  that  the  stored  energy  is  not 
entirely  elastic; 

d.  Anelastic  creep  is  evident  for  cyclic  loading 
particularly  for  large  mean  loads;  the  creep  tends  to 
reach  a  limit  after  many  cycles  so  as  to  form  a  closed 
hysteresis  loop; 

e.  The  amount  of  hysteresis  energy  loss  per  cycle 
definitely  depends  upon  the  range  of  cyclic  load; 

f.  The  material  can  experience  a  time  dependent 
relaxation  of  strain  (negative  creep)  upon  sudden 
release  of  load; 

g.  Under  cyclic  loading,  the  apparent  spring 
coefficient  (slope  of  the  principal  axis  of  the  hysteresis 
loop  for  small  range  of  load)  increases  with  increasing 
mean  load  (Fig.  4) . 

The  latter  effect  can  also  be  seen  in  terms  of  the 
curved  axis  of  the  large  hysteresis  loops  shown  in  Fig.  ? 
(dashed  curves  in  upper  panels).  The  rope  or  cable 
acts  like  a  stiffening  spring  as  well  as  displaying 
significant  hysteresis. 

Conclusive  measurements  are  lacking  in  regard  to 
the  influence  of  cycling  rate  on  the  apparent  spring 
coefficient  and  hysteresis  for  a  given  mean  load. 


PTUmmm  CITfUPO^  M 


Repeated  load-elongation  tests  for  mooring 
line  materials:  (a)  and  (b)  steel-wire 
ropes;  (c)  and  (d)  coir  fiber  ropes  (after 
Wilson,  1967). 
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2.  A  Simple  Maxwell  Model 

A  possible  approach  for  representing  the  stress- 
strain  processes  in  a  rope  or  cable  is  to  adopt  a 
simple  three-parameter  Maxwell  model  to  simulate  its 
visco-elastic  properties  (Freiberger,  i960,  p.  579  and 
77^).  Such  a  model  is  shown  schematically  in  Fig.  5. 

In  this  model  Kx  and  K,,  represent  spring  coefficients 
with  units  of  force  per  unit  strain  and  N  is  a 
frictional  coefficient  with  units  of  force  per  unit  time 
rate  of  strain.  Let  e  be  the  overall  strain  for  the 
system;  this  is  the  same  as  the  strain  in  segment  DE 
or  in  segment  AC  (Fig.  5).  Let  ex  and  £  be  the 
partial  strain  in  segments  between  AB  and  BC, 
respectively.  Let  Tx  be  the  common  tension  in  the 
series  segments  AB  and  BC  and  let  T0  be  the  tension 
in  segment  DE.  The  total  load,  T,  is  given  by 
T  =  T0  +  Tx  ;  moreover  e  =  ex  +  §.  Finally  we  take 
Tx  =  Kj  ex  =  N  d§/dt  and  T0  =  Kg  e .  Eliminating 
T0  ,  Tx ,  ex  and  e2  between  the  above  five  relations 
leads  to  the  following  constitutive  equation  for  the 
system 


=  Kq  e  +  (Kg  +  iq  )  t  de/dt. 


T  +  t  dT/dt 


(32) 
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where  t  =  N /K*  .  The  latter  parameter  represents  a 
characteristic  relaxation  time  for  the  system.  In 
the  above  relation  it  is  presumed  that  the  three 
parameters  K0 ,  ,  and  t  are  constants.  A  relation 

essentially  equivalent  to  (32)  was  deduced  by  Meixner 
(1965)  based  upon  an  entirely  different  approach  in¬ 
volving  the  thermodynamic  non-equilibrium  processes 
in  anelastic  materials. 

The  above  model implies  that  fcr  very  slow  changes 
(32)  reduces  to 

T  *  Ko  «,  (33) 

while  for  rapid  changes 

at  =  (Ko  +  iq  )  As.  (34) 

Thus  the  system  is  characterized  by  two  limiting  elastic 
stress-strain  relations  which  we  will  refer  to  as  static 
and  dynamic. 

Consider  the  following  special  time  dependent 
process  illustrated  in  Fig.  6: 

a.  Starting  at  zero  load  and  zero  strain  apply 
a  sudden  load  T‘A  resulting  in  elongation  eA  =  TA/(Ko  +  ) . 
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b.  Under  the  sustained  load  T^,  (32)  predicts 

that  e  will  creep  ultimately  to  the  value  =  T^/Kq 
after  an  elapsed  time  of  many  t  units. 

c.  Having  reached  point  B  now  suppose  the 
load  is  suddenly  removed;  the  strain  immediately  after 
removal  should  be  ec  =  eB  “  eA  s^-nce  curve  BC  is 
parallel  to  curve  OA. 

d.  Again  at  constant  (zero)  load,  the  system 

is  not  in  equilibrium  and  (32)  predicts  a  relaxation  of 
the  strain  with  nearly  complete  recovery  after  an  elapsed 
time  of  many  t  units. 

The  above  process  is  an  idealized  hysteresis  loop; 
in  fact  it  gives  the  maximum  hysteresis  compared  with 
any  ether  cycle  having  the  same  range  of  T.  For 
simple  harmonic  cycling  of  T,  the  hysteresis  predicted 
by  this  model  will  be  frequency  dependent. 

Suppose  T  is  cycled  according  to  the  relation 

T  =  T0  +  a™  cos  u)t  ,  (35) 

then  a  periodic  solution  for  s  satisfying  (32)  is 
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T 

e  =  ^  +  [Ko  +  (Kq  +  iq  )  (udt)2  ]  cos  cut 

+  Kj  cut  sin  cut  [V  +  +  Kj  )2  (wt)*]-i  ,  (36) 

where  both  T  and  e  are  of  period  2ir/cu.  Since  e  is 
out  of  phase  with  T,  a  plot  of  T  vs  e  as  in  Pig*  7 
shows  a  characteristic  hysteresis  loop  Had  we  allowed 

for  an  arbitrary  initial  condition  then  the  resulting 
plot  of  T  vs  e  would  also  show  some  creep  effect 
as  in  Fig.  4,  with  T  vs  e  ultimately  approaching 
the  stable  hysteresis Joop  implied  by  (35)  and  (36). 

The  hysteresis  for  a  closed  cycle  will  be  defined 

by 


H  * 


T  de 


(37) 


For  the  idealized  process  represented  schematically 
in  Fig.  6,  the  hysteresis  is  simply 

Hi  =  TA  (eB  "  £a) 
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or 


H. 

l 


-  K,,  K  +  %) 


(38) 


where  AT  =  T./2.  For  simple  harmonic  cycling  we  can 
compute  H  by  rewriting  (27)  in  the  form 


H 


dt  . 


Using  (35)  and  (36)  this  gives 


(39) 


TT  Kl  { AT)8  mot _ 

(V  +  (Ko  +  K,)*  (u>T)=a 


(40) 


A  non-dimensional  plot  of  H  versus  ou  according  to 
this  relation  is  shown  in  Fig.  8.  The  latter  has  a 
maximum  value  when  idt  =  KQ  /(Kg  +  ) ,  the  maximum 

value  of  H  being 

Hm  =  2  ^  [CV  iq )  * 


(41) 


^3 


Comparison  with  (38)  indicates  that  the  maximum 
hysteresis  for  sinrnle  harmonic  cycling  is  tt/8  times 
that  for  the  idea  :ea  process. 

The  slope  of  the  major  diagonal  of  the  hysteresis 
ellipse  (line  AB,  Fig.  7)  represents  an  apparent  spring 
coefficient  K*  whose  value  can  be  shown  to  be  given  by 


f  V  +  (Kp  4  Kx  )a  (orrja  '  & 

1  1  +  (iut)2  J 


(42) 


For  out  «  1,  K*  is  nearly  equal  to  Kq  while  for 
•ut  »  1,  K*  approaches  K0  +  ^  .  Thus  the  apparent 
spring  coefficient  for  this  model  is  frequency  dependent. 

The  elastic  energy  per  unit  length,  Ee>  for  the 
Maxwell  model  is  the  sum  of  that  stored  in  segments  AB 
and  DE  (see  Fig.  5)-  This  can  be  expressed  in  the 
form 

Ee  =  I  +?Ki  d2  •  (43) 

The  elastic  strain  ex  can  be  rewritten  as  the  difference 
e  -  §  where  §,  it  will  be  recalled,  is  the  anelastic 
strain  for  segment  BC.  Thus  the  elastic  energy  can  be 
rewritten  as 
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Ee  =  |  *0  |  (e  -  S)2  •  W 

Now  consider  the  partial  derivatives  of  the  function 
Ee  (e,  5): 

dE  Vde  =  (Ko  +1^)  c  -Kj  5  , 

(*5) 

=  -\e  +  Kj?  . 

The  above  two  relations  can  be  shown  to  be  consistent 
with  (32)  if  and  only  if 

dE  /de  *  T, 

(46) 

dE  =  -N  d§/3t. 


where  N  =  KjT.  These  relations  were  employed  as  a 

starting  point  in  the  analysis  by  Meixner  (1965).  They 

are  asserted  to  be  more  general  than,  the  special 

Maxwell  relation  since  E  need  not  have  the  form  (44) 

e 

for  more  general  materials.  Indeed  relations  (46)  will 
serve  as  .a  useful  point  of  departure  in  the  considerations 
of  the  following  section. 
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3.  Generalization  of  the  Maxwell  Model 

Although  the  above  simple  model  allows  for 
hysteresis  and  associated  relaxational  phenomena,  it 
is  obviously  deficient  in  at  least  two  respects.  First, 
there  is  no  allowance  for  permanent  elongation.  Second, 
it  assumes  constant  values  of  elastic  coefficients. 

Both  of  these  shortcomings  can  be  overcome  by  adopting 
a  scheme  similar  to  that  suggested  by  Wilson  (19^7) . 

In  essence  this  scheme  distinguishes  between  the  quasi¬ 
elastic  strain  e  -  eQ  and  the  permanent  strain  e0 , 
in  which  the  latter  is  a  function  of  the  maximum  tension 
which  the  material  has  previously  experienced.  Morec 
for  given  eo,  the  tension  is  regarded  as  a  non-linear 
function  of  the  quasi-elastic  strain.  The  data  further¬ 
more  indicate  that  the  form  of  this  dependence  is 
independent  of  e0. 

The  proposed  generalized  stress-strain  model  which 
incorporates  these  features  can  be  obtained  in  a  sys¬ 
tematic  manner  by  starting  with  a  rational  generaliza¬ 
tion  of  the  elastic  energy  function  Eg.  It  vail  be 
understood  that  Eg  is  the  elastic  energy  per  unit 
material  coordinate  (s).  As  before  we  will  regard  this 
as  a  function  of  the  total  strain  e  and  the  anelastic 
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strain.  The  latter,  however,  consists  of  a  passive 
permanent  strain  e0  and  an  active  anelastic  strain  §. 
These  are  distinguished  by  the  fact, that,  under  normal 
conditions,  d§/dt  can  be  positive  or  negative  while 
de0/dt  cannot  be  negative  unless  the  material  is 
thermally  or  mechanically  reworked.*  Specifically  e0 
represents  the  equilibrium  value  of  e  at  zero  load. 

The  elastic  part  of  the  strain  is  (e  -  e0  -5)  and 
the  counterpart  of  e  in  the  previous  model  is  now 
(e  -  eQ).  With,  this  in  mind  the  proposed  generalization 
for  Ee  is  the  following: 

Ee  =  f(c  -  e0)  +  (iq/2)  (e  -  e0  -  ?)2  ,  (47) 

in  which  is  a  constant  with  dimensions  of  force 

and  f(a)  denoces  a  non-linear,  positive  function  of 
a  e  -  e0 ,  dependent  upon  the  material  and  con¬ 
figuration.  Relation  (44)  is  a  special  case  in  which 
%=  0  and  f(e)  =  (Xr)/2)e9 .  Using  (46)  yields 


*In  essence,  this  is  a  manifestation  of  the  second 
law  of  thermodynamics  as  applied  to  real  materials  in 
which  e0  is  related  closely  to  the  entropy  of  the 
material,  at  least  for  adiabatic  conditions. 
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T  =  f  «(e  -  e0)  +  \  (e  -  e0  -  |)  (48) 


and 


Ti  =  (e  -  e0  -  §)  ,  (49) 

where  f*(a)  s  df/dct.  and  tx  =  N/i^  .  The  latter 
time  constant  now  carries  a  subscript  since  we  will 
introduce  at  least  one  more  time  constant  later. 

For  very  slow  rate  of  strain,  (49)  reduces  to  the 
equilibrium  relation 

e  -  ec  -  §  =  0  (50) 

and  the  associated  quasi-static  load-strain  relation 
becomes 


T  =  f’(e  -  e0).  (51) 

On  the  other  hand,  for  very  rapid  rate  of  straining  § 
is  small  compared  with  e  -  e0 ,  although  tx  d§/dt 
is  not.  In  this  case,  (48)  reduces  to  the  dynamic 
load- strain  relation 
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Td  =  f'(e  -  e0)  +  K,  (e  -  e0).  (52) 

The  associated  dynamic  spring  coefficient  is 

("5^)  =  fH(e  "  eo)  +  ^  .  (53) 

oe  d 

Wilson  (1967)  has  suggested  a  simple  power  law 
dependency  of  T  on  (e  -  e0 )  for  dynamic  conditions; 
this  is  of  the  form 

Td  =  »  (e  -  %)n  ,  (54) 

where  Y  is  a  constant.  He  ?inds  that  n  =  3  for 
nylon  rope  and  3/2  for  stranded  steel  cable.  However, 
this  ^elation  appears  to  give  a  reasonable  fit  to  the 
data  only  for  moderate  to  large  values  of  e  -  e0 .  At 
small  (e  -  e0 )  the  data  pertinent  to  the  dynamic  tests 
tend  to  indicate  a  linear  dependency  of  on  (e  -  e0 ) 

since  (dT/de)d  has  a  finite  slope  at  e  =  e0  (see  Fig.  3)- 

An  alternative  is  to  take  f(a)  proportional  to 
an  such  that  reduces  to  a  linear  relation  at  small 

e  -  e0 .  However,  as  will  be  shown,  the  data  of  Paquette 
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and  Henderson  (1965)  indicates  that  an  exponential  form 
for  f(a)  is  more  appropriate,  at  least  for  nylon  rope. 
They  evaluated  the  dynamic  spring  coefficient  at  ten 
mean  tensions  from  the  slope  of  the  stabilized  hysteresis 
loops  (of  small  AT)  for  the  test  illustrated  in  Fig.  4. 

Such  measurements  arc  potentially  more  accurate  than 
inferences  drawn  from  the  large  hysteresis  loops 
analysed  in  Wilson’s  study.  Their  results  are  shown 
as  a  function  of  tension  in  Fig.  9*  These  seem  to  indicate 
a  linear  dependency  of  {dT/de)d  on  T.  A  least  squares 
fit  of  the  form 

(H)d  -  K  +  w. 


yields 


(55) 


K  =  1.4  x  104  lbs 

b  =  26.0  , 

with  a  standard  error  of  estimate  of  (dT/de)d  of 
about  0.7  x  104  lbs  for  the  1/2  inch  diameter  nylon  rope. 


Apparent  spring  coefficient  (10,000  lbs) 
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Fig.  9  Apparent  spring  coefficient  for  one-half 

inch  diameter  nylon  rope  versus  mean  tension, 
based  upon  tests  by  Paquette  and  Henderson 
(1965)  at  +  180  lb  cycling  amplitude;  full 
line  represents  a  least  square  linear 
regression. 
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If  we  accept  the  empirical  fit  (55)  then  from 
(52),  (53),  and  (55)  we  get 

f  "(a)  -  b  f'(a)  =  (K-Kj)  +bK  a,  (56) 

where  a  s  e  -  e0 .  This  leads  to  a  solution  for  f‘(a) 
of  the  form 


f '(a)  =  ce1*-f-S1  a,  (57) 

where  c  is  a  constant  of  integration.  Thus  from  (52) 
we  get 

Td  =  c  eb  <£  -  £o)  .  K/b,  (58) 

however,  Td  must  vanish  when  e  =  e0  (by  definition 
of  e0)  and  hence  we  require  that  c  =  K/b.  Consequently 

Td  =  f  t  e”  (e  ”  -  1  1  •  (59) 

The  limiting  slope  at  zero  e  -  e0  (  zero  Td  )  is 

(dTd/de)0  =  K  .  (60) 
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Moreover  for  the  static  case,  ror  which  T  is  given 
by  (51),  we  get 

T  =  Td  -  K*  (e  -  Co)  .  (61) 

This  has  a  limiting  slope  at  zero  e  -  e0  given  by 

(dT/de)0  =  (K-K2),  (62) 

which  clearly  requires  that  £  If.  . 

Finally  from  (47)  and  (57)  with  c  =  K/b  we  get 
for  the  elastic  energy 

Ee  =  §3  [  e  b<e  "  c«)  .  „  (C  -  £o)  ] 

-  K,  (e  -  c0)§  +  (K./2)  5=  .  (63) 

For  very  small  (e  -  e0 )  this  relation  reduces  to  that 
for  the  simple  Maxwell  system  where  K  is  equivalent  to 
Kq  +  Kj  for  that  system. 

As  an  application  of  (59) >  with  the  coefficients 
given  by  (55),  we  note  that  a  dynamic  tension  of  7000  lbs 
will  give  rise  to  a  quasi-elastic  elongationai  strain 
(e  -  e0)  of  about  ten  per  cent.  However,  under  quasi¬ 
static  conditions  the  elongation  would  be  significantly 


larger,  particularly  if  is  nearly  equal  to  K. 


For  the  general  straining  conditions,  (48)  and  (57) 

yield 

T  =  |  [  e  b<e  -  e°>  -1]  -  iq  §  ,  (64) 

which  is  consistent  with  cSE&/&e  from  (63).  This 
together  with  the  prognostic  relation  for  §,  Eq  (49), 
represents  a  generalized  model  for  the  case  of  constant 

eo . 


4.  The  "Saturation"  Relation 

To  relations  (64)  and  (49)  we  must  add  a  third 
relation  governing  .  Figs .  3  and  4  indicate  that 
for  a  given  e0  there  apparently  exists  a  sort  of 
"saturation"  curve  in  the  T,  e-diagram.  This  corresponds 
to  the  path  which  would  result  for  very  slow  undirectional 
straining  from  the  original  state  to  the  ultimate  point 
(path  0  S  U  of  Fig.  10).  It  is  strictly  an  irreversible 
one-way  path,  1  nr  if  the  load  is  removed  at  some  inter¬ 
mediate  point  S,  even  very  slowly,  the  strain  will 
return  to  some  value  e0  along  a  different  path  SB. 
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Thus  in  principle 

«o  =  *  (Ts)  -  (65) 

the  form  of  the  function  being  readily  deduced  empirically. 
The  data  summarized  by  Wilson  suggest  that  a  power  law 
form 


P  (Ts)  a  Ts“  (66) 

leads  to  a  reasonable  fit,  at  least  for  tensions  less 
than  75  per  cent  of  the  ultimate  value.  For  ropes  m  <  1 
while  for  steel  cable  m  >  1. 

It  is  not  implied  by  the  above  discussion  that  T 

s 

is  an  upper  bound  for  T  for  given  e0 .  The  actual 
T  versus  e  near  the  saturation  level  depends  upon  the 
rate  of  strain.  For  example,  referring  again  to  Fig.  4, 
the  relation  Tg  versus  e  is  shown  by  the  dashed  line 
drawn  through  the  uppei  ends  of  the  final  hysteresis 
loops.  There  is  an  ’’overshoot”  of  this  saturation  curve 
which  occurs  when  the  mean  load  is  suddenly  increased 
after  having  reached  a  stable  hysteresis  loop.  In  order 
to  simulate  this  effect,  the  following  relation  is 
suggested  for  e0  : 
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To  de0/dt 


e0  ^ 


if  e0  <  F(T) 
if  e0  *  F(T) 


(6?) 


where  t0  is  a  second  characteristic  relaxation  time 
for  the  rope  or  cable.  If  T  is  initially  at  the 
saturation  level  for  given  e0  and  if  dT/dt  >  0, 
but  of  very  small  magnitude,  then  e0  will  change 
nearly  in  accord  with  F  (T). 

Suppose,  on  the  other  hand,  that  T  is  initially 
at  some  saturation  level  T  (Fig.  11),  corresponding  to 
e0a,  but  experiences  a  sudden  increase  to  the  value 
after  which  it  remains  constant.  Then  since  e„a  <  F(Tb), 
eQ  must  increase  with  time  approaching  the  final  value 
e0f  =  F(T^)  asymptotically.  Specifically  relation  (67) 
implies  that 


=  Gc 


"  (e< 


eoa) 


,-t/Tc 


(68) 


for  this  case.  The  change  of  e  will  have  essentially 
the  same  behavior.  As  a  test  of  this  relation,  the  data 
of  Fig.  4  was  analysed  for  the  case  where  the  mean  tension 
was  suddenly  increased  from  about  6300  lbs  to  7000  lbs. 

In  this  test  T  was  also  varied  periodically  over  the 
range  7000  +  180  lbs  after  the  increase,  thus  giving 
a  progressive  series  of  loops  which  provide  a, time 
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Pig.  12  Relative  strain  (log  scale)  versus  elapsed 
time,  illustrating  '•/he  creep  phenomenon 
near  the  "saturation"  curve  (based  upon 
data  of  Paquette  and  Henderson,  1965), 
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scale  for  the  creep  phenomenon.  The  period  of  cycling 
for  this  case,  as  deduced  from  the  data  of  Paquette  and 
Henderson,  is  about  ten  seconds.  A  semi-log  plot  of 
e0f  “  eo  (in  relative  units)  versus  time  is  shown  in 
Fig.  12.  This  clearly  confirms  the  exponential  re¬ 
laxation  and  yields  a  value  of  t0  for  the  nylon  line 
of  about  40  seconds. 

5.  Hysteresis  for  Nylon  Rope 

The  measurements  of  Paquette  and  Henderson  (1965) 
for  the  one-half  inch  diameter  nylon  rope  included 
evaluation  of  the  hysteresis  for  five  different  ranges 
of  aT  all  at  a  common  mean  tension  of  2000  lbs  and 
;t  cycling  periods  from  12  to  about  60  seconds.  Their 
results  indicated  that  the  hysteresis  energy  absorption 
per  cycle,  H,  was  nearly  independent  of  the  period  of 
cycling.  The  results  are  summarized  in  Table  1. 

In  terms  of  the  present  load-strain  model,  the 
lack  of  dependency  of  H  on  the  period  of  cycling  should 
imply  that  H  is  near  its  maximum,  so  that  it  is  linearly 
independent  of  id. 
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TABLE  1 

HYSTERESIS  IN  HALF-INCH  NYLON  ROPE 
AT  MEAN  TENSION  OF  2000  lbs 

(After  "aquette  and  Henderson,  1965) 


Tension 
Variation  AT 
(lb) 


Hysteresis 
per  cycle 
(lb) 


Period  of  Cycling 
(sec) 


+  i4o 

0.12 

18,  36 

+  280 

0.52 

12,  16,  40,  37 

+  560 

2.3 

16,  30 

+  1120 

15.2  • 

30,  50 

+  1400 

29.7 

30,  55 

6o 


For  small  amplitude  cycling  about  a  mean  tension 
T0,  the  linear  theory  of  section  2  should  be  applicable 
provided  that  (Kg  +  X*  )  is  taken  as  the  dynamic  spring 
coefficient  for  the  mean  tension  T0 .  From  Fig.  9  with 
T0  =  20001bs  we  get  (Kq  +  Kx  )  =  6.8  x  104  lbs.  From 
the  results  of  section  3  it  was  found  that  Kx  should 
not  exceed  1.4  x  10*  lbs  for  this  same  case.  If  we 
consider  that  H  is  at  its  maximum  given  by  (41)  and 
that  =  1.4  x  10*  lbs  then  we  find  that  at  most 

H  =  6.0  x  10-s  (AT)2  ,  (69) 

for  small  aT,  in  the  case  of  the  one-half  inch  nylon 
line  tested  by  Paquette  and  Henderson.  A  plot  of  the 
data  from  Table  1  is  shown  in  Fig.  13 .  The  full  curve 
is  simply  a  smooth  interpolation.  The  dashed,  curve 
is  a  plot  of  relation  (69).  The  fact  that  it  is 
asymptotic  to  the  full  curve  for  small  AT  seems  to 
be  a  verification  of  the  assumptions  leading  to  relation 
(69) ;  namely  that  H  is  at  its  maximum  for  given  AT 
and  that  is  the  maximum  permissible  value  consistent 

with  the  data  of  Fig.  9*  If  this  is  correct  then 
should  be  equal  to  the  value  for  which  H  has  it  maximum, 
namely  ^/(Kq  +  )  which  has  the  magnitude  of  0.79* 


Fig.  13 


>  500 
AT  (lbs) 


1000 


Hysteresis  (H)  versus  amplitude  of  cyclic 
loading  (AT)  for  one-half  inch  diameter 
nylon  rope;  circled  points  represent 
neasuremerits  by  Paquette  and  Henderson 
(1965);  dashed  line  is  based  upon  the 
model  discussed  in  tkr  text. 
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The  overall  mean  period  of  cycling  for  these  tests  was 
about  30  seconds  and  hence  t1  is  presumably  about  four 
seconds. 

There  are,  however,  two  disturbing  features  of  the 
data  which  are  not  adequately  explained  by  the  present 
model.  First  of  all,  the  appreciable  departure  of  H 
at  large  aT  from  that  given  by  (69)  is  not  adequately 
explained.  Second,  that  data  of  Paquette  and  Henderson 
indicate  that  the  dynamic  spring  coefficient  decreases 
with  increasing  aT  for  fixed  value  of  mean  tension. 

The  non-linear  aspect  of  the  load-strain  relation  does 
not  adequately  explain  these  observations.  A  possible 
explanation  may  lie  in  the  frequency  dependence  of  H 
and  the  apparent  spring  coefficient.  However,  adequate 
tests  of  the  frequency  dependence  are  lacking. 

6.  Summary  and  Possible  Further  Generalization 

At  least  many  of  the  observed  properties  of  ropes 
and  cables  seem  to  be  in  accord  with  the  foregoing 
model  which  is  summarized  here: 

T  =  |  [  e  b<e  “  e°)  *  10  -  K*S  ,  (69) 

Ti  It  =  (€  -  eo  -  §) 


7 


(70) 
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T 

0 


f F(T)  -  e0  ,  if  e0  <  F(T) 
|  0  ,  if  e0  a  F(T) 


(71) 


where  F(T)  =  e0  corresponds  to  the  saturation  relation 
appropriate  to  very  slow  unidirectional  straining. 

It  must  be  emphasized  that  further  critical 
measurements  are  needed  to  check  the  validity  of  the 
model,  particularly  in  respect  to  the  effect  of  -.me 
rate  of  strain.  Measurements  of  hysteresis  over  a 
wide  range  of  cycling  periods  are  required  for  an 
accurate  evaluation  of  the  parameters  t,  and  Kj. 
for  a  particular  rope  or  cable.  Indeed  it  is  quite 
conceivable  that  the  true  behavior  of  ropes  and  cables 
requires  a  more  complex  model,  incorporating  a  discrete 
spectrum  of  relaxation  times  and  associated  K  values. 

In  order  to  clarify  the  above  remark  we  note  that 
relations  (69)  and  (70)  can  be  recast  in  the  following 
integral  form  in  which  §  has  been  eliminated 

00 

T  =  |  (eb<X  -  1)  -  j  a(t  -  \)e“x/Ti  dX  ^  (72) 

o 

where  a  s  e  -  e0 .  The  integral  is  a  particular  con¬ 
volution  of  the  function  a(t).  A  natural  generalization 
of  this  relation  is 
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T  -  |  (eb<1  "  X'  '  [  a(t  -  x)  G(X)  dx  ,  (73) 

in  which  the  generalized  function  G(x)  can  be  repre¬ 
sented  in  the . form 

*. 

n 

0(x)  =  l  (KjA  )  .  (74) 

d=l 

Relation  (72)  is  a  special  case  in  which  n  =  1.  Clearly 
(73)  has  greater  flexibility  in  respect  to  fitting  the 
observed  frequency  dependence  of  hysteresis  in  the  case 
of  cyclic  loading.  For  small  amplitude,  simple  harmonic 
variation  of  T  or  e  at  fixed  e0  and  frequency  uj, 
relation  (73)  can  be  shown  to  yield  the  following 
hysteresis 


H  =  t  h. 


(75) 


in  which  h  (<u)  is  the  transfer  function 


h  (tu)  =  J  G(x)  sin  u»X  dX 


(76) 


and  K*  is  the  dynamic  spring  coefficient.  Using  (74) 


gives 
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n 


h  (uj)  =  £ 


UIT 


K. 


1 


3=1 


[1  + 


(iDTj)8 


(77) 


As  an  example,  suppose  n  =  2  with 

=  Ka  and  t2  =  10  tx. 

In  this  case  the  function  h.^)  possesses  a  very  broad 
range  of  u>  over  which  it  is  nearly  a  constant  (Pig.  14). 
There  is  some  indirect  evidence  that  this  behavior  may 
be  closer  to  reality  than  that  implied  by  the  simpler 
(n  =  1)  model.  However,  lacki^~  really  definitive 
evidence  for  this,  we  will  adopt  the  model  delineated 
by  (69)  to  (70)  in  the  remainder  of  this  study. 
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IV.  ENERGY  CONSIDERATIONS 


1.  Energy  Equation  for  a  Mooring  Line 

The  equation  of  motion  of  a  mooring  line  in  vector 
form  was  given  by  (15)*  If  we  form  the  scalar  produc- 

A 

of  V  with  this  equation  we  obtain 


a 


[-2  (m  +  ma)  |Vja  oa  VT2  ] 


=  V 


dT 

Ss 


A  A  A  A 

w  k-V  +  R*V  , 


(78) 


where  for  simplicity  we  will  regard  m  and  w  as 

A  A 

constants.  The  term  w  k*V  can  be  written  as  d(wz)/dt 
where  z  is  the  elevation  of  a  material  point  of  the 
line  at  s,t  and  wz  is  an  effective  potential  energy 
in  the  water  per  unit  of  s. 

A  ~ 

Now  consider  the  term  V*dT/ds.  This  can  be 
written  as 


(79) 
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Furthermore,  since  V  =  3r/3t  we  can  write 

dV  3 
■3s  =  ‘St 

or 

$  =  4-  (t  |2.) 
os  ot  '  OS7 

using  (22).  Now  dcr/ds  =  1  +  e  where  e  is  the 

elongational  strain  at  s,t.  Hence  we  find 


(32) 
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where  V  is  simply  the  magnitude  of  V.  It  will  be 
recalled  from  the  previous  chapter  that  the  elastic 
energy  Eg  is  a  function  of  at  least  three  different 
strains  (e,  and  e0).  It  follows  therefore  that 

dE  dE  de  dE  d§  dE  de 

e  _ e  e  ,  _ e  _ °  r^\ 

sr  =  srst  +-5§“^t  +  -^7^t"  • 

However,  since  e  and  e0  enter  only  in  the  form 
(e  -  e0 )  then  dEe/de0  =  -  dEg/de .  Also  using  (46) 
we  find 

dE  de  d§  de. 

sr  -  T  si  -  N  -  T3t-  •  t84) 

Consequently  the  energy  equation  finally  takes 
the  form 

It  <ek  +  Ep  +  V  “I?  ti-i)  +  R.y 

-  »  (||)2  - *  •  <85> 


in  which 
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E,  =  i  m  Vs  +  i  in  V  2 
k  2  2  a  n 


(86) 


w  z  , 


where  V  2  =  V®  -  V„ a .  Each  of  these  energy  functions 

n  t 

represents  energy  per  unit  material  coordinate  (s). 

The  kinetic  energy  includes  that  of  the  mooring  line  plus 
that  of  the  effective  added  mass  of  water,  which  is 
related  to  the  normal  component  of  motion  of  the  line. 


2.  Energy  Supply  end  Loss  for  the  Line 

If  ;ve  integrate  (85)  with  respect  to  s  from  the 
anchor  point  to  the  point  of  attachment  on  the  surface 
vessel  we  obtain 


(87) 


L 

(* 

where  |  E  ds  is  the  total  mechanical  energy  of  the  line 
(kinetic,  potential,  and  elastic  energy) .  The  term 

A  A 

(V*T)r  represents  the  rate  of  work  done  on  the  mooring 

A  A 

line  by  the  ship.  The  integral  involving  R*V  is  the 
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net  work  done  by  the  water  on  the  mooring  line  per  unit 
time.  The  remaining  two  terms  represent  rates  of  energy 
loss  (or  more  correctly,  rates  of  energy  conversion 
from  mechanical  to  unavailable  thermal  form).  The’'  are 
related  to  the  non-equilibrium  processes  occurring  within 
the  line  (hysteresis  and  permanent  strain).  The  3atter 
terms  are  always  of  one  sigi:  for  clearly  neither  (da;/dt)2 
nor  de0/dt  can  be  negative.  However,  de0/dt  can  be 
zero  if  the  range  of  T  for  given  e0  ?s  suitably 
limited.  On  the  other  hand  (d§/dt)2  can  vanish  only 
under  static  conditions. 

If  the  water  is  nominally  ac  rest  (U  =  0),  then 

A  A 

the  term  R*V  is  negative  and  represents  a  rate  of 

> 

energy  loss  to  the  water  by  form  drag.  In  this  case 
a  steady  state  of  ^namic  energy  of  the  line  can  be 

A  A 

maintained  only  if  (V*T)L  is  positive.  Clearly  in  the 

A  A  A 

deep  water  where  U  is  negligible  R*V  will  indeed 
tend  to  be  negative  and  hence  lead  to  damping  of  trans¬ 
verse  oscillations  in  that  region  of  the  line.  The 
energy  loss  terms  related  to  (d§/dt)2  and  den/dt  , 
on  the  other  hand,  are  effective  in  damping  longitudinal 
waves  in  the  line. 
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It  can  be  shown  that  for  cyclic  variation  of  T 
the  mean  value  of  N(d§/dt)a  is  simply  H/P,  where 
P  is  the  period  of  cycling  and  H  is  the  hysteresis 
defined  by  (37). 


V.  MOORING  LINE  EQUATIONS  IN  CHARACTERISTIC  FORM 


1.  Geometrical  Considerations 

Let  x,  y,  z  be  a  fixed  Cartesian  coordinate  system 
with  origin  at  the  anchor  point  and  with  z  representing 

A  A  ^ 

elevation.  Let  i,  j,  k  represent  the  associated  unit 
vectors  for  this  fixed  coordinate  system.  The  natural 
coordinates  are  defined  in  terms  of  the  local  orientation 
of  the  mooring  line  at  material  point  s  at  time  t 

A 

(Zajac,  1957).  It  will  be  recalled  that  t  is  a  unit 
vector  tangent  to  the  axis  of  the  mooring  line  at  s,  t 

A 

and  in  the  direction  of  increasing  s.  Now  let  p  be  a 

A  A 

unit  vector  normal  to  k  and  r  with  the  direction  of 

A  A  A  A 

k  x  t  «  Moreover,  let  v-  be  a  unit  vector  normal  to  t 

A  A  A  AAA 

and  \x  with  direction  n  x  t.  Thus  v,  n,  t  form  a 
mutually  orthogonal,  right  handed,  set  of  local  coordinate 
vectors. 

Let  e  denote  the  zenith  angle  and  $  the  azimuth* 

A 

angle  of  the  local  unit  vector  t  relative  to  the  fixed 
Cartesian  coordinates  (see  Fig.  15).  Note  that  these 
angles  are  not  coordinste  angles  of  the  position  vector 

*  More  precisely  is  the  azimuth  of  the  vertical 
plane  containing  vector  t,  as  measured  from  the  x  axis. 
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Fig.  15  Schematic  illustrating  the  local  coo: 

system  for  a  mooring  line  and  its  re: 
to  a  fixed  Cartesian  frame  of  refere] 


of  P  but  rather  of  the  orientation  of  the  tangent  at 
point  P.  From  the  geometry  it  can  be  shown  that  the 
relations  between  the  local  coordinate  vectors  and  the 
fixed  reference  coordinate  vectors  are  given  by 


v  =  i  cos  6  cos  $  +  j  cos  e  sin  $  -  k  sin  e. 


p  =  -i  sin  §  +  j  cos  $  , 


(88) 


t  =  i  sin  6  cos  $  +  j  sin  8  sin  $  +  k  cos  8  . 


The  orientation  of  v,  d,  t  in  general  vary  with  s  along 
the  line  at.  given  t  or  with  t  at  fixed  s.  Using 
relations  (88),  it  is  easily  established  that 


=  p  cos  8  d§  -  r  d0  , 


A  A 

=  -v  cos  0  d$  -  t  sin  8  d$  , 


(89) 


dt  =  v  do  +  4  sin  8  d$ 


2.  Kinematical  Relations  Involving  the  Dependent  Variables 

A 

The  vector  V,  which  represents  the  velocity  of  a 
material  point  of  the  mooring  line,  can  be  represented 
in  the  form 


=  Vv  v  +  p  +  Vt  t  , 


V 


(90) 


in  which  V  and  V  are  normal  components  and  V 
v  \i  t 

is  the  tangential  component  of  the  velocity.  The 
acceleration  of  a  material  point  is  given  by 

dV  dVv  .  dV  .  dV  . 

•St  =  •ST'v+Wi^  +  -StLT 

+  vv  +  vp  +  vT  .  (91) 

A 

However,  the  turning  vectors  like  d^/dt  can  be  evaluated 
in  terms  of  d0/dt  and  d#/dt  using  (89).  The  resulting 
relations  for  the  acceleration  is 

dV  dV  *  dV  *  dV  . 

+  Vv  (cos  e  ||  *  -  ||  t) 

+  \  (-  cos  0  ||  V  -  Sin-  0'||  T  ) 

+  Vt  (H  v  +  sin  0  It  »)  *  (92) 


This  relation  will  be  useful  for  recasting  the  vector 
equation  of  motion  into  the  natural  component  form. 
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A  relation  entirely  similar  to  (92)  can  be  obtained 

A 

for  dV/ds,  in  which  the  time  derivative  terms  are  re¬ 
placed  by  derivatives  with  respect  to  s.  However,  we 

A 

have  already  found  that  bV/bs  is  related  to  be/bt  and 

A  A 

br/bt  by  relation  (80).  If  we  now  write  dr/dt  using 
the  last  of  relations  (89)  and  equate  coefficients  of 
corresponding  unit  vectors  in  the  two  different  expressions 

A 

for  bV/bs ,  we  find  the  following  (prognostic)  relations 
for  e,  $,  and  e: 


bVs 


(x  +  e)  3t  =  ”55”  +  V* 


t  ds 


-  \  cos  e  &  ,  (93) 


(1  +  e)  sin  9  ||  =  -g—  +  (Vv  cos  e 


dV 


+  Vt  sin  6)  ||  , 


(94) 


5e 


bV_ 


”5t  =  ”55”  ’  sin  e  -  Vv  ds 


&  • 


(95) 


These  are  purely  kinematical  compatability  relations 
among  the  dependent  variables  9,  $,  e,  V,V,  and  V^. 
Moreover  they  are  relevant  only  for  the  transient  state; 
in  the  case  of  equilibrium  conditions  they  are  satisfied 
trivially. 
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Relation  (95)  can  be  interpreted  in  the  following 
we.y.  Suppose  we  solve  for  dVT/ds  and  integrate  over 
the  full  range  of  s.  This  yields  (since  Vt  (o,t)  =  0): 

L  L. 

Vt(l,  t)  =  |  |f  ds  +  J  (V^  sin  e  |f 

+  Vv  |f  )  ds*  (96) 

As  a  special  case*  suppose  dg/ds  and  dj/ds  vanish. 

This  implies  that  the  line  is  straight  and  hence  (96) 
implies  that  a  positive  tangential  velocity  can  exist 
at  the  end  of  the  line  only  if  the  line  stretches.  On  the 

other  hand  if  the  line  is  curved,  then  part  of  Vt 

at  the  end  of  the  line  can  be  due  to  the  effect  of  "talcing 

up  slack"  in  the  line$  this  is  represented  by  the  second 
integral  in  (96). 

Relations  (93)  and  (94)  can  be  interpreted  in  a 
somewhat  analogous  manner.  Specifically  the  line  can 
have  a  normal  component  of  motion  at  its  upper  extremity 
if  the  line  has  a  progressive  tilt  with  time  or  if  it 
experiences  a  flexing. 
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The  remaining  kinematical  relations  of  importance 
are  those  relating  the  Cartesian  coordinates  of  a 
material -point  of  the  line  to  the  natural  coordinates. 
First  of  all  we  note  that 


dr 

*5t 

A 

where  V 


=  V  ,  (97) 

is  expressed  in  natural  coordinates  by  (90)  and 


r  =  xi  +  y]  +  zk.  (98) 

Using  (88),  (90),  (97),  and  (98)  we  find 


dz 

Zt 


=  cos  9  cos  $  -  sin  i  +  sin  0  cos  $,  (99a) 
=  Vv  cos  0  sin  i  +  cos  $  +  Vt  sin  e  sin  $,  (99b) 
=  -Vv  sin  0  +  VT  cos  e  .  (99c) 


We  can  also  employ  (22)  to  obtain  relations  for  the 
derivatives  of  x,  y,  and  z  with  respect  to  s.  Using 
(22)  and  (88)  yields 
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dx 


(1  +  e)  sin  0  cos  § 


(100a) 


(1  +  e)  sin  0  sin  $ 


(100b) 


=  (1  +  e)  cos  e 


(100c) 


For  the  case  of  steady  equilibrium  conditions,  the  latter 
set  of  relations  are  useful.  Either  set  of  relations  can 
be  used  in  the  transient  problem. 


3.  Equations  of  Motion  in  Natural  Coordinates 

The  vector  equation  of  motion  (15)  contains  the 


dT/ds 

which 

can  be  expressed 

A 

dT 

3s  = 

dT  : 
3Ft 

+  T&. 

(101) 


where  br/bs  is  the  local  vector  curvature  of  the  line. 
Using  the  last  of  relations  (89)  we  find 


H  - 


Also  from  relations  (88)  we  get 

A  A  A 

k  =  -  sin  0  v  °+  cos  0  t  . 


(102) 


(103) 
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The  residual  force  R,  given  by  (18),  contains  the 
terms  dUT/dt  which  can  be  rewritten  in  the  form 


auT 

ST 


au 

St 


T  +  U 


A 

8t 

•St  * 


(104) 


It  will  be  convenient  to  denote  by  G  that  part  of  R 

a 

which  does  not  depend  upon  dT/dt.  Accordingly 


R 


G 


(105) 


and 

G  =  x  (Un  -  Vn)  +  r  ;  X  (Un  -  Vn) 

r  \  au  m  A  ,au 

+  (md  +  ®a)  "  ma  T  (“St  *  )  ' 

where 

x  *  (p/a)  Cd  d  (l  +  c )  |un  -  Vn|  , 

r  =  X  T)  cos  (u)s  t  +  6)  . 

Finally  we  note  in  analogy  to  (102)  that 

A 

aT  8a  .  a  $  a 

-5t  =  #  v  +  sin0^t^- 


(106) 


(107) 


(100) 
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Consequently  using  (15),  (92)  and  the  above  relations 
of  this  section  we  obtain  the  following  component  equations 
of  motion 


dV, 


<ra  +  raa)  |  ST  +  Vt  It  "  Vn  cos  0  ll  \ 

-  T  ||  =  w  sin  0  +  Gv  , 

fdV  v  ' 

(m  +  raa)  +  (Vv  cos  8  +  VT  sin  e)  -^|- 


-  T  sin  9 


G  , 


(109) 


(110) 


and 


m 


dV 

_ T 

"ST 


-  [■  Vv  +  ma  <Vv  -  V  1  ll 

-  t“  \  +  ma  -  V  1  s“  9  ll 


dT 


-  *3s  =  -  w  cos  e  +  Gt  . 


(Ill) 


The  components  of  G  in  the  natural  coordinate  system 
are  as  follows: 
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Gv  -  X  (Uv  -  Vv)  -  r  %  -  V,x)  +  (md  +  » a)^.f 
GH  =  *  (Up  -  V  +  r  (Uv  -  Vv)  +  (md  +  na)  ,  (112) 


where  Uv>  U^,  UT  are  the  natural  coordinate  components 
of  the  vector  dU/dt. 

a 

It  is  supposed  that  the  U  is  prescribed  in  terms 
of  its  Cartesian  components,  Ux,  U^,  Uz.  Using  (88) 
we  find  that 


U 


U 


UT 


U  cos  e  cos  i  +  U  cos  e  sin  $  -  U  sin  e 
^  y  z 

-  Ux  sin  i  +  Uy  cos  i  ,  (113) 

U  sin  0  cos  $  +  U  sin  e  sin  $  +  U  cos  0. 
a  y  z 


Moreover  it  is  supposed  that  U  ,  U  ,  U  are  prescribed 

X  y  z 

functions  of  x.  v.  z.  ^nd  t.  Hence  they  are  implicitly 
related  to  s,  t  through  the  dependent  position  coordinates 
of  the  mooring  line.  The  functions  U  ,  U  ,  and  XL 

V  \i  T 

are  related  to  their  Cartesian  counterparts  U  ,  U  and 

x  y 

U  through  equations  similar  to  (113). 
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Finally  we  note  in  passing  that  the  evaluation  of  the 
coefficient  \  requires  the  relation 

l«„  -  v„|  =  [  (  Uv  -  Vv)2  +  (0„  -  v/  ]1/2  .  (ll4> 


4.  Characteristic  Form  of  the  Mooring  Line  Equations 

in  Natural  Coordinates 

In  the  natural  coordinate  representation  we  have 
the  following  dependent  variables:  Vv,  V  ,  VT,  e,  $,  T, 
and  e.  Each  "  these  is  a  function  of  s,  t.  They 
are  governed  by  the  kinematic al  relations  (S3)*  (94),  (95) 
and  by  the  dynamical  relations  (109)*  (HO)*  (111).  To 
these  relations,  we  add  the  diagnostic  equation  of  state 
(69)  relating  T  to  «;  plus  the  two  additional  variables 
§  and  e0,  the  latter  two  being  governed  by  the  prog¬ 
nostic  relations  (70)  and  (71).  This  gives  a  total  of 
nine  dependent  variables  and  nine  equations.  We  can 
also  add  to  ...ese  the  three  relations  (99a,  b,  c)  or 
(100a,  b,  c)  governing  the  three  position  coordinates 
x,  y,  z. 

Thus,  in  summary,  the  genera.l  three-dimensional 
transient  mooring  line  problem  as  formulated  here  (aside 
from  upper  end  conditions )  involves  eleven  prognostic 
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relations  involving  the  twelve  variables 

9 ■*  §»  ®o-»  y>  ^ 

plus  one  diagnostic  relation  involving  T,  e,  §,  and  e0  . 
Each  of  the  prognostic  equations  is  of  quasi-linear  form. 

For  example,  referring  to  (109)  to  (111)  we  note  that  all 
derivatives  of  the  dependent  variables  enter  linearly, 
their  coefficients  being  functions  of'  the  dependent 
variables  but  not  of  their  derivatives.  Likewise  the 
terms  Gv,  G^,  and  CGt  are  functions  of  the  dependent 
variables  but  not  of  their  derivatives.  Any  quasi-linear, 
hyperbolic  system  of  n  equations  in  n  dependent 
variables  but  involving  only  two  independent  variables 
(s,  t  in  this  case)  can  be  recast  in  a  characteristic 
form,  which  facilitates  integration  of  the  system  (see 
for  example.  Freeman,  1951). 

The  transformation  leading  to  the  characteristic 

© 

form  of  the  equations  is  outlined  in  Appendix  B.  The 
resulting  equations  are  as  follows : 

dV  ± 

dt^  +Bl  I  =  G-/(m  +  ma)  (H5a) 

along 


1 


9 
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where 


along 


where 


H  =  ±  [T/(m  +  ma)  (1  +  «)  31/a  ,  (115b) 


B,  =  (V  cos  0  +  V  sin  e) 
1  '  v  t  ' 


+  [  (1  + 

e)  T/(m  ma)  ]1^2  sin  8  5 

(115c) 

dVv 

.  B±  dQ- 

(Gv  +  w  sin  e) 

dt 

+  *2  dt 

(m  +  maJ 

+  V  cos 
U 

9  dt  ' 

( ll6a) 

ff  =  +  t  T/(m  +  ma)  (1  +  e)]1/2  ,  (H6b) 


B2  =  Vt  +  [  (1  +  e)  T/(m  +  ma)]1'/2  j'  (116c) 


x  oa/a  dT  -  w  C0S  8)  ± 

rH-.  H+ .  “  m  +  r 


+  dV  +  dV 


+  B3  StT  + 


+  CIV  +  +  ia- 

b4  St*  +  B5  df  +  B6  dt  ’ 


(117a) 


along 


ds  /v  /  \ 1/2 

«  ■  ±  <Y/m) 


(117b) 


where 


Y  =  K  e 


b  (e  -  e0) 


(117c) 


~t~  ~t~  +  + 

and  B^,  Bg,  P  are  functions to  be  defined 

later;  finally 


(e  -  e0  -  §) 


(ll8a) 


along 


at  ■ 


(li8b) 


-  -  .1  -- 


-<  1  -I,.  «(8»e(i6ftOT'yi2a 
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and 


dep 

dF 


S(T,  c0)  , 


(H9a) 


along 


ds 

dt 


0  , 


(119b) 


where 


S(T,  e0)  = 


(F  (T)  -  c0)/ 't0,  if  e0  <  F  (T) 

0  ,  if  ec  *  F  (T).  (119c) 


Moreover,  from  (69) 


e 


K  +  bT  +  biq  § 

- K -  ), 


(120) 


which  implies  that  Y  has  the  alternative  form 


Y 


K  +  bT  +  biq  l  . 


(121) 


•f 

Moreover,  P  is  given  by 

I  y 

P-  =  +  (y/m)1/2  j  ^  (e  -  €0  -?)  +  S(T,  e„) 

-  s!  [  (\  -  V  («  6  +  eg  +  (vu  -  ug  G|i] .. 

* 

(127) 

To  these  relations  we  add  equations  (99a,  b,  c)  for 
prediction  of  the  Cartesian  position  coordinates  of  the 
line  at  s,t.  An  alternative  is  to  employ  (100a,  b,c) 
to  evaluate  these  quantities  at  a  given  instant  from  q 
and  i . 

Relations  (115)  and  (116)  govern  the  transverse 
wave  modes  in  the  line,  while  relations  (117)  govern  the 
longitudinal  wave  modes,  the  damping  of  which  are 
governed  by  (118)  and  (119).  The  nominal  signal  speeds 
for  the  two  modes  are: 
for  transverse  modes, 

C  *  =  [T/(m  +  ma)  (1  +  e)  ]1/2  ,  (128) 

for  longitudinal  modes, 

C”  =  [Y/m]1''2  . 


(-29) 
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Note,  however,  that  these  refer  to  speeds  in  the  sense  of 
ds/dt.  The  true  speed  in  the  vertical  direction  can  be 
evaluated  from  the  relation 

§  -  It  +  §c>  (^o) 

where  C  is  the  value  given  by  (128)  or  (129)  depending 
on  the  wave  mode.  Using  relation  (99c),  (100c)  and  (130) 
gives 


=  (1  +e  )  C  cos  6  -  Vv  sin  6  +  VT  cos  e  ,  (131) 

along  ds/dt  =  C . 

The  transverse  and  longitudinal  wave  modes  are 
coupled  through  certain  non-linear  effects.  In  particular 
the  longitudinal  waves  are  coupled  to  the  transverse  waves 
through  the  terms  involving  the  time  rates  of  change  of 
V^,  V  ,  6  and  $  which  appear  on  the  right  hand  side 
of  (117a).  These  coupling  effects  are  significant  only 
if  d@/dt  and  d$/dt  are  large.  Now  since 

S§  -  H  +  CH- 


9 


(132) 
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where  C  =  +  (Y/m)'1'^  for  this  case,  it  follows  that 

d8/dt  can  be  large  if  the  curvature,  d0/ds,  is  large 
even  though  the  local  rate  of  change  of  0  may  be  small. 
In  a  similar  manner,  df/dt  can  be  large  if  the  pitch, 
dj/ds,  is  large.  It  is  also  possible  for  the  angular 
velocity  d#/dt  to  be  large  for  rotary  motion  of  the 
line.  We  note  as  well  that  the  v-component,  transverse 
wave  modes  can  be  coupled  to  the  p-componer.t,  transverse 
modes  if  d$/dt  is  large. 

5.  Magnitude  of  the  Signal  Speeds 

The  upper  limits  on  C’  and  C"  are  determined 
respectively  by  the  ultimate  values  of  T  and  Y.  Wilson 
(1967)  gives  the  following  values  for  6  x  37  galvanized 
steel  cable  of  nominal  diameter  DQ  (inches): 

m  =  4.9  x  10"2  Dq2 

Tu  =  7.3  x  104  D02 

Yu  =  6  X  106  D02 

e  =  0.03  • 

u 


(slugs/ft) 

(lbs) 

(lbs) 
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For  three  strand,  regular  lay  nylon  rope: 

m  =  0.9  x  10-2  DQ2  (slugs/ft) 

Tu  =  3.1  X  104  Do2  (lbs) 

Yu  =  8  x  105  Dq2  (lbs) 

su  =  0.52  . 

The  latter  three  values  are  taken  from  Paquette  and 
Heiderson's  data. 

The  value  of  m  is  given  by 

cl 

ma  “  ca  P  5  C1  +  £)  * 

2  2 

If  we  assume  that  D  (1  +  e)  =  D  and  that  C  =  1, 

'  •  o  a 

v;e  get 


ma  =  1.1  X  10-2  D02  (slugs/ft), 

where  D  is  in  inches, 
o 

Using  (128)  and  (129)  with  the  above  data  yields 
the  maximum  values  of  the  signal  speeds  as  presented 


in  Table  2. 
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TABLE  2 

MAXIMUM  SIGNAL  SPEEDS 


Steel  Cable 

Transverse  (C^)  1010  ft/sec 

Longitudinal  (C^)  11,000  ft/sec 


Nylon  Rope 
1000  ft /sec 
9400  ft/sec 


These  speeds  are  independent  of  the  diameter  of  the  rope 
or  cable.  Also  the  maximum  speeds  for  nylon  rope  are 
seen  to  be  very  close  to  the  corresponding  speed  for 
steel  cable.  Finally  we  note  that  the  longitudinal 
speed  is  about  a  factor  of  ten  larger  than  the  speed 
for  transverse  waves. 


6.  Characteristic  Form  of  the  Mooring  Line  Equations 
in  Cartesian  Coordinates 

The  characteristic  relations  in  natural  coordinates 
as  given  in  (115),  (116),  and  (117)  have  a  serious  potential 
pitfall  which  can  lead  to  difficulties  in  application. 

For  those  conditions  where  0=0,  the  solution  for  $ 
becomes  indeterminate.  Note  that  B^  =  Bf  when  9=0, 
so  that  only  one  relation  is  obtained  from  (115).  Even  if 
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we  were  to  apply  these  relations  only  to  those  cases 
where  the  mean  value  of  0  differs  from  zero,  it  is  still 
possible  for  the  instantaneous  value  of  0  at  a  particular 
point  s  to  be  zero  due  to  the  presence  of  transverse 
waves. 

This  difficulty  can  be  circumvented  by  rewriting 
relations  (115)  to  (H7)  in  Cartesian  form.  The  details 
of  ohe  transformation  are  lengthy,  but  straightforward. 

This  is  outlined  in  Appendix  C.  In  principle  we  use 
relations  (88)  plus  the  definitions 


a  =  cos  5  sin  0  , 

3  =  sin  §  sin  9  ,  (132) 

Y  =  cos  0  . 


Note  that  a,  3,  y  are  simply  the  Cartesian  components 

A 

of  the  vector  t  (see  (88)  ).  Moreover,  they  must 
satisfy  the  relation 


a2  +  32  +  v2  =  1 


(133) 


and  hence  only  two  of  these  components  are  independent. 

Let  Vv,  Vv,  V7  represent  the  Cartesian  components 
~  •» 

of  V.  Clearly  these  are  related  to  x,  y,  z  by 


(134) 


dx/at  =  Vx, 

ay/at  =  vy, 

az/at  =  vz  . 

Also,  a,  p,  Y  are  related  to  x,  y,  z  and  e  by  the 
relations 

ax/as  =  (1  +  e)  a  , 

ay/as  =  (1  +  e)  0  ,  (135) 

az/as  =  (i  +  s)  Y  - 

Note  that  these  satisfy  (133)  automatically  since 


(i  +  e)3  =  (dx/6s)a  +  (dy/as)2  +  (az/as)2  .  (136) 


The  transformed  versions  of  (115)  to  (117)  in  terms 
of  Vx,  Vx,  Vz,  a,  p,  Yj  e,  and  T  can  be  shown  to  have 
the  following  form: 


(1  -  of ) 


dV 

_ _ > 

dt 


-  ap 


dV 

__Z  - 

dt 


aY 


dV 

_ 2 

dt 


+  pL+_dT]V2  da 

l,  ®  j  dt 


(°nx  +  w  a  y)/“' 


(137a) 
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along 


ds 

dt 


=  +  [  T/m  ’  fl  -f  e)  ]1//2  ; 


(137b) 


JV„  dV  dV  _ 

*  ?  +  (1  "  p2)  dt  “  Py  dT-  + 


r.a.+..?.)T 

I  m» 


1  1/2 


la 

dt 


=  (Gny  +  w  3  y)/m’  ' 


(138a) 


along 


ft  -  +  [?/*'  (1  +  E)]1/2  ; 


(138b) 


and 


a 


dV 
_ _ x 

dt 


+  0 


dV 

__Z 

dt 


dV 


ZT 


V  \  -L/c:  UJ. 

_+  (rnY)  hT 


-1/2 


aT 

dt 


=  (Gt  -  w  Y)/m  +  (Y/m)1/2  [Kx  (e  -  €o  -  ?)Atx 
+  S  (T,  e0 )  ]  +  6—  ,  (139a) 


along 


ds 

dt 


±  (Y/m) 


1/2 


(139b) 


Here  we  have  let 


aft 


m*  =  m  +  mQ  ,  (l4o) 

cL 

for  convenience.  The  terms  G  „  and  G_„  denote  the  x 

nx  ny 

A  A 

and  y  components  of  the  vector  G  -  T  GT  •  The  term  5— 

denotes  a  residual  of  secondary  importance  which  will 

hereafter  be  ignored;  it  is  given  in  Appendix  C  and  the 

justification  of  neglecting  it  is  discussed. 

The  quantities  ,  G  „  and  G  in  terms  of 

nx  ny  t 

Cartesian  variables  are  as  follows 

Gnx  -  X  [  (1  -  a‘)(Ux  -  VK)  -  aS(Uy  -  Vy)-ay(Uz  -  V,)] 

+  r  te(trz  -  Vz)  -  Y(uy  -  Vy)] 

+  (md  +  ma)  [  (1  -  aa )  Ux  -  afi  Uy  -  ay  «Z1  ,  (Mia) 
Gny  =■  \  [-  a9(Ux  -  Vx)+  (1  -  f|=)(Uy  -Vy)  -  SY(UZ  -  Vz)] 

+  r  [Y(UX  -  Vx)  -  a(G2  -  VZ)1 

+  (»>d  +  «a)  [-  o3  Ux  +  (1  -  0=)  6y  -  »Y  0.)  ,  (141b) 


(141c) 
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In  the  evaluation  of  X  the  term  |U  -  V  !  now  must  be 

1  n  n' 

determined  from  the  relation 

I'n-’nl  ■  {[(l-a=)(Ux-Vx)-ae(Uy-Vy) 

-  av(Uz  -  Vz)]2+  [-ap(Ux  -  Vx) 

+  (1  -  e2)(Uy  -  Vy)  -  By(Uz  -  Vz)]2|1/2.  (142) 

Equations  (137)  to  (139)  plus  (118)  and  (119)  involve 

the  prognostic  variables  V  ,  V ,  V  a,  p,  T,  §,  and  e  . 

a  y  z  u 

All  other  dependent  variables  like  e,  y,  Y  are  related  to 
these  through  the  appropriate  diagnostic  relations  already 
stipulated.  The  Cartesian  form  of  the  relations  have  the 
advantage  that  the  relations  remain  determinate-  when  a 
or  p  vanish  (i.e.,  0=0).  Moreover,  they  do  not  require 

the  evaluation  of trigonometric  functions,  which  can  be 
of  economic  advantage  in  sequential  numerical  integration 
operations.  However,  (137)  to  (139)  do  have  the  disadvantage 
that  the  time  derivatives  of  Vx,  V  ,  and  Vz  enter  in 
all  relations.  This  implies  that  the  numerical  integration 
procedure  based  on  (137)  to  (139)  is  more  complex  than  that 
based  on  (115)  to  (117).  Clearly  this  is  the  price  one 
must  pay  to  be  assured  of  a  determinate  system  under  all 
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conditions.  The  numerical  method  of  solution  based  on  the 
method  characteristics  for  the  Cartesian  variables  is 
outlined  in  the  subsequent  section. 

It  vdJl  be  remarked  in  passing  that  relations  (137) 
to  (139)  can  be  derived  from  the  original  equations  of 
motion  in  Cartesian  form  by  a  procedure  similar  to  that 
outlined  in  Appendix  B.  This  procedure  in  fact  was  em¬ 
ployed  by  the  writer  as  a  check  on  the  transformation 
discussed  above,  but  for  a  less  general  case. 

7.  Numerical  Procedure  Using  the  Method  of  Characteristics 

We  consider  in  this  section  a  numerical  method  of 
evaluation  of  the  prognostic  variables  which  is  adaptable 
to  high  speed  digital  computers. 

Let  Q(<j,k)  denote  any  of  the  dependent  variables 
at  discrete  values  of  s  and  t  defined  by 

s j  —  J  j  ~  0,  1,  2.  .  .  .N 

(143) 

ti^  =  k  ^t,  k  =  0,  1,  2.  ... 


where  AS 


L/N  . 
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The  set  of  prognostic  variables  includes  Vx,  V  , 

V  ,  a,  T,  e0  •  It  will  tie  supposed  that  each 

of  these  is  stipulated  for  each  j  at  k  =  0  as  initial 
conditions.  The  strain  e  can  of  course  be  obtained  from 
(120)  and  the  initial  position  coordinates  x,  y,  z  for 
each  j  can  be  obtained  by  simple  numerical  quadrature 
using  (135)  with  the  condition  x  =  y  =  z  =  0  at  j  =  0. 
Alternatively  we  could  supply  x,  y,  z  for  each  j  at 

k  =  0  and  evaluate  a,  3,  e  using  (135)  then  get  T 
from  (69).  The  main  point  is  that  only  eight  of  the  above 
variables  can  be  specified  arbitrarily  at  time  zero. 

As  end  conditions  at  the  sea  bed  we  consider  that 
x(0,k),  y(0,k)  z(0jk)  are  zero  for  all  k.  For  the  upper 
end  conditions  we  will,  for  the  time  being,  suppose  that 

V  (N,k),  V (N,k),  and  V  (N,k)  are  specified  for  all  k 

a  y  Z# 

as  fluctuating  values  with  zero  mean  value  (no  net  motion 
over  an  extended  period  of  time). 

The  discussion  of  the  computational  procedure  at  a 
general  time  step  is  facilitated  by  reference  so  Fig.  l6. 
We  assume  that  V  ,  V  ,  V  ,  a,  3,  T,  §,  and  e0  have 

A.  J  Z, 

already  been  evaluated  for  each  j  at  time  level  k. 

We  wish  to  evaluate  these  same  variables  at  the  next 

time  level  k  +  1  for  a  common  value  of  j  (point  B) . 

The  lines  A  B,  n  =  1,  2.  .  .  .8,  denote  the  character- 

n 

istic  paths  with  slope  ds/dt  =  (-1)  +  C 


associated 
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with  the  eight  different  prognostic  relations  obtained 


from  (137a), 

(138a), 

(139a),  (ll8a)  and  (119a). 

Specifically 

Cn  =  < 

~°n 

°n 

0 

for  n  =  1,  2,  3,  4 

for  n  -  5,  6 

for  n  =  7,  8 

(144) 

where 

C '  = 

(T/m*  (1  +  e))1/2 

(145a) 

c"  = 

(Y/m)1/2  , 

(145b) 

Y  = 

K  +  bT  +  bKx  §  . 

(145c) 

The  subscripts  on  C'  and  C"  imply  appropriate  average 
values  of  these  variables  for  the  path  in  question,  the 
first  estimates  being  interpolated  values  at  the  points 

An* 

The  numerical  computation  scheme  described  below 
involves  a  successive  approximation  procedure  similar 
to  that  employed  by  O'Brien  and  Reid  (1967).  This 
scheme  requires  appropriate  interpolational  relations  for 
the  variables  at  points  An.  Let  an  represent  the  values 
of  s/as  corresponding  to  points  An.  If  the  are 

appropriate  average  values  then 


an  =  i  +  (-l)n  Cn  At/As. 


(146) 


For  values  of  J  near  the  end  points  it  is  possible  for 

an  to  fall  outside  the  admissible  range  of  j.  In  this 

case  special  procedures  are  required  and  are  allowed  for 

in  the  following  discussion. 

Let  J  be  the  integer  truncation  of  a^  such  that 
n  n 

if  an  >  0, 


n 


n 


jn  +  i. 


(147) 


Moreover  if 


Os  an  <  N  , 


(148) 


then  a  given  variable  Q  can  be  estimated  for  points  An 
by  linear  interpolation  as  follows: 


«  ( V  k)  A  (an  -  Jj  9  (Jn  +  1,  k) 


!•  (Jn  +  1  -  »n)  »  (Jn.  k). 


(149) 
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i 

C  o 

i 

i  i 

k  k~l 


<  0  or  a^<  0 


a2  >  N  or  a..  >  K 


Fig.  17  Schematic  illustrating  special  conditions 
for  the  characteristic  paths  near  the  end 
points  of  the  sysxem. 


If  (148)  is  not  satisfied,  this  implies  that  the 
characteristic  line  in  question  intersects  the  boundary 
j  =  0  or  j  =  N  at  k  +  bn  where  0  <  b^  <  1  (see 
Fig.  17).  The  value  of  bn  can  be  evaluated  from 

bn  “  <a„  -  Je>  <an  *  ■»  •  U50) 

where  J  =  0  if  &  <  0  or  J  =  N  if  a  s  N  .  The 
e  n  e  n 

value  of  variable  Q  at  point  An  (Fig.  17)  can  be 
estimated  by  the  interpolational  relation 

Q  (Je,  k  +  bn)  &  bn  Q  (Je,  k  +  1) 

+  (1  -  V  Q  (Je,  k).  (151) 

Clearly  this  requires  that  Q  be  known  at  the  end  points 
for  the  time  level  k  +  1.  Accordingly  the  calculations 
of  all  eight  prognostic  variables  should  be  carried  out 
for  the  end  points  first,  making  use  of  the  end  conditions 
to  replace  those  relations  corresponding  to  characteristic 
paths  for  which  aR  falls  outside  the  allowable  range. 

The  discussion  of  the  numerical  versions  of  the 
prognostic  relations  is  greatly  facilitated  by  employing 
the  following  index  notation  to  represent  the  prognostic 


variables  : 
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Sn  =  -V  V  V  a'  T'  5'  e°}  '  (152) 

where  m  =  1,  2,  ...  .8.  Thus  Qn  =  V  ,  Q,  =  V 

x  -  y 

and  so  on.  The  prognostic  equations  all  contain  terms 
of  the  type  B  d  Q^/dt  which  can  be  approximated  by 


2  [B  (<3>  k  +  1)  +  B  (An)]  [Q^  (j,  k  +  1)  -  (An)]/6tn  , 


where  B(An)  or  Qm(An)  denotes  the  value  of  the 
variable  concerned  at  point  '  An  arid  §t  denotes  the 
time  step  appropriate  to  path  n.  Specifically  for  n  odd 


(An)  - 


(an<  k) 


if  an  s  0 


S,(0.  6n)  if  an  <  0  ’ 


(153a) 


while  for  n  even 


«m  (An>  * 


(a  ,  k} 


if  an  <  N 


^  (N,  k  +  bn)  if  an  2  N  , 


(153b) 
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while 


6tn  =  ^ 


At 


(1  -  bn)At 


if  0  s  a  <  N 
n 


if  a  <  0  , 

n  ’ 


(153c) 


or  an  ;>  N  . 


By  employing  the  above  approximations  in  each  of  the 
prognostic  equations  (137a),  (138a),  (139a),  (118a)  and 
(119a)  we  obtain  eight  algebraic  relations  for  the 
eight  unknown  at  J,  k  +  1.  These  can  be  ..ritten 

in  the  systematic  form 

8 

k)  (j,  k  +  1)  *  Pn(3,k)  ,  (154) 

m=l 

for  n  »  1,  2,  .  .  .8.  The  matrix  of  coefficients  M 

n,m 

and  the  terms  Pn  depend  upon  the  variables  at 

level  k  and  level  k  +  1.  These  functions  are  given 

in  Apnendix  D.  We  note  that  the  matrix  M 

n,m 

following  structure: 


has  the 
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3 

II 

H 

2 

3 

4 

5 

6 

7 

8 

n  =  1 

M11 

M12 

M13 

-1 

0 

0 

0 

0 

2 

M21 

M22 

m23 

1 

0 

0 

0 

0 

3 

M31 

«32 

M33 

0 

-1 

0 

0 

0 

4 

m41 

m42 

M43 

0 

1 

0 

0 

0 

5 

M51 

M52 

M53 

0 

0 

-1 

0 

0 

6 

m6i 

cvj 

s° 

M63 

0 

0 

1 

0 

0 

7 

0 

0 

0 

0 

0 

0 

1 

0 

8 

0 

0 

0 

0 

0 

0 

0 

1 

Thus  (15*0  can  he  reduced  to  the  simpler  system: 

3 

£  (M2n,m  +  M2n-l,m)  Sa  =  (P2n  +  P2n-1^  * 
m=l 

3 

Si+3  =  ^P2n  P2n-l)  ~J^M2n,m  "  M2n-l,m^  %i’ 
for  n  =  1,  2,  3j  and 

Q7  -  P7  ,  *8  -  p8  *  (155°) 


(155a) 

(155b) 
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delation  (155a)  represents  three  equations  in  the  three 
velocity  components  at  j,  k  +  lj  this  set  is  easily 
solved  by  elimination  or  by  Cramer’s  rule.  Once  Qi , 

Qa  are  obtained,  then  a,  p,  T  are  readily  evaluated 
from  (155b)  and  §,  e0  are  given  by  (155c). 

In  carrying  out  the  calculations  it  is  understood  that 

first  estimates  of  the  C  ,  P  ,  and  J4  _  are  obtained 

n  n  n,m 

by  taking  k  +1)  on  which  they  depend  (Appendix  D) 

as  the  values  0^(0, k).  Relations  (155a,  b,  c)  then  give 
first  estimates  of  (^(j,  k  +  1);  the  latter  are  in  turn 
used  to  re-evaluate  the  terms,  C  ,  P  .  and  M  .  Re- 
peated  use  of  (155a,  b,  c, )  with  the  new  coefficients, 
then  yields  better  estimates  of  the  0^(0,  k  +  1).  This 
process  is  repeated  until  the  difference  of  successive 
estimates  is  less  than  some  acceptable  fractional  error. 
Three  or  four  successive  approximations  usually  suffice 
for  sufficiently  small  time  step  At  (O'Brien  and  Reid, 
1967). 

At  the  end  points  j  =  0,  N  we  have  assumed  that  the 
velocity  components  are  stipulated.  At  j  =  0  we  employ 
only  those  prognostic  relations  for  which  an  2  0.  These 
correspond  to  the  partial  set  n  =  2,  4,  6,  7,  8  of 
relation  (154).  Since  the  velocities  Qx,  Qg ,  Qg  are 
given,  one  can  easily  solve  for  the  remaining  unknowns. 


•iwfsivjfc'WtP'Hjtii  mwwwmu 
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At  o’  =  N  we  must  employ  the  prognostic  r  lations  for 
which  an  <  N;  these  correspond  to  the  pa:  .al  set 
n  =  1,  3,  5,  7,  8.  As  pointed  out  previou  y,  the 
calculations  for  j  =  0,  N  should  preceed  zi.  •'se  for 
intermediate  points. 

When  all  eight  prognostic  vari^vles  have  been 
calculated  at  time  level  k  +  1  for  all  j,  then  the 
position  coordinates  x,  y,  z  for  each  s  at  this  time 
level  can  be  evaluated  by  simple  numerical  quadrature 
using  the  relations 


bx/bs  =  (1  +  e)  a 

by/bs  =  (1  +  e)  p 

dz/ds  =  (1  +  e)  y 

Y  =  (1  -  c?  -  $>  )1/Z 


(156) 


with  the  conditions  x  =  y  =  z  =  0  at  s=0. 

The  whole  procedure  is  then  repeated  for  the  next 
time  level.  This  process  is  carried  out  for  some  pre¬ 
selected  number  of  iterations,  depending  upon  the  time 
span  of  interest  and  the  selected  value  of  At.  The  only 
limitation  on  At  and  as  in  this  method  is  that  of 
accuracy  and  resolution;  the  method  if  inherently  stable. 


112 


Accuracy  is  improved  in  principle  by  selection  of  small 
values  of  both  as  and  At.  However,  economy  of  computation 
plus  recognition  of  accumulation  of  round-off  errors 
implies  some  compromise  value  for  these  parameters.  It 
is  recommended  in  such  computations  that  as  be  taken  as 
about  100  feet  and  that  At  be  of  the  order  of  1/10 
second.  These  values  are  based  upon  the  considerations 
of  the  possible  modes  of  motion  treated  in  Chapter  VI 
via  linearized  perturbation  approach. 
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8.  Generalized  Upper-End  Condition 

In  the  preceding  discuss:*.^;-  it  was  assumed  that 

V  ,  V  ,  and  V  were  stipulated  as  a  function  of  time  at 
x  y  z 

the  upper  end  of  the  line.  A  more  general  condition  was 
implied  in  Chapter  II  (section  5),  in  which  the  dynamics 
of  the  surface  vessel  is  taken  into  account.  The  purpose 
of  this  section  is  to  show  how  such  conditions  can  be 
imposed  in  place  of  the  simple  upper  conditions  of  the 
preceding  section. 

Let  q.J  (i  =  1,  2,  3)  denote  the  three  Cartesian 
components  (x,  y,  z  respectively)  of  the  anomoly*  cf 
translational  velocity  of  the  center  of  mass  of  the 
surface  vessel,  and  let  denote  the  Cartesian  components 

of  the  anomoly*  of  rotational  velocity  of  the  vessel.  It 
will  be  assumed  that  the  x-axis  is  taken  such  as  to  be 
parallel  to  the  mean  orientation  of  the  longitudinal  axis 
of  the  vessel.  This  will  correspond  to  the  mean  orientation 
of  the  surface  current  for  normal  mooring  with  a  single  line. 

Then  in  the  notation  of  Section  3,  Chapter  II,  we 
require,  for  the  case  of  single  mooring  line,  that  along 
ds/dt  =0  at  j  =  N  ( s  =  L) : 


*  These  are  anomolies  from  the  velocities  which  the 
vessel  would  have  in  the  absence  of  the  mooring  line. 
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dR£ 

dt”  =  qi  ’  (157) 


d*i 

dt”  “  Qi  *  (158) 


dq.1 

dT"  =  t"  T  Ti^i  “  2  % 

-  a.a  R£  +  To  t±oM±]  ,  (159) 


do; 

dt”  =  ^  ^  Pi 

-  ai3  *i  +  Jt0/L^  *  (160) 


where  i  =  1,  2,  35  and 


™  A  a  A 

V  =  q*  +  Q»  x  P°  +  F(t)  .  (161) 

In  (159)  Tj_  are  the  Cartesian  components  of  the  unit 

a 

vector  T  (previously  designated  a,  g,  Y)«*  The  term 

A 

F(t)  is  a  prescribed  vector  function  of  time  denoting  the 
potential  velocity  of  the  point  of  attachment  of  the 
mooring  line  on  the  vessel,  in  the  absence  of  the  mooring 


*  The  reader  is  cautioned  that  there  is  a  potential 
ambiguity  in  the  semantics  here;  the  scalars  a,  g,  X 
have  a  meaning  quite  distinct  from  the  vector  components 
ai>  of  Section  5,  Chapter  II. 


115 


line.  In  terms  of  the  notation  of  Section  5,  Chapter  II, 

A  A  A  A 

F(t)  =  x+axpP.  The  terms  are  the  Cartesian 

components  of  a  vector  torque  which  depends  upon  T,  T 
and  the  angular  displacement  of  the  vessel  (see  Eq.  (31)  ). 

One  can  readily  write  finite  difference  counterparts 
of  relations  (157)  to  (160)  analogous  to  those  for  the 
prediction  of  § .  These  allow  the  evaluation  of  the 
vector  velocity  at  the  new  time  level  via  relation  (l6l). 
Having  the  velocity  components  at  the  new  time  level  for 
0  =  N  then  makes  the  system  equivalent  to  that  described 
in  the  preceding  section.  It  must  be  noted,  however,  that 
the  prediction  of  q£  and  fjf  depend  on  the  new  T  and 

a 

t  at  j  =  N.  Hence  a  successive  approximation  procedure 
is  required  similar  to  that  described  in  the  preceding 
section.  The  initial  conditions  on  the  components 
Ri’  *1*  and  must  be  prescribed.  Bear  in  mind, 

however,  that  these  variables  are  pertinent  only  to  the 
single  position  j  =  N. 

The  generalization  of  the  above  procedure  for  the 
case  of  multiple  mooring  lines  is  straightforward  following 
the  system  indicated  in  Section  5,  Chapter  II. 
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VI.  LINEARIZED  PERTURBATIONS  RELATIVE  TO 
A  COPLAN AR  EQUILIBRIUM  CONFIGURATION 

I.  Dynamic  Equilibrium  State 

A 

The  oceanic  current  U  in  general  consists  of: 

(1)  a  quasi-permanent  horizontal  flow  which  decreases 
slowly  with  increasing  depth;  (2)  horizontal  tidal  flow 
which  is  nearly  uniform  with  depth  but  changes  slowly 
in  direction  and  speed  with  time;  and  (3)  rapidly 
fluctuating  three-dimensional  motion  associated  with 
surface  gravity  waves.  The  latter  phenomena,  which  has 
a  spectral  band  width  from  about  0.3  to  6  rad/sec,  is 
confined  tc  a  relatively  thin  surface  layer.  As  a 
practical  approximation  we  will  suppose  that  the  wave 
motion,  while  directly  governing  the  motion  of  the  vessel, 
does  not  directly  influence  the  mooring  line.  In  effect, 
we  are  implying  inat  the  fluctuating  part  of  the  current 

A 

U  is  confined  to  a  layer  comparable  to  the  draft  of  the 
vessel  and  that  over  the  bulk  of  the  mooring  line  the 
current  can  be  considered  steady,  at  least  for  durations 
of  the  order  of  one  hour. 


On  the  other  hand  the  line  itself  is  not  static 
but  is  in  a  continual  fluctuating  state  of  motion  in 
response  to  the  motion  of  the  vessel  to  which  it  is 
secured.  Now  the  drag  force  is  non-linear ly  related  to 

A  A 

the  relative  motion  (U  -  Vn)  of  the  fluid  with  respect 
to  the  line.  Accordingly  the  average*  configuration 
of  the  line  (or  lines),  for  given  steady  current  U  and 
a  stochastically  stationary  time  sequence  of  ship  motion 
with  given  variance,  will  not  in  general  be  the  same  as 
the  steady  state  configuration  of  the  line  with  the  same 
current  U  but  with  the  ship  fixed  in  its  mean  position 
(i.e.,  no  wave  action).  The  effect  of  the  fluctuating 
motion  of  the  line,  even  though  its  mean  value  is  zero, 
is  to  produce  a  greater  net  drag  than  that  due  to  U 
alone.  This  can  be  quite  marked  if  the  amplitude  of 

a  a 

V  exceeds  thp  magnitude  of  U. 

Since  the  average  configuration  depends  upon  the 
perturbation  of  the  line  through  the  non-linear  drag 
force,  we  will  refer  to  this  configuration  as  the  dynamic 
equilibrium  state.  It  will  serve  as  the  reference  state 
in  the  definition  of  the  perturbation  variables  of  the 
system.  In  particular  we  will  denote  all  reference 

*  If  the  motions  of  the  line  are  stochastically 
stationary  then  time  averages  are  equivalent  to  expected 
values  in  the  probability  sense  (Wiener,  1950). 


This  can  be 


n  O 


XXO 


state  variables  by  a  bar  over  the  symbol, 
interpreted  as  the  average  value  of  that  variable  at 
material  point  s  for  an  ensemble  of  possible  con- 

A 

figurations  associated  with  the  current  U  and  a 
stipulated  ensemble  of  ship  motions. 

A 

In  regard  to  the  current  U,  we  will  suppose  that 
(below  the  vessel)  this  is  not  only  steady  but  also 
horizontal  and  coplanar  at  all  depths.  Accordingly,  if 
the  statistical  distributions  of  the  component  velocities 
of  the  line  have  zero  mean  value  and  are  symmetric  about 
their  mean,  then  we  should  expect  that  the  dynamic 
equilibrium  configuration  will  lie  in  a  vertical  plane 

A 

parallel  to  the  current  U.  Actually  this  also  assumes 
that  the  mean  force  exerted  on  the  vessel  by  the  combined 
action  of  currents,  waves,  and  wind  is  in  the  sam, 
direction  as  the  current.  We  will  take  the  x, z-plane 
to  coincide  with  the  plane  of  this  dynamic  equilibrium 
configuration . 

Specific  relations  for  dT/ds  and  Sg’/ds  for  the 
mean  state  are  given  in  the  following  section  in  which 
the  perturbation  equations  are  derived. 
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2.  Quasi-Linear  Perturbation  Equations 

The  equations  of  motion  expressed  in  a  natural 
coordinate  system  based  on  the  mean  configuration  of  the 
line  are  the  most  convenient  for  analysis  of  perturbations 
of  the  line.  These  can  be  derived  in  a  manner  similar 
to  that  employed  in  section  3  of  Chapter  V.  In  formulating 
these  relations  we  must  be  careful  to  distinguish  between 

a 

an  instantaneous  vector-  like  t  and  its  mean  value  for 
material  point  s.  Since  the  overbar  notation  is  comber- 

a 

some  for  vectors,  we  will  denote  the  average  vector  T 

A  A  A 

as  T(J  at  given  s  and  similarly  for  v  ana  The 

A  A 

mean  unit  vectors  v0  and  t0  lie  in  "the  x,z-plane 

A 

while  p0  is  parallel  to  the  y-axis. 

Let  V,’.  VJ,  V*  hereafter  denote  the  components 

A 

of  V  in  the  fixed  local  coordinate  system  defined  by 

AAA 

the  vectors  v0j  u*,,  t0  at  point  s.  With  this 
understanding  the  vector  acceleration  is  simply 


av  dv«  A  A  3V'  A 

St  =  ST  v°  +  St  +  ST  T°  * 


(162) 


The  vectors  v0  s  Pq  ,  t0  are  related  to  the  Cartesian 
reference  vectors  by  the  following  relations  in  analogy 
to  (88)  with  e  replaced  by  J  and  T  =  0: 


p  '  J  '  I 

I  -l  i 


v0  =  i  cos.  0  -  k  sin  0  , 


(163) 


=  i  sin  0  +  k  cos  e  . 


Now  let  vS  u’,  t’  denote  the  perturbations  of  the  vectors 

A  A  A  A  A 

v>  |ij  t  in  sense  that  v  =  v0  +  v’j  etc.  It  will  be 
understood  that  the  perturbations  are  small,  i.e., 

A  A  A 

|v’|  «  1,  ||i* |  «  1,  and  |T*|  «  1.  If  we  also  let 

0 1  and  i 1  represent  associated  perturbations  of  0 
and  i  where  { 0  *  J  «  1  and  |$*|  «  1  then  we  obtain 
the  following  approximations  based  on  the  use  of  (89): 


v  *  =  no  cos  e  *  1  -  t0  0  *  , 


ul  =  -  v0  cos  0  $  1  -  T  sin  0  $  1  ., 


(164) 


-  Vo  6  *  +  Mo  £*in  8  $  1 


Now  referring  to  (13)  we  note  that  the  term 

A  A 

T  dV  /dt  can  be  approximated  by  To  dV’/dt  neglecting 

T  w  T 

products  of  perturbation  quantities.  Moreover,  referring 

A  A 

to  (105)  we  see  that  the  quantity  U*dT/dt  is  of  second 
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order  if  U  is  of  order  V  or  less.  Thus  (15)  reduces 
approximately  to 


(m  +  raa>  ~  Ea  To  St" 

Now  consider  the  first  term,  on  the  right-hand  side  of 
(165).  Letting  T  =  T  +  T1  where  T*  «  T  and 

AAA 

T  =  To  +  t 1  we  obtain  after  neglecting  quantities  of 
second  order: 


(165) 


*0  )  +  (T  T  '  )  +  (T  '  T0  )  . 


(166) 


Now  from  (163) 


we  find  that 

*  a? 

T°  "5s  9 

0  , 


(167) 


Hence  using  (164)  and  (167),  relation  (166)  takes  the  form 


12? 
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0  =  *5s  ”  w  cos  6  »  (170c) 

a  a 

the  average  of  G  ♦  t0  being  zero.  These  are  the  relations 
defining  the  dynamic  equilibrium  configuration. 

a 

The  vector  G  was  defined  by  (106)  and  (107).  We 
note,  first  of  all,  that  this  simplifies  somewhat  with 

a 

dU/dt  =  0.  In  order  to  put  the  remaining  part  in  a 
tractable  form  we  will  take 

Dn  -  Vn  i  (0  cos  8  -  Vv')  ;o  +  (-  y)  .  (171) 

A  A  A 

Furthermore,  the  term  r  t  x  (Un  -  Vn)  will  be  approxi- 
*  A  . 

mated  by  T  j0  x  Vo  U  cos  q.  This  implies  that  “strumming” 
of  the  line  by  vortex  shedding  is  allowed  only  in  association 
with  the  steady  fl  "'w  U,  while  that  associated  with  the 
oscillatory  motion  of  the  line  is  neglected*.  Accoraingly 
we  have 

G  ‘  v0  =  f[  (l-  -os  0  -  Vv*)2  +  (Vu*)2]1/2 

x  (U  cos  ?  -  V  *)  (172a) 

*indeed  the  phenomenon  of  vortex  shedding  is  some¬ 
what  confused  in  respect  to  oscillatory  relative  flow. 


U 
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Ho 


=  f[  (U  cos  0  - 


V  •) 
v 


+  V 


,)2]l/2 


+  T)  u  cos  e  cos  (u,  t  +  5)], 


G  *  To  =  0  , 


(372b) 

(172c) 


where 


f  =  (p/2)  Cd  D  (1  +  *)  . 

The  latter  coefficient  can  be  expressed  in  the  alternative 
form 


,  2  cd  nd 

f  "  TT  D  (173) 

where  md  is  the  displaced  mass  of  water  per  unit 
material  length. 

The  residual  equations  of  motion  are  obtained  by 
subtracting  the  corresponding  relations  in  (169)  and 
(170).  For  convenience  we  will  use  the  abbreviated 


notation 
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g  *  =  sin  e  *  »  , 


(174) 


(175) 


The  term  p '  in  fact  represent  the  perturbation  of  the 

* 

y  component  of  T  as  defined  by  (132).  The  resulting 
dynamic  perturbation  relations  are: 


m* 


dVv* 

ot 


St 


+  T> 


Se 

•Ss 


+  G  « 
v 


(176a) 


m 


-zr 


ar  p» 

“3s 


+  g, 


(176b) 


m 


dV  « 

T 

“5t“ 


ii1.  a.  *  * 

^  0  1  ^s 


(176c) 


where  m*  =  m  +  ma.  The  relations  (176a,  b)  are  only 
quasi-linear  since  the  terms  G^’  and  G  ’  are,  in 
general,  non-linear  functions  of  the  normal  components 
of  velocity  of  the  line.  All  other  terms  are  linear 
functions  of  the  perturbation  variables. 

To  these  relations  we  add  the  following  linearized 
versions  of  the  kinematical  relations  (93)  to  (95): 
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^  +  £)  It  = 

(!  +  7)  ||~  = 


dVv* 

-7T-+  VT 


(177a) 

(177b) 


(177c) 


Finally  the  linearized  versions  of  the  load- strain  relations 
(69)  xo  (71)  are  taken  as  follows 

Ti  It'  -  (*'  -  5')  -  (ITS) 

T>  =  Y0  e'  -  Kx  5>  ,  (179) 

where  e0  Is  regarded  as  a  constant  ( the  maximum  tension 
being  assumed  such  that  e0  >  F  (Tmax))-  In  the  above 
relations  e  1  =  e  -  7,  §  ’  =  §  -  7  where  7  =  7  -  e0 

and 


Y0  =  K  +  bT  +  biqs,  (180) 

which  represents  the  dynamic  spring  coefficient  for  the 
mean  conditions  at  s. 
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i  c 

3.  Associated  Perturbation  Energy  Equation 

If  we  multiply  (176cLy.  b,  c)  by  V  *,  V  ‘  and 

v  |i  > 

VT '  respectively,  then  multiply  (177a,  b,  c)  by  T  0», 
T  p‘,  and  T*,  then  multiply  (178)  by  -iq  (e>  -  §*)/Tl 
and  add  the  restating  relations  we  fini  the  following 
energy  relation  after  using  (179)  and  collecting  terms: 


It  {  £-5  m'  (V*  t  Vs)  +  5»  VT'»] 

* 

+  [|  Y0  ««  -  iq  e*  §•  +  |  iq  g«] 

+  [|  (i  +  7)  t  ce«  +  p**)]  i 

=  Is-  [  t  (e«  vv«  +  p*  v^)  +  t *  vT>] 

iq 

+  Gv>  Vv*  +  G^»  V^*  -  —  (c  *  -  S’)2  , 


(181) 


Comparing  this  energy  relation  with  (85),  it  appears 
that  we  have  introduced  a  new  form  of  energy  term  involving 
(0i2  +  p»s).  Actually  the  last  two  sets  of  terms  in 
brackets  on  the  left  sidj  of  (l8l)  represent  dynamic 
elastic  energy  associated  with  the  perturbation.  Only 
part  of  the  anomaly  of  elongat tonal  strain  is  represented 
in  the  term  e>;  it  is  that  pa.-rt  related  to  motion 
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parallel  to  the  reference  configuration.  If  transverse 
waves  exist  relative  to  this  reference,  then  second  order 
elongational  strain  occurs  which  is  proportional  to 
(0i2  +  3«a).  Thus  the  term  (1/2)  (1  +  7)  T  (6  12  +  p'2) 
is  the  elastic  energy  associated  with  transverse  perturbation, 
while  the  term 

[•|  Y0  e’2  -  Kj.  e'  §*  +  -|  K*  §'2] 

is  the  elastic  energy  associated  with  the  longitudinal 
perturbations.  Note  that  this  will  always  be  positive 
as  long  as  Y0  >  iq  which  is  certainly  the  case  (see 
Chapter  III). 

We  note  that  terms  involving  the  curvature  of  the 
mean  configuration  (djf/ds)  which  were  present  in  (176a,  c) 
(177a,  c)  have  no  contribution  in  the  above  perturbation 
energy  equation.  These  terms  lead  to  coupling  of,  and 
hence  energy  exchange  between,  the  transverse  and 
longitudinal  perturbations  but  do  not  account  for  any 
loss  or  gain  for  the  total  dynamic  energy  budget.  On 
the  other  hand,  the  last  three  terms  involving  drag  effect 
and  hysteresis  indeed  do  significantly  influence  the 
total  energy  budget. 
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If  the  perturbations  represent  a  stochastically 
stationary  process  then  the  average  total  energy  is 
constant  in  time  (or  at  least  only  slowly  varying) 
and  hence  the  average  of  the  terms  on  the  right  side 
of  (l8l)  should  produce  a  balance.  The  drag  terms  are 
particularly  significant  in  this  balance  of  energy  loss 
and  energy  gain  per  unit  time.  This  will  be  considered 
in  a  statistical  sense  in  the  subsequent  section. 

v 

4.  Statistical  Considerations  of  the  Drag  Force 

In  order  to  arrive  at  a  rational  approximation 
of  the  highly  npn-linear  relations  for  the  drag  force, 
we  will  consider  certain  critical  statistics  associated 
with  this  force  for  a  simple,  yet  realistic,  stochastic 
model  of  the  line  motion.  The  motions  of  the  sea 
surface  and  of  the  surface  vessel,  in  response  to  the 
surface  waves,  can  be  considered  to  represent  a  Gaussian 
process  in  the  first  approximation  (St.  Denis  and  Pierson, 
1955).  We  will  suppose,  moreover,  that  V  1  and  Va* 
are  also  of  Gaussian  character.  The  principal  additional 
assumptions  are  that  these  variables  are  statistically 
independent  and  that  they  have  the  same  variance  a3 
and  zero  mean  value.  In  general  a3  will  vary  with 
position  s  on  the  line.  The  assumption  of  equal 


130 


is  perhaps  the  most  stringent  of*  the  assumptions, 
but  considerably  simplifies  the  analysis  of  the  statistics 
related  to  the  drag  force. 

For  simplicity  of  notation  in  the  present  analysis 
we  will  let 


=  x, 

a  x*  (182) 

s  y  . 

Since  V  1  is  assumed  to  be  Gaussian  with  zero  mean  and 
variance  a2 ,  then  x  is  Gaussian  with  mean  x  and 
variance  a2.  Moreover  the  joint  probability  density  of 
x  and  y,  in  view  of  the  foregoing  assumptions  is 


U  cos  0 


U  cos  0  -  V  * 
v 


V  1 
u 


F(x,  y) 


1  ~[(x  -  x)2  +  y2]/  2  a2 

2m  a2  e 


(183) 


In  the  above  notation,  the  component  of  G  in  the 

A 

direction  Vo  is 


G  .  0,  =  t  (x>  +  y=)-1/2 


;i84) 


X  . 
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132 


Q  =  x/a  =  U  cos  e/o  , 


2 

i 


(188)  1 

3 


and  I0(z),  Ij (z)  are  modified  Bessel  functions  of 
order  indicated.  Prom  the  properties  of  the  latter 
functions  it  is  readily  shown  that  the  mean  drag  force 
has  the  following  limiting  forms 


r  _ 

f  (U  cos  0)3,  for  U  cos  q  »  a 

* 


i  M 


f  a  U  cos  0,  for  U  cos  q  « 


(189) 


In  the  latter  case,  it  is  evident  that  the  standard 
deviation  (a)  of  V  *  or  V  ’  has  a  very  definite 

V  ji 

control  on  the  curvature  of  the  equilibrium  configuation 
through  relation  (170a). 

Another  quantity  of  considerable  concern  is  the 
mean  rate  of  energy  dissipation  by  the  drag  force.  This 
is  given  by 


3D  =  f[(U  cos  0  -  V  ')*  +  (V  t)*]1/2  [(v  »)a+(V  ’)2] 


4 


(190) 
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In  terms  of  the  notation  (182)  this  can  he  written  in  the 
form 


B  =  f  r  [  (x  -  x)*  +  y»] 

=  f  (r3+  x8  r  -  2  x  rx  ),  (191) 

1/2 

where  r  =  (x8  +  y8)  '  .  Again  using  (183)  to  ev  equate 

the  averages  in  the  above  relation  we  find 

I  =  f  &  e'CaA  [C»  +  3)  I„(C2A) 

+  C2  Ii  (c=A)l  .  (192) 

where  c,  is  as  defined  in  (188) . 

We  can  now  define  an  effective  dynamic  drag  coefficient 
f  v  such  that 


3D  =  f  v  (Vv18  +  Vj»2)  .  (193) 

Novr  =  V^'8  —  °s  i*1  the  foregoing  analysis 

and  hence 


13^ 

V  -  f  B  -C2A  [(C3  +  3)  X0(cb/4) 

+  Ca  McaA)].  (194) 

This  has  the  limiting  values 

for  U  cos  *9  »  a 

for  U  cos  ?  «  a  .  (195) 

A  plot  of  v/a  versus  C  is  shown  in  Fig.  18.  This 
indicates  a  very  rapid  transition  from  one  limiting 
form  to  the  other,  centered  at  q  =  =  1.88 

(i.e.,  U  cos  7  =  1.88  a).  The  behavior  of  the  function 

A  A 

G  •  jQ/f  a2  £  is  indistinguishable  from  v/a  for  all 
practical  purposes. 

Based  upon  the  above  analysis  we  will  adopt  the 

A 

following  approximation  of  the  components  of  G  : 

G  •  v0  ==  f  \  (U  cos  7  -  V  *)  ,  (196a) 

5  ’  Ho  =  ■  f  v  V  '  +  t]  f  v  U  cos  7  cos  (u)  t  +  6)  .  (196b) 
0  H  s 


r 


v  =--  < 


U  cos  Q, 


3  '<  it  „ 

2  V  2  0  ' 
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Thus 

G  *  =  -  f  v  V  »  (197) 

V  V 

where  v  is  at  most  a  function  of  s.  This  approximation 
with  v  as  given  by  (19^)  leads  to  the  appropriate  mean 
drag  and  also  the  appropriate  rate  of  energy  dissipation 
by  drag  force. 

5.  Perturbations  of  a  Nylon.-  Mooring  Line  about  a 
Straight  Equilibrium  Configuration 
In  the  case  of  nylon  line,  the  net  weight  per  unit 
length  in  water  is  negligible  and  can  be  ignored  for  all 
practical  purposes.  From  (170c)  this  implies  that  T 
is  essentially  uniform  over  the  whole  line.  Now  from 
(170a)  we  find  that  the  curvature  of  the  equilibrium 
configuration  is  given  by  the  mean  drag  force  per  unit 
length  divided  by  T.  If  the  latter  is  large,  the  mean 
configuration  will  tend  to  be  nearly  straight.  In  this 
limiting  case  of  negligible  cie/os,  there  is  no  coupling 
between  transverse  and  longitudinal  oscillations  of  the 
line.  Hence,  these  modes  of  oscillation  can  be  analysed 
independently.  We  will  consider  this  situation  in  the 
present  section. 
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In  the  considerations  of  the  transverse  oscillations 

we  will  examine  two  limiting  cases  in  respect  to  the  drag 

force.  These  are  the  cases  where  U  cos  J  >  5  a  and 

U  cos  0<O.5o  (see  Pig.  18).  In  the  first  case  the 

term  v  is  independent  of  a  and  the  problem  is  entirely 

linear.  In  the  second  case  v  is  directly  proportional 

to  a  which  in  turn  is  a  function  of  s  dependent  upon 

the  statistical  behavior  of  the  solution  for  V  '. 

v 

Clearly  the  latter  case  is  a  non-linear  problem  which  must 
be  treated  quite  differently  from  the  first. 

Case  A.  Transverse  Modes  (U  cos  6~  >  5c):  In 
this  case  v  i  U  cos  0*,  which  will  be  treated  as  a 
constant,  at  least  over  that  portion  of  the  line  where 
the  transverse  oscillations  are  significant.  For  realistic 
values  of  f  v,  we  will  find  that  the  transverse  modes 
are  not  significant  at  depths  greater  than  about  1000  ft 
below  the  vessel  and  hence  the  approximation  of  uniform 
v  seems  justifiable  in  respect  to  the  analysis  of  these 
perturbations . 

With  uniform  T  and  zero  mean  curvature,  relations 
(176a)  and  (177a)  reduce  to 


m 


'sf-  =  T  fr  -  f  v  V  * 


(1  +  -)  m:  ,  **  . 

St  3s 


Eliminating  0 1  between  these  relations  yields 


S2V  1 

V 

-zt?~ 


&V  •  ~  „  s  v  * 

c  2  _  V _ UL 

!  3  s  m  ’  3t 


wnere 


Ci  =  (T/m>  (1  +  e))1/2 


(198) 

(199) 


(200) 


(201) 


This  has  elementary  solutions  of  the  form 

V  1  =  A  e  1(U)t  +  ks)  ,  (202) 

v 

where 

k2  =  (uj2  -  i  f  v  uu/m 1 ) /Cj 2  .  (203) 


If  uj  is  regarded  as  real  then  k  has  the  roots 
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k  =  ±  (h  -  i  u)  (204) 

where 

x  =  |k|  cos  a/2 

u  =  |  kj  sin  a/2  (205) 


in  which 

|k|  *  «,[  1  +  (f  v/m«  u))2]1/4/Cx  (206) 

and 

a  =  tan"1  (f  v/m*  u>)  -  (207) 

Note  that  *  is  an  odd  function  of  u,  while  jj.  is  even. 

A  more  general  solution  is  obtained  by  forming  a 
linear  combination  of  solutions  of  the  type  (202)  for  the 
two  admissible  values  of  k  for  given  ou.  Such  a 
solution  satisfying  the  condition  1  =  0  at  s  =  0 

for  all  t  is 


1*10 


V  '  =  A  [euS  cos  (ujt  +  Hs  +  p) 

-  e”^s  cos  (uut  -  xs  +  p)  ]  t  (208) 


in  which  we  have  allowed  for  an  arbitrary  phase  angle  p. 
The  first  term  in  the  brackets  represents  a  progressive 
wave  traveling  downwards  along  the  mooring  line,  its 
amp?itude  decreasing  exponentially  in  the  direction  of 
travel.  The  second  term  is  a  reflected  wave  traveling 
upwards  along  the  line  and  attenuated  exponentially  in  the 
direction  of  travel,  if  the  attenuation  is  sufficiently 
large,  or  the  line  sufficiently  long,  then  the  reflected 
wave  may  not  even  exist.  The  amplitude  of  V\> ’  at 
the  position  of  the  vessel  (s  =  L)  can  be  shown  to  be 
given  by 


1  /p 

A  (2  cosh  2  u  L  -  2  cos  2  *  L;'  . 


(209) 


UL 


If 


L 


» 


1 


then  this  reduces  simply  to  A  e 


14 1 


Using  (173)  we  find  that 

f  v  h  Cd  md  )  v 

mru)  \7r  m1  /  (ij  D  * 


(210) 


For  nylon  line,  the  factor  in  parentheses  is  approxix  itely 
1/tt.  Suppose  v  is  only  1  ft /sec  and  u)  =  2r/10  rad/sec 

corresponding  to  a  10  second  period  excitation.  Then 
for  a  2  inch  diameter  line 


f  v  . 
m*oo 


3  . 


Using  relations  (205)  to  (207)  and  taking  =  700  ft/sec 
(corresponding  to  T  equal  to  about  half  of  the  ultimate 
strength)  yields 

x  =  1.3  x  10-3 /ft 
m  =  0.9  x  10-3/ft  . 

The  wave  length  is  2f/x  or  about  4800  ft;  on  the  other 
hand  an  e-fold  attentuation  of  amplitude  occurs  in  an 
arc  length  of  1/p.  or  about  1100  ft. 


The  attenuation  of 
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amplitude  in  one  wave  length  is  about  99  per  cent.  Thus 
for  all  practical  purposes,  the  transverse  waves  are 
critically  damped  3n  one  wave  length. 

The  above  example  serves  to  demonstrate  the  extreme 
importance  of  drag  effect  on  the  transverse  oscillations. 
Typical  values  of  fv/m'Jiul  exceed  unity.  In  the 
limiting  case  of  very  large  values  of  this  parameter  we 
find 


M  ~  M 


f  V.  (-1  +  e) 
T 


'I  1/2 


(211) 


for  f  v/m’  J  u)|  »  1.  This  situation  corresponds  to 
negligible  inertial  effect  (the  term  m>  dV^'/dt  being 
negligible  in  (198)). 

The  solution  for  V  '  given  by  (208)  is  valid 
provided  that  the  excitation  at  the  upper  end  is  periodic 
in  time  (with  p>erlod  2ir/m) .  A  more  realistic  condition 
at  the  upper  end  would  allow  for  the  representation  of 
V  «  as  a  general  time  sequence  with  stipulated  statistical 
properties.  We  will  employ  here  the  Fourier-Sfielt jes 
integral  representation  as  follows: 

09 

V  <L>t)  =  J  («.) 

—03 


3 


(212) 
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where  dFx  (uo)  is  a  stochastic  complex  variable  having 
the  following  statistical  property 

dF7*U7  dF.  („,»)  =  6,(d)  -  u>')  (w)  dw  ,  (213) 

where  6*(x)  is  the  Dirac  delta  function,  dPj *  is  the 
complex  conjugate  of  dF! ,  and  (uu)  is  the  variance 
spectrum  associated  with  * (L,  t) .  The  variance  of  * 
at  the  upper  end  of  the  line  is 

CD 

°i2  (I*)  =  J  *xU)  da.  •  (214) 

—oo 

The  function  dFj (w)  also  has  the  following  symmetry 
property 


dFj  (-a)  =  dFx*(»)  ,  (215) 

which  assures  that  the  representation  (212)  leads  to 

real  values  of  V  *  . 

v 

For  other  positions  on  the  line,  the  generalization 
of  (202),  taking  into  account  that  the  waves  are  essentially 
attenuated  completely  before  they  reach  the  anchor 
(liL  »  1),  is  as  follows 


ihk 


oo 

Vv'(s,t)  =  f  e^s  "  L)  ei^U)t  +  dP^oj)  .  (216) 

—  oo 


This  reduces  to  (212)  at  s  =  L.  From  (216)  and  (213) 
we  find  that  the  variance  of  V  *  at  any  position  s 
is  given  by 


00 

v( S)  =  J  e2^s  ~  L)  ^(«))  do,  ,  (217) 

.CO 


where  u  is  nearly  proportional  to  jtuj1^2,  see  (211). 

In  the  case  of  transverse  perturbations  in  the 
lateral  (y)  direction,  the  appropriate  equations  (under 
the  conditions  of  constant  T)  are 


dV  ' 
— U 
6t 


v  V  1 


+  f  v  P 


(218) 


(!  +  e)  ||- 


av  * 

u 

'6's  * 


(219) 


where  P  represents  a  random  function  associated  with 
the  vortex  shedding.  We  can  represent  this  as  follows 
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P 


r  e1”*  dft(»)  , 


where 


(220) 


dQ*(uo)  dQ(tu') 


6  '  (u)  -  a)  ’)  <fs(uu)  duu 


(221) 


and 


00 

J  ^s(u>)  dm  -  (n  U  cos  8)a/2  . 
—00 


(222) 


If  the  vortex  shedding  were  exactly  periodic  then  its 
spectral  function  would  be 

ijrs(u))  =  -  ujg)  (ri  U  cos  e)2/4  .  (223) 

However,  the  above  representation  for  the  vortex  shedding 
will  allow  some  smearing  over  a  finite  band  width. 

Eliminating  p1  between  (218)  and  (219)  and  using 
(220)  yields  the  non-homogeneous  relation 


3.46 


1  «  c*V  * 

“3^2 '  C1  2  "5i2  +  iT  ’5tU" 


fv 


J  i  0)  elu)t  d  Q(u)) 


(224) 


An  appropriate  solution  is 


V  » 
U 


(s,t)  =  r  e^8  ■  L)  e1  (a,t  +  xs>  dP,(o)) 


+  j  (1  +  i  ra’uu/fv)"1  [  1  -  e^<s  "  L>  +  1X3 


-e'»S  '  lxS]  e1^  dQ(uj)  , 


(225) 


where 


dFa(u>*)  =  6(a)  -  uu1)  (m)  dm  .  (226) 


14? 


Moreover,  dFs  and  dQ  are  regarded  as  statistically 
independent  for  all  w;  this  assumes  that  the  input  by 
the  surface  vessel  is  completely  indepdent  of  the  vortex 
shedding  and  vice  versa.  Considering  chat  »  1, 
the  solution  (225)  reduces  essentially  to  zero  at  s  =  0 
while  at  the  upper  end  we  have 

00 

Y  (I,t)  =  f  e  imt  U)  ,  (227) 


which  is  governed  by  the  motion  of  the  vessel. 

The  variance  of  ‘  at  any  s  is  given  approximately 
by  the  relation 
00 

022(s)  =  |  e2  «<s  -  L)  M  dw 


00 

■  +  J  [1  +  (m«  uo/f v ) 3 3 ”*1  {  1  -  2[(eu  (s  "  L) 

—00 

+  e  ^S]  cos  xs  +  e 

+  e~  2^s  |  (uj)  dm  .  (228) 
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Again  considering  that  »  1,  the  above  relation  gives 
essentially  zero  value  at  s  =  0  and  at  the  upper  end 
we  get  simply 


V  (L) 


J  fe  (m)  duo 

—OO 


(229) 


Case  B.  Transverse  Modes  (U  cos  ~e  <  0.5a): 
Relations  (200)  and  (224)  are  still  valid  for  this  case, 
however,  we  now  have 


at  least  under  the  condition  that  the  variance  of  V  *  and 

v 

V  1  have  the  common  value  a  at  position  s.  We  have 
u 

seen  in  the  previous  case  that  for  intermediate  and  large 
depths  below  the  vessel,  the  oscillations  in  the  lateral 
direction  may  actually  dominate.  Hence  the  present  analvsis 
is  limited  to  the  upper  reaches  of  the  line  where  Oj2 
remains  comparable  to  a22, 
equal  at  the  vessel. 


assuming  that  these  are  nearly 
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As  before  we  will  employ  the  representation  (212) 
for  the  input  sequence  V^(L,  t).  For  other  s  we  then 
propose  a  solution  of  the  form 

00 

V  »  (s,t)  =  f  M  (s,  w)  el(u,t  +  N(s'*>)  dFx  (u))  ,  (231) 

V  v 


where  functions  M(s,U))  and  N(s,ty)  remain  to  be  determined. 
In  order  for  this  form  to  be  consistent  with  (212)  we  re¬ 
quire  +hat 


M  (L,  oj)  =  1,  (232a) 

N  (L,  cu)  =  0.  (232b) 

It  will  be  observed  that  the  solution  (2l6)  for  case  A 
is  simply  a  special  case  of  (231)  in  which  M  is  an 
exponential  function  of  s  and  the  phase  function  N 
is  linear  in  s. 

Using  (231)  and  the  statistical  property  (213) 

we  find  that  the  variance  of  V  ’  is  given  by 

v 

CO 

o12( s)  =  J  IP  (s,  id)  (w)  dm 


9 


(233) 
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where  the  spectral  function  ^  (w)  is  presumed  to  be 
stipulated. 

Relation  (231)  will  satisfy  (200)  provided  that 

-  uu2M  =  Ci2  (M"  +  i  2  M‘  N*  +  i  M  N" 

-  M  (N’)2]  -  i  f  v  o»  M/m«  ,  (234) 


where  M*,  N’  imply  derivatives  with  respect  to  s.  The 
functions  M  and  N  are  considered  real.  Hence  upon 
separating  real  and  imaginary  parts  of  (234)  we  obtain 
the  following  two  relations: 

o  wll  .  1/2 

»'  -  ±{(^)  +  5  }  (?35) 

and 

(M3  N * )  '  =  M3  f  v  uj/m’  C  2  .  (236) 

1 

Moreover  from  (230)  and  (233)  we  have 

o  / - r  00  ,1/2 

=  2  vilJ  M(s'w)  tx  (®)  da)  j 


V 


—00 


(237) 
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The  set  of  equations  (235)  to  (237)  represents  a 
highly  non-linear,  integro-dj ff erential  system  for  the 
unknown  functions  M  and  N.  An  exact  general  solution 
is  not  known.  However  it  is  possible  to  obtain  approximate 
solutions  for  two  limiting  cases.  These  are  based  upon 
the  magnitude  of  the  non-dimensional  parameter  f v/ujm  ’ . 

For  fv/tum’  «  1,  then  inertia  dominates  over  drag 
effect  and  we  should  expect  that  M  will  be  a  slowly 
varying  function  of  s.  This  implies  that  (uj/Cj  )2 
dominates  over  the  term  M”/M  in  relation  (235).  Thus 
N1  =  +  uu/C!  as  a  limiting  case.  Using  this  in  (236) 
leads  to 

(M2)’  =  +  M2  f  v  /ml  Cx  ,  (238) 

whi-'  h  does  not  involve  u>  explicitly.  Thus  M  is 
independent  of  ^  for  this  case  and  consequently  (237) 
and  (214)  yield 

v  =  1  \%  ®  (T-)  M(s)-  (239) 

Using  (238)  and  (239)  leads  to  the  differential  equation 

M"3  (K2)«  =  |  \p[  f  a(L)/ra»Cx  , 


(240) 
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in  which  our  choice  of  sign  confines  attention  to  waves 
propagating  downwards  along  the  line.  This  assumes  that 
the  line  is  sufficiently  long  that  the  transverse  modes  are 
completely  damped  before  reaching  the  anchor  point  (a 
rather  stringent  assumption  for  this  case)  A  solution 
of  (240)  satisfying  condition  (232a)  is 

M  =  [1  +  *,  (L  -  s)]"1  ,  (241) 


where 


*i  »  £  t  oflO/a'C,  .  (242) 

Finally,  from  (233),  (214),  and  the  above  result  we  get 
for  any  ^  («j)  : 

°2  (s)  =  [1  +  Ax(L  -  s)]"2  o»  (L).  (243) 

Note  that  the  attenuation  coefficient  %x  is  proportional 
to  a(L)  according  to  (242). 

For  the  case  f  v/uj  m1  »  1,  then  drag  effect 
dominates  over  inertia  and  hence  we  expect  in  this  case 

p 

that  M"/M  greatly  exceeds  (uj/C1  )  .  This  is  the  normal 
situation  which  would  be  expected  for  nylon  line.  This 
implies  from  (235)  that 
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N«  =  (M"/M)1//2  ,  (244) 

where  our  attention  is  again  confined  to  waves  propagating 
downwards  along  the  line  (and  attenuated  before  reaching 
the  anchor  point) .  As  a  generalization  of  the  solution 
(241)  we  will  consider  in  this  case  that 

M  =  (1  +  *,  (L  -  s)j-n  ,  (245) 

where  -ft2  and  n  are  to  be  determined  for  this  case. 
Relation  (244)  gives  the  following  associated  N* 
function 

N*  =  [n  (n  +  1)]1/2  *3  (1  +  *a(L  -  s)  )-1.  (246) 

Substituting  the  last  two  relations  into  (256)  and  using 
(237)  leads  to 

[n  (n  +  1)]1/2  (2n  +1)  V  U  +  (L  -  s)]"2 

00 

=  1  Vi  (!■  I»|A'  c;»  I  J  [l 

—00 


+  (L  -  s)] 


(247) 
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This  relation  can  be  satisfied  for  all  s  if  n  =  2, 
but  only  if  ^  (w)  is  a  very  narrow  spectrum  centered 
at  some  frequency  uu0.  As  a  rough  approximation  for  the 
case  of  a  realistic  spectrum  ^  (w)  we  will  take  w 
in  the  coefficient  on  the  right-hand  side  of  (247)  as  the 
mean  value,  u^,  defined  by 

00 

«%  =  '<F(l)  J  M  ti  (®)  d{U  •  (248) 

“CO 

Then  from  (247)  with  n  =  2 

**  ={lo(?5>1/2  JVS  ■  (249) 

The  phase  function  can  be  obtained  from  (246)  as 

N  (s)  =  ]/6  m  (1  *  (L  -  s)  ).  (250) 

Finally,  the  variance  of  V  ’  at  any  oosition  s  is 

v 

given  approxi_.  .tely  by 

a2  (s)  =  [1  +  (L  -  s)]“4  a*  (L)  (251) 


for  this  case,  where  -fts  is  given  by  (249). 
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The  results  for  the  approximate  relations  for  a  (s) 
for  transverse  waves  are  summarized  in  Table  3.  All 
cases  assume  that  the  line  is  sufficiently  long  that  the 
variance  is  attenuated  to  a  negligible  value  at  the 
anchor  point.  Also  the  approximation  is  made  in  each  case 
that  the  attenuation  rate  is  based  on  the  mean  frequency 
cum,  characterising  the  spectrum  of  motion  at  the  top 
of  the  line.  The  term  r>  defined  in  Table  3  is  a 
non-dimensional  quantity  in  which  v  is  either  U  cos  e" 
or  1.88a,  depending  on  the  case.  The  usual  conditions 
in  application,  at  least  for  the  upper  reaches  of  the 
line  are  for  r  -»  o  >  2  U  cos  J.  An  example  plot 
of  a  vs  (L  -  s)  is  given  in  Fig.  19  for  the  following 
conditions: 


D  =  1  inch  ,  iUjjj  =  2ir/10  rad/sec  , 

a  (L)  =  3  ft/sec. 
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Any  r 


r  « 


r  ^ 


TABLE  3 

APPROXIMATE  STATISTICS  FOR  TRANSVERSE  MODES 
IN  THE  EQUILIBRIUM  PLANE 


U  cos  e  <  0.5  a(s) 
v  =  1.88  a(s) 

|  a(s)  =  e~^  "^a(L)  a(s)  =  [1  +  *(L  -s)] 

! 

■ 

i  |  u  =  r  ! 

i  • 

i 
» 

i 

I _ _ _  _  j 

i  i 

l  !  U  =  0.707  r1/2  I  A  =  0-286  rfl A1^2  ^ 

'  ‘  |  "  ' 

. . . i  n:2  . 

*  r  H  f  v/m’  0Jm 


*  -  |  r  (l)  “>m/c, 

n  =  1 


U  cos  0  >  5  a(s) 
v  =  U  cos  o’ 


_1./T 


a(L) 


(ft/  sec) 


1^8 


where  U  cos  e  is  assumed  to  be  less  than  a/2  at  all 
(L  -  s).  This  requires  virtually  zero  value  of  U  for 
depths  greater  than  3000  ft  if  the  results  of  Fig.  20 
are  to  be  valid 


Case  C.  Longitudinal  Modes:  Relations  (17oC)  and 
(177c)  under  the  condition  of  zero  mean  curvature  reduce  to 


m 


av  * 

T 

aT« 

at 

=  TTs~ 

Se« 

av  * 

_  _ i— 

at  as 

(232) 

(253) 


which  together  with  relations  (178)  and  (179)  govern  the 
longitudinal  waves,  in  the  absence  of  mean  curvature. 
Also  for  T  uniform  then  Y0  is  independent  of  s. 

These  relations  lead  to  the  following  differential 
equations  for  '  : 


a 

"St 


f  8=  vT  '* 

{  — 


Y„  d=V' 
m  Ss2  - 


i  ffV 

T1  1  St2 


(y0  -  iq ) 


- 1 =  0. 

3s2  J 


/’PRZl^ 
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A  solution  of  this  relation  is 


V  >  =  A  +  ks)  , 


where 


k3 


For  most  applications  Y0  greatly  exceeds  iq  . 
this  condition,  k  has  the  approximate  roots 


k 


+  [  (uu/Cs  )  -  i  v]  > 


where 


v  =  (u>/C2)(iq/2Ye)  ojtiAI  +  Ur*)*] 


(  -1  1 
C8  =  |[Y0  -  iq  (1  +  (u,TX)2)  x3/m  j 


Under 


and 


i6o 


The  attenuation  coefficient  v  is  a  small  fraction  of  the 
wave  number  for  these  modes,  hence  we  must  definitely  allow 
for  reflection  of  the  longitudinal  waves  at  the  anchor 
point  so  as  to  assure  V  *  (o,  t)  =  0. 

The  generalized  solutions  for  an  ensembe  of  possible 
1  at  L  is  accordingly 

Vt*  (s,t)  =  J  [evS  e  1  u)(t  +  S/Cs) 


-  e 


-vs  i  uj(t  -  s/Cs 


)]B_1/2 


dPs(e>),  (260) 


where 


d  F3*(w)  dF3  (u>«)  =  6’(u>  -  U)J)  <ra(«)  (26l) 


and 


B  =  4[  sinh2  (VL)  +  sin3  (ujL/Ca)].  (262) 

The  function  '^(ou)  corresponds  to  the  variance  spectrum 
of  V  *  at  s  =  L,  its  integral  over  all  w  being  the 
variance  of  V  ’  at  the  vessel.  The  variance  of  V  " 

T  T 

for  any  position  s  is  given  by 


(iw)  d(U  .(263) 


(Trp 


sinh3  (ys)  +  sina  (ms/Ca  ) 
sinh2  (vL)  +  sin2  (uuL/Ca ) 


Clearly  this  vanishes  at  s  =  0  as  required  by  the  lower 
boundary  condition. 

From  (252)  and  (260)  it  follows  that 


T* 


[evS  e1  +  S/Ca) 


+  e"vS  e1  "  S/Cs)]  B ~X/2  d  F3  (oj)  .  (264) 

Thus  the  variance  of  T*  using  the  statistical  property 
(26l)  is,  at  any  position  s  , 


(T  *  )2  = 


_ UJ2 _  (  sinh2(ys)  +  COS2 (mS/Cg )  ? 

['(uj/C2)2  +  v2]  (  sinh2  (vL)  +  sin^L/c"  )  j 


V3  (u))  dli)  •  (265) 


Since  v  «  m/C2  and  yL  «  1  for  the  normal  range  of 
frequencies  (even  with  L  ~  104  ft),  the  above  relation 
can  be  approximated  by 
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wn 


2  =  m2 


r  cag  c°s2 u  m  do, .  {266) 

J  r  /..T .^2  !  o  -i  r\2  I  ...T  /r*  \2 


[  (VL)2  +  sin2  (u)L/Cs  ) 


As  a  check,  we  note  that  if  o3  (&)  is  confined  xo 
very  small  frequencies,  implying  slow  rates  of  stretching, 
then  the  above  relation  with  the  use  of  (259)  reduces  to 
approximately 


WY  =  (Y0  -  Kx  )2 


I 


iu 


,-2 


$3  (uu)  d(jj  , 


(267) 


for  any  s.  Now  since  *3  (<n)  represents  the  variance 
spectrum  for  V  ’  at  s  =  L,  then  uf 2  v3 (w)  is  the 
variance  spectrum  for  elongational  displacement  of  the 
end  of  tne  line.  Consequently 


IT7! 


s  _ 


i  r 

jjz  \  X*  ^3  (uu)  doi 


(268) 


and 


\I^F  =  (y„  -  iq )  yf772  , 


(269) 
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as  we  should  expect  for  slow  rate  of  stretching.  Here 
(Y0  -  Kx )  corresponds  to  the  quasi-static  spring 
coefficient,  whereas  Y0  is  the  dynamic  coefficient  for 
the  given  mean  tension  T. 

On  the  other  hand  for  a  broad-band  spectrum,  reJ  ation 
(266)  indicates  the  possibility  of  resonant  conditions 
at  frequencies  corresponding  to  zero  value  of  the  sin 
term.  These  frequencies  are  uu  =  +  a)n  where 

=  n  7 r  Ga  /L ,  n=l,2,  ...  ( 267 ) 

For  C3  ~  7000  ft/sec  and  L  ~  14,000  ft  this  gives 
<un  ~  n  7t/2  rad/sec.  Near  a  resonant  frequency 

sin2  ( uiL/Cs  )  =  [(uj  -  u,n)  L/C3]2 

(268) 

cos2  (u)L/Ca )  =  1  . 

The  contribution  to  the  variance  of  T*  at  s  =  L 
due  to  a  resonant  response  centered  at  +  u)n  is  given 
approximately  by 
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where  it  is  understood  that  (-m)  =  (w) .  Now 

0,n  tx  <-  2  to  if  n  ~  4  sec  as  found  in  Chapter  III 
for  nylon  line.  Hence  near  the  resonance 

V  ~  (iq/2  Y0  c2  Tl) 

and  (2?0) 

C2  ~  (Y0/m)1/2  . 

Since  the  resonant  response  has  a  very  narrow  band  width 
(of  the  order  of  \/Y0  tx  rad/sec)  we  can  take  at 

the  value  uun  in  (269)  as  a  suitable  approximation.  We 
then  obtain  upon  carrying  out  the  integration 

4tt  Y  3  t 

=  - - - -  *3  (Oh)  .  (271) 

n  jq  l2  n 

for  resonant  response  at  s  =  L. 
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6.  Perturbations  of  a  Nylon  Mooring  Line  Relative  to  a 
Curved  Equilibrium  Configuration 

It  was  noted  previously  that  the  longitudinal  modes 
and  the  transverse  modes  in  the  plane  of  the  equilibrium 
configuration  are  coupled  to  some  extent  if  the  equilibrium 
configuration  possesses  curvature.  In  the  following  analysis 
we  retain  the  approximation  of  uniform  T  and  for  simplicity 
we  will  also  suppose  the  curvature  d'e’/ds  has  a  uniform 
value.  The  pertinent  coupled  equations  are  (176a,  c), 

(177a,  c),  (178),  and  (179)  in  which  G  ’  is  given  by 
(197).  In  the  latter  relation  we  will  consider  only  that 
case  where  v  is  uniform.  The  above  relations  then  possess 
solutions  of  the  form 

<V,  V*  0‘,  e*,  ?*)  = 

provided  that  the  coefficients 
algebraic  relations : 


(iuum’  +  f 

V  )A:  : 

=  ikT  Ag  +  (de/ds)  A%  . 

(273a) 

ia>m  Ag 

=  iK 

A4  -  T  (de/ds)  A3  , 

(273b) 

itu  ( 1  + 

e ■  ^  A- 

~  /  0 

=  ik  A,  +  (aTT/ds)  A->  , 

(273c) 

iu)  As  - 

=  ik 

4  -  ( d’s/ds )  Aj,  , 

(273d) 

(1  +  i«J7j 

)  a6 

-  As  > 

(273e) 

11 

-  K,  A6  . 

(27  If) 

(  A*  t  1  ^3  »  A4  >  As,  i  Ag  ) 


x  e 


i(ks  +u)t) 


(272) 


An  satisfy  the  following 
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The  last  two  of  these  relations  yield 


As 


1  +  Iuitt  )  A 
[Y0  (1  +  xmi)  -  Kj 


(274) 


The  remaining  four  relations  using  (274)  will  lead  to 

non- trivial  values  of  A  if  k  satisfies  the  characteristic 

n 

relation  obtained  by  requiring  the  determinant  of  the  co¬ 
efficients  in  (273)  to  vanish.  This  leads  to  the  following 
possible  roots  for  k2 : 


k2 

=  j  (k,2  +  kg  2  )  +  (dQ/ds)2 

i{|  -  V)s  +  (V  +  V  +  K*  +  V 

VA>® 

7  1/2 

x(d9/ds)2V  , 

(275) 

V 

■  m  u)2  (i  +  7)/r, 

(276a) 

V 

=  uj(  1  -f  7)(u)  m*  -  i  fv)/T  , 

(276b) 

ka2 

m  U72  (1  +  i  .JU  Ti  ) 

“  [Y0  (1  +  i  w  Tj  )  -  Kj 

(276c) 
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As  a  special  case  we  note  that  for  dg/ds  =  0  the 
two  roots  for  k2  are  simply  kx  2  and  kg2  which 
correspond  to  relation:  (203)  and  (236)  in  the  previous 
analyses  for  transverse  and  longitudinal  wave  modes 
respect' ;e.y.  For  typical  conditions: 


IVAa|  ~  1  * 

I k;!2 Ai2  |  ~  10"2  ,  (277) 

j  (dl/ds)/!^2  |  ~  10“ 2  . 


Accordingly,  in  the  case  of  relatively  small  curvature 
(275)  reduces  approximately  to 


(278) 


for  transverse  and  longitudinal  modes  respectively.  One 
effect  of  the  curvature  is  to  increase  the  rate  of 
attenuation  of  the  longitudinal  waves.  For  the  case  of 
fv/mu)  »  1 


K  =  ±  U/C2  -  i  v’) 


(279a) 


where  v'  is  the  attenuation  coefficient 


(2?9b) 


v  being  given  by  (258)  and  C2  by  (259).  The  longitudinal 

» 

waves  can  lose  energy  not  only  by  hysteresis,  but  also  by 
net  transfer  to  the  transverse  waves  by  coupling  due  to 
the  curvature.  The  transferred  energy  is  ultimately 
dissipated  by  form  drag. 

Returning  to  relations  (273)  we  find  the  following 
approximate  relations  among  Ax ,  Ag ,  A3 ,  A4 ,  making  use 
of  (277): 


For  transverse  modes: 


As  1  =  -  (i/kjJJdl/ds)  Ai 

A3'  =  (kx/a,  (1  +  D)  A,  *,  (280) 

A*  *  =  -  (i  m  <*Ai2)(l  +  k12/k02)(d'e/ds)  Ax  «  ; 


For  longitudinal  modes: 
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A*"  =  -(i/ka)(koVk12)(d0/ds)  4"  , 

A3  =  --(i  muj/kj  2 )  ( 1  +  k12/k02)(d'e/ds)  As"  ,  (28l) 

A4"  =  (muj/ks)  Ag51  . 


Note  that  if  dj/ds  =  0  then  P^',  A4  and  Ax  ”,  A3” 

would  vanish. 

The  general  solutions  for  V  1  and  V  which 

v  T 

are  accurate  to  first  order  in  d@/ds  and  yield  essentially 
zero  values  at  s  =  0,  are 

00 

Vv«(s,t)  =  J  eiuit  +  iki(s  "  L)  dFx  (iu) 

.00 

*  i  [*'iM  - cos 


(282a) 


and 


r^’(s,t)  =  J  sin  kg s  eitut  dF3(u>) 


j  1  H  eili,t  +  iki  (s  -  L)  aFi  («>)  ,  (282b) 
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where  dFx (w)  and  dF3 (®)  remain  to  be  determined  in 

terms  of  the  repi esentaticns  of  V  *  and  V  *  at 

v  T 

s  =  L.  In  the  above  relations 


=  h  -  iii, 

K,,  =  a»/Ca  -  iv*  , 

where  jjL  »  1. 

Let 

CO 

vv'(L,t)  =  J  eiBt  dH,(«) 

—«o 

00 

VT,(L,t)  =  J  e1®*  dH3  (,) 

—oo 

where 

dHj  (a)  dHx,(aj)  =  6  1  ( uo  -  («u)  > 

dH3*U)  cffisC®')  =  6 1  (a)  -  in' )*»(«>)  . 


(283) 


(234a) 


(284b) 


(28?) 


In  (285)  (uj)  ana  $3  ( uu)  are  the  variance  spectra  for 
V  ’  and  V  1  respectively  at  s  =  L.  Moreover  we 

V  T 

will  suppose  that  V  *  and  VT 1  at  s  =  L  possess 
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a  cross-spectrum  ty13(u))  such  that 


dHi  ( ou)  dH3(u}’)  =  ft'  (cu  -  id’)  $13(u>)  duu. 


(286) 


If  we  now  equate  (282a,  b)  at  s  =  L  with  v284  a,  b) 
respectively  we  can  solve  for  dFj  (uj)  and  dF3  (m)  in 
terms  of  dl^  (m)  and  dH3(u)).  Inserting  the  resulting 
relations  in  (282a,  b)  and  neglecting  terms  of  order 
(d"e/ds)2  yields 

CO 

Vv‘(s,t)  =  J  elkl^S  “  elu,t  dH^cu) 

—00 

*!  j^{  (-'ik*S  -  cos  *.) 


_  eiki(s  ~  L)(p”ikiL  _ 


_iujt 

^sfn^L  (287> 


VT'(s,t) 


r  sin  koS  _iu>t  JTT  ,  v 
=  J  sin  gl  e  ■“(■(•) 


00 

I  i{ 


eiK  (S  -  L)  .  ,i*t  ^  Ll  (a88) 


*nl 
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The  solution  for  T*  to  the  same  order  «_f 
approximation  can  be  shown  to  be  given  by 


T 1 (s, t) 


00 


mm 


—00 


cos  ka  s  i^t 
s' in  gL 


«iHa  (a») 


00 

,  dg_  P  nijj,  r  k^_  cos  kas 

+  ds  J  kj 8  Lkg  sin  kgL 
—00 


-  i(l  +  kj2/!^2 )  elkl^s  ■  L^"|  e1U}L  dl^  (»)  . 


(289) 

The  primary  effects  of  the  curvature  are:  (a)  to 
cause  the  transverse  oscillations  to  reach  a  greater  depth 
and  (b)  to  produce  a  somewhat  greater  variance  of  T' 
at  the  upper  end  of  the  line  compared  to  that  v/hich 
would  exist  without  curvature  of  the  line.  At  inter¬ 
mediate  depths  along  the  line,  the  contribution  to  the 
variance  of  V^1  by  the  first  integral  of  (287)  becomes 
negligible  (as  seen  in  the  previous  section)  and  the  second 
integral  gives  a  residual  variance  given  approximately  by 
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(V  «)2 
'  V  ' 


(H)2  ]  \&*  la  «.(.) 

—CO 


(290) 


at  mid- depths . 

It  may  be  remarked  in  closing  that  the  effect  of 
curvature  is  of  more  significance  in  the  case  of  steel 
cables  where  the  net  weight  in  water  produces  a  relatively 
greater  value  of  curvature  than  for  nylon  line,  other 
conditions  being  comparable.  In  fact  for  nylon -line,  the 
simpler  theory  given  in  the  pieceding  section  in  which 
curvature  is  ignored  is  probably  adequate  for  most  conditions. 

Further  application  of  the  methods  discussed  herein 
would  be  profitable,  particularly  in  respect  to  the  use 
of  the  generalized  upper  end  conditions  discussed  in 
Chapter  II.  Unfortunately  this  is  well  beyond  the  scope 
of  the  present  study,  but  it  is  hoped  that  the  present 
study  will  serve  as  a  useful  guide  in  future  research  on  the 
problem  of  mooring  line  dynamics. 
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APPENDIX  A.  EQUATIONS  OF  MOTION  OF  VESSEL 

We  wish  to  establish  relations  connecting  the 
motion  of  a  moored  vessel  with  the  potential  motion  of 
the  sane  vessel  in  the  same  seaway  if  it  were  not  moored 
and  not  underway. 

First,  consider  the  equations  of  motion  for  the 

unmoored  vessel.  Let  X^  (j  =  1,  2,  3)  denote  the  three 

Cartesian  components  of  the  positi.on  of  the  center  of 

the  mass  of  the  vessel  (relative  to  a  fixed  equilibrium 

position).  Let  a  .  ( j  =  1,  2,  3)  represent  the  com- 

0 

ponents  of  angular  displacement  about  the  mutually 
orthogonal,  principal  axes*  of  the  vessel.  Assuming 
the  vessel  is  symmetrical  about  its  mid-section,  one 
principal  axis  is  normal  to  the  plane  of  symmetry  and 
the  other  two  lie  in  the  plane  of  symmetry.  Under- 
equilibrium  conditions  it  will  be  assumed  that  the 
plane  of  symmetry  is  vertical  and  that  the  longi* udinal 
principal  axis  is  horizontal.  It  will  be  assumed  more¬ 
over,  that  the  motions  are  sufficiently  small  such  that 
the  non-linear  coupling  effects  of  the  translational 
and  rotational  components  of  motion  can  be  ignored. 


*  The  principal  axes  are  those  for  which  she  moment 
of  inertia  matr  x  can  be  represented  in  diagonalized  form. 


FoJ  ' owing  Si.  Denis  and  Pierson  (1955)>  the  linearized 
equations  of  motion  for  an  unmoored  vessel  are  taken 
as  follows  (j  =  1,  2,  ): 

M.  d2\./at2  +  H.  o’X  ,/dt  +  K_.  X,  =  F , '  ,  (A-l) 

J  J  j  J  <J  U  <J 

I.  a3  a,  ./at2  K.*  da  ,/dt  +  K  *  a.-  =  H,f  ,  (A-2) 

JO  J  u  J  J  J 

where  the  M  are  the  effective  masses  for  surge,  sway, 

J 

and  heave  (including  the  effective  added  mass  of  water), 

the  N  .  are  damping  coefficients  associated  with 
J  * 

viscous  drag  and/or  radia'ion  of  ship  induced  waves, 

the  K.  are  restoring  force  coefficients  (force  per 
J 

unit  displacement)  and  the  F  *  are  the  components  of 
exciting  force  due  to  wave  motion  and/cr  steady  motion 
of  the  water  relative  to  the  vessel.  Similarly  the 
I  are  the  effective  moments  of  inertia  about  the 

.j 

orinciual  axes  of  the  vessel,  tne  N '  are  damping 

u 

coefficients,  the  K  !  are  res’ oring  coefficients 

(torque  per  radian)  and  the  H.'  are  components  of 

0 

exciting  torque  due  to  waves.  It  will  be  ’understood 
specifically  that  Xl  ,  X2 ,  X3  are  respectively  the 
translational  displacements  in  surge,  sway,  and  heave. 
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moreover,  ax ,  ,  a3  are  respectively  the  rotational 

displacements  in  roll,  pitch,  and  yaw.  Accordingly 
Ki ,  ,  and  *  are  zero;  the  only  restoring  forces 

being  those  for  heave,  roll,  and  pitch. 

We  can  rewrite  (A-l)  and  (A-2)  in  the  form 

d2X  ./dt2  +  23  4  dX./dt  +  02  .  X,  =  f  ,  (A-3) 

J  J  J  j  0  j 


u2a  ./dt2  +23.'  da./dt  +  a  »2  X  .  «  h . »  ,  (A-4) 

J  J  J  j  .j  u 

wnere 


sj 

5  N  ,/2M.  , 

J  « 

j 

5  Y/2Ij  • 

o . 

J 

s  Vk./m.  , 

J  J 

o . ' 
J 

/ -  (A-3) 

:i  Vk-'A  .  , 

J 

f 

J 

< 

•o 

ill 

and 

o 

h.* 

0 

E  • 

The  four 

sets  of  coefficients  p  .,  3.-’,  a  and  a.’ 

J  j  j  J 

all  have 

units  of 

frequency:  the  3  terms  govern  the 

damping  rate  for  free  motion. 


The  natural  frequencies 


i  r>  ^ 

xuv 


( rau^ans/seeond) ,  for  tne  case  where  p3  «  o?  ,  pj  1  «  cx  ' 
and  p2  1  «  as  are  given  approximately  by  a3 ,  ax  '  and 
o2  '  (for  heave,  roll,  and  pitch  respectively).  Again  note 
that  ox  ,  a2  and  o3  '  are  zero.  It  will  be  presumed 
that  all  six  damping  rates  and  the  three  natural  frequencies 
for  the  unmoored  vessel  are  known. 

For  the  moored  vessel  we  must  include  the  effect  of 
the  tensile  force  for  each  mooring  line.  Let  T,  (t) 

K 

denote  the  tensile  force  exerted  by  mooring  line  k 

(k  =  1,  2,  .  .  .)  at  its  point  of  attachment  to  the 

vessel.  Let  p  represent  the  vector  separation  between 
K. 

the  vessel's  center  of  mass  the  the  kth  attachment 
point  (see  Fig.  2  ).  Recalling  the  convention  for  T, 
we  note  that  the  force  applied  to  the  vessel  by  line 
k  is  -T  .  Likewise  line  k  exerts  a  torque  given 

xC 

^  A 

by  -p  x  Tk.  Let  T-  -  represent  the  jth  Cartesian 
component  of  the  vector  Ti  and  let  pkj  be  the  jth 
component  of  the  separation  vector  p,  ,  the  components 
being  referred  to  the  same  coordinate  system  as  that 
for  and  a,.  Then  pk  x  Tk  has  the  components 

Jkj-  given  by: 


l8l 


\  =  (Pk,  Tk,  '  P^  Tk» ) 

\  *  K  \  -  Pk,  Tka ) 

Jkg  -  (Pk,  Tks  -  Pka  Tk,  )  • 


If  we  now  let  R.  denote  the  translational  dis- 
J 

placements  of  the  center  of  mass  for  the  moored  vessel 
and  ijr  .  the  angular  displacements  then  the  equations  of 

t) 

motion  for  the  moored  vessel  take  the  form 

d2RJ./dt2  +  20 j  dR  /dt  +  of  R^  =  f  j  -  ^  (A-7) 

d2  i  j/dt2  +  23j-«  dijf  j/dt  +  afs  =  hj  JkJ^Lj  (A~8) 

vhere  J,  .  is  given  by  (A-6)  and  the  terms  f .  and 
KJ  J 

h.  are  the  counterparts  of  f,s  and  h.1  for  the 
J  J  J 

moored  vessel. 

Suppose  (A-3),  (A-4),  (A-7),  and  (A-8)  are  applied 

for  the  same  vessel  and  for  the  same  sea  state.  Moreover 

suppose  for  the  unmoored  condition,  that  the  vessel  is 

not  underway.  It  is  reasonable  to  suppose  then  uhat  the 

mean  values  of  the  exciting  terms  f .  and  h  are  zero. 

0  J 
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This  is  not  true  for  the  moored  situation, 
wila  be  assumed  that  the  departures  of  f 
from  their  mean  values  are  the  same  as  f 

0 

i. e. , 


However,  it 

and  of  h . 

3 

and  h  . 1  : 

3 


and 


(A-9) 


n . 

o 


=  h  .  ’ . 
3 


We  expect  that  the  mean  values  of  the  time  derivatives 

of  R  and  of  a  .  will  vanish  and  hence 
3  3 


a.2  R. 
J  3 


\  A 


a . J  +  i  J  A  . 
J  3  K  kj  j 


(A-10) 


Consequently  from  (A-3),  (A-4),  (A-7),  (A-8),  (A-9),  and 
(A-10)  we  obtain  the  following  relations  for  the  -  anomalies 
of  ship  motion  induced  by  the  mooring  line: 


(d2  /dt2  +  2p  d/dt  +  R/ 


=  -  L  T  ’/&  , 

k  "J  J 


(A-ll) 
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(d2/dt2  +  2B  .  *  d/tdt  +  a  . 12  )  ^  .  * 
«J  J  J 


..l 


k  J 


(A-12) 


where 


R.»  =  R 

J  s 


(A-13) 


It  is  understood  that  \  and  a  are  prescribed 

0  0 

functions  of  t;  these  are  the  potential  displacements 
which  vould  exist  in  the  absence  of  the  mooring  lines 
for  the  same  sea  state. 

To  these  relations  we  must  add  a  kinematical 
relation  governing  the  orientation  of  the  vectors  p  . 

For  a  rigid  vessel  the  separation  distances  Jpk| 
are  constant,  however,  the  orientation  is  governed  by 
the  vector  angular  velocity  of  the  vessel,  Q, 
according  to  the  relation 

A  * 

dp  /dt  =  n  x  p  (A-l4) 

K  K 

A 

The  Cartesian  components  of  Q  are  given  by  d$  ./dt. 

<3 

H«nce  (A-14)  can  be  expressed  in  the  following  component 


form: 


(d/dt)  pk2 
(d/dt)  pk2 

and 

(d/dt)  pk? 
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(d'V/dt)  pk3  -  (d\jj3/dt)  pk2 
(d^/dt)  pkl  -  (d^/dt)  pk3 
(d^/dt)  pk2  -  (d^/dt)  pkl 


(A-18) 


The  restriction  has  already  been  imposed  that  the  motions 

of  the  vessel  are  small.  Consequently  we  may  approximate 

relations  (A-15)  by  taking  the  mean  values  of  the  terms 

p  .  on  the  right  hand  sides  and  obtain  by  integration, 
kj 

using  (A-13): 


Pkl‘  =  +  Pk3°  "  (*3*  +  as)  Pk2°  » 

Pk2 '  =  (♦s*  +  aa)  Pkl°  "  (*1  '  +  ai  )  Pk3°  »  (A-16) 

Pk3'  =  (*1  '  +  <h)  Pk2°  -  (V8  ’  +  a8)  Pkl°  > 


wnere 


pw'  -  pkJ  -  pk.i 


and  p.  .  denotes  the  mean  values  of  p  .. 
K,l  kj 


(A-17 ) 
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Finally  to  the  above  relations  we  added  the  following 
geometrical  relation  for  each  line 

rlc  -  R  +  pk  +  °k  <A-l8> 

A 

where  r^  is  the  radius  vector  from  the  anchor  point  of 
line  k  to  the  point  of  attachment  on  the  vessel,  while 

a 

Ck  is  a  constant  vector  representing  the  radius  vector 

from  the  anchor  point  for  line  k  to  the  fixed  reference 

* 

point  for  the  displacement  vector  R.  The  above  relation 
can  also  be  written  in  the  form 

A  A  A  A  A 

r,  =  r,  °  +  R‘  +  p.  ’  +  x 
k  k  k 

A 

A 

where  r, 0  is  the  mean  value  of  r  . 

K  k 

We  note  in  passing  that  the  vector 
written  in  the  form 

Pk  =  P^0  +  U:  +  a)  x  pko  .  (A-20) 

This  is  consistent  with  (A-l6)  and  (A-17). 
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APPENDIX  B.  DERIVATION  OF  THE  CHARACTERISTIC 
FORM  OF  THE  MOORING  LINE  EQUATIONS 


The  equations  governing  the  local  time  rates  of 
change  of  the  dependent  variables  Vv,  Vu,  V  ,  0,  §,  e,  %, 
and  e0  are  given  by  (109),  (HO),  (111),  (93),  (9^), 
(95),  (70),  and  (71),  respectively.  These  equations  also 
contain  the  dependent  variable  T,  however,  we  have  the 
diagnostic  equation  (69)  relating  T  to  e,  §  and  e0 . 
We  wish  to  recast  these  equations  in  the  form 


8 


P, 


(B-l) 


along  the  path 


ds 

at 


(B-2) 


Here  the  Q^  signify  the  eight  prognostic  variables  and 
the  quantities  B  ,  P  and  C  are  functions  of  these 
variables.  In  general  there  will  be  eight  possible 
values  cf  the  speed  C  (not  necessarily  all  distinct) 
and  eight  associated  relations  of  form  (B-l).  All 
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of  the  admissible  values  of  C  will  be  real  if  the 
system  of  differential  equations  is  hyperbolic,  which 
is  true  in  the  present  case. 

First.  of  all,  we  obtain  from  (69) 


dT 


3e  de 


*  <®-^>  -  § . 


(B-3) 


where 

Y  h  K  eb(€  “  e°)  .  (B-4) 

This  relation  can  be  employed  in  (111)  so  as  to  eliminate 
any  dependency  on  derivatives  of  T. 

In  order  to  put  the  system  of  equations  in  the 
form  (B-l)  with  condition  (B-2),  we  form  a  quasi-linear 
combination  of  the  original  set  in  the  following  sense. 
Multiply  (109),  (110),  and  (111)  by  bj,  b2,  and  b3 
respectively;,  multiply  (93),  (9*0,  ^  (95)  by  b4,  bB , 
and  b8  respectively;  multiply  (70)  and  (71)  by  b7  and 
ba,  respectively;  and  add  the  resulting  equations.  The 
combination  is  quasi-linear  in  the  sense  that  the  co¬ 
efficients  bx  through  b8  do  not  involve  derivatives 
of  the  dependent  variables.  After  collecting  terms,  the 
result  is: 


b 
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bv  bv 

«>i  m’  'Jt2  ‘  b«  ’5T"  ]  +  [b°  m 

bv  bv  n 

•f  [  b3  m  -  be  5  +  [bx  m*  VT  -  b3  Xa  +  (1  +  e)b4]|S- 

+  (-  \  T  -  b4  +  b6  V  .)  ||  ] 

+  [  (-  bj  m'  V  cos  0  +  b8  m'  Xi  -  b3  X3  sin  6 

+  b5  (1  +  e)  sin  0)  +  (-  b2  T  sin  8  +  b4  sin  0 

-  bs  Xj  +  b6  sin  0)  |f  ] 

+  lbe  ||  "  b3  Y  ||  ]  +  [b7  ||  +  b3  K,  ||  ] 
de0  8e0 

+  [b6  'gy ~  +  b3  Y  =  bi  (w  sin  0  +  Gv) 

+  bs  Gn  +  b3  (~  w  cos  9  +  GT)  +  '°7  (e  -  c0  -  §)Ax 

+  b8  S  (T,  eQ)  .  (B-5) 


bV 

_1±  - 
dt 


3vm  , 
■SJ"  1 
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In  the  above  relation 


(B-6) 

if  £0  <  F(T) 

(B-7) 

if  e0  2  F(T) 

and 

m’  =  m  +  ma.  (B-8) 

The  coefficients  bj  through  b8  are  to  be 
selected  such  that 


=  -  C  m*  bx  , 

(B-9) 

b5 

=  -  C  m '  ba  , 

(B-10) 

b6 

=  -Cm  b3  , 

(B-li) 

b3 

Y  =  -  C  b6  , 

(B-12) 

*>3 

Ka  =  C  b7  , 

(B-13) 

b3 

Y  =  C  b8 

(B-14) 

t- 
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(*Tb1  -  VT  b4  +  V.  b6)  = 

C  (m1  Vt  bi  -  Xa  b3  4-  (1  +  e)  o4),  (B-15) 


and 


(-  T  sin  0  ba  +  cos  0  b4  -  Xx  bg  +  sin  0  b5 ) 


=  C  (  -  m1  V  cos  S  bj  +  m1  5^  b  -  X3  sin  0  ba 


+  (1  +  e)  sin  0  b5 ) . 


(B-16) 


each  of  the  dependent 
in  the  form 


which  is  the  same  as  dQ^/dt  along  the  line  ds/at  =  C. 

Using  (E-9),  relation  (B-15)  can  be  simplified  to 
the  form 


The  latter  relations  assure  that 
variables,  in  (B-5)  enters 


(■5 r  +  0  sr  > 


(-  T  \  +  Vv  b.) 


C  (-  Xe  b3  +  (I  +  e)  b„).  (B-17) 


-tyj. 

Mcreover,  using  (B-9)  and  (B-10),  relation  (B-l6)  can 
be  reduced  to  the  form 

(-  T  b2  +  b6)  C  (-  Xa  b3  +  (1  +  e)  bs),  (B-18) 

provided  that  sin  0  does  not  vanish. 

Since  the  eight  homogeneous  relations  (B-yy  through 
(B-14)  plus  (B-17)  and  (B-l8)  contain  eight  parameters, 
bn,  plus  an  eigenvalue  C,  it  follows  that  at  least 
one  of  the  bR  for  given  C  is  arbitrary.  The  eigenvalues 
C  can  be  evaluated  as  the  roots  of  the  eighth  degree 
polynominal  in  C  found  from  the  requirement  that  the 
determinant  of  the  coefficients  in  the  simultaneous 
system  of  equations  for  b^  vanishes.  This  determinant 
in  the  present  case  contains  many  zeros  and  can  be  shown  to 
lead  to  the  following  eighth  degree  polynominal: 

C2  (Y  -  Cam)  (C*  -•  (1  +  e)  -  T)»  =  0.  (B-19) 

Therefore,  C2  can  have  the  possible  roots 

0,  Y/m  or  T/m1  (1  +  e)  , 


the  last  having  a  multiplicity  of  two. 


For  the  case  C2  =  0  we  find  that  the  relations  for 
are  satisfied  if  bx  =  =  b3  =  b4  =  b5  =  bs  =0 

for  arbitrary  'o7  and  be .  Thus  for  one  of  the  two  roots 
C  0,  take  b7  =  1,  b8  =  0  and  for  the  other  take 
b7  =  0,  ba  =  1.  This  implies  that  we  merely  recover 
relations  (70)  and  (71)  in  their  original  form,  as  might 
have  been  anticipated,  since  they  involve  only  derivatives 
of  §  and  eQ  respectively. 

For  the  case  C2  T/r. ’  (1  +  e),  one  possible  set  of 
bfi  are  found  to  be  b^  =  =  bs  =  b6  =  b7  =  bs  =  0 

with  bj  arbitrary  and  b4  given  by 

b4  =  -C  m1  bj  . 

Thus  if  we  take  bx  =  1,  then 

b4  =  +  l  m«  T/(l  +  e)]1/2 

for  (B-20) 

C  =  +  [T/m*  (1  +  e)]1/2. 

A  second  possible  set  of  bfi  for  the  case 
C2  =  T/m1  (1  +  e)  is  found  to  be  bj  =  b3  =  b4  =  bs  = 
b7  =  b8  =0  with  ba  arbitrary  and  bs  given  by 
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bs  =  -  C  m  *  ba . 

Thus  with  bg  =  1  we  get 

bs  =  7  [ni*  T/(l  +  e)]x/2 

for  (B-21) 

C  =  +  (T/m*  (1  +  e)]1/2. 

Fina.lly  for  the  case  Ca  =  Y/m  we  get  from  (B-9) 
to  (B-i4)  with  (B-17)  and  (B-l8)  the  following  relations 
for  the  bn  : 

*.  -  -c  »a  (Vv  -  V'2  ’ 

ba  =  -0  ma  (Vu  -  Uu)/2  , 

bs  =  1  , 

fc»  =  Y  m'  ma  (Tv  -  Uv)/mZ  ,  (B-22) 

bs  =  Y  m<  ma  (V^  -  D^/mZ  , 
b6  =  -Y./C, 

b,  Kx/C  , 
be  =  Y/C 


i 
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where 

C  =  +  (Y/m)1/2 

and 

N. 

Z  h  Si  (1  +  e)  Y  -  T.  (B-23) 

Since  the  spring  coefficient,  Y,  is  generally  much 

greater  than  T,  the  quantity  Z  is  positive  (even 

when  T  is  at  or  near  its  ultimate  value). 

The  relations  of  form  (B-l)  follow  immediately 

from  the  above  results  by  substituion  of  the  b  into 

n 

(B-5)  for  each  particular  C.  In  the  case  C  =  +  (Y/m)1^2, 
the  terms  involving  de/dt,  d§/dt,  and  de„/dt  can  be 
recollected  using  the  relation 


The  resulting  characteristic  form  of  the  mooring  line 
equations,  making  use  of  definitions  (B-6),  are  those 
presented  in  Section  4  of  Chapter  V. 


(B-24) 
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\ 

APPENDIX  C.  TRANSFORMATION  OF  THE  CHARACTERISTIC 
EQUATIONS  TO  CARTESIAN  FORM 

Using  (88)  gives  the  following  relations  between 
the  natural  and  Cartesian  velocity  components 

V  —  -V  sin  $  +  V  cos  $  , 

U  x  y 

V  =  V  cos  e  cos  *  +  V  cos  e  sin  *  -  V  sin  g  , 

v  x  y  z 

(C-l) 

V  =  V  sin  0  cos  4  +  V  sin  0  sin  4  +  V  cos  $  . 

t  x  y  z 

These  are  employed  in  (115),  (116),  and  (117)  to  replace 

dV  /dt,  dV  /dt,  and  dV  /dt  by  the  time  derivatives 
U  v  t 

of  Vx'  V  V  311(1  #* 

Now  from  the  definitions  given  in  (132)  we  find 

< 

&  =  cos  4  cos  e  ||  -  sin  4  sin  e  ||  , 

(C-2) 

-I-  =  sin  i  cos  e  +  cos  $  sin  0  ■§+*  • 
dt  dt  dt 


These  relations  are  employed  to  eliminate  dg/dt  and 
d$/dt  by  appropriate  combination  of  the  transformed 
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versions  of  (115)..  (116),  and  (117).  Finally  using  (132) 

*  the  coefficients  involving  sin  0,  cos  0,  sin  $  ,  cos  $ 

can  be  put  in  terms  of  <x>  0,  y  The  resulting  relations 

are  those  given  by  (137)  to  (1?9)-  The  last  of  these 

was  given  in  abbreviated  form  involving  a  residual  term 
4- 

6  .  The  complete  form  1'or  (139)  is  as  follows: 


dV  dV  aV  -j  /0 


,1/2 


(Gt  -  wY)/m  ^  (Y/m)  [Kj  (e  -  e0  -§)ATi 


+  S(T,  O]  +  (Y/m) 1/2  j|  l  (ov  +  »/a *-T-s»)  (Vv  -  uv) 


,  9  maT 

+  <VM  -  V1  -  W-  '■ 


<Vx  -  Ux>  It  +  <Vy  -  V  f 


+  <Vz  -  V  &  ]  +  (V®)172  ^  (Vn  -  0n)  •  •!}  (0-3) 


along  ds/dt  =  +  (Y/m)  '  . 

The  lengthy  set  of  terms  in  braces  represents  the 
4- 

residual  6  .  In  this  relation 


«  f—  (i  +  c)  y  -  t. 


(c-H) 
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Fortunately  under  normal  conditions  the  residual 

6—  is  of  very  minor  consequence  and  can  be  neglected 

for  practical  purposes.  To  justify  this  we  note,  first 

of  all,  that  Z  is  of  the  order  of  magnitude  of  Y. 

Thus  the  terms  involving  G  and  G  are  of  the  order 

v  u 

of  magnitude 


(V  -  U  ) 
v  n  ny 

— r — 


9 


1/2 

where  C”  =  (Y/m)  '  .  Under  normal  conditions  we 

would  expect  that  (Vn  -  Un)  is  of  the  order  of  magnitude 
of  10  ft/sec  or  less.  However,  C”  is  of  the  order  of 
10,000  ft/sec.  Thus  the  above  terms  are  negligible  com¬ 
pared  with  (G^  -  wy)/m. 

The  terms  involving  dc/dt,  dp/dt,  and  dY/dt 
are  of  the  ordex' 


(c  VC")*  (Vn  -  0n)  $  , 

where  C'  is  the  signal  speed  for  transverse  waves. 
The  ratio  (C'/C”)2  is  of  the  order  1/100.  If 
(V  -  U  )  dS/dt  is  of  the  order  /dt,  then  the 

above  terms  can  be  neglected. 
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Finally,  the  terms  involving  dV/dt  in  the  residual 
6—  are  of  the  order 

(V.  -  U  )  dV 
v  n  n'  _ n 

C"  at 

The  ratio  (V  -  U  )/C"  is  of  the  order  10-3,  hence 
'  n  n 

these  terms  are  negligible  compared  with  a  8  V  /dt, 

A 

0  d  V^/dt,  and  Y  d  Vz/dt. 

It  will  be  noted  furthermore  that  all  of  the 
terms  in  the  residual  are  proportional  to  m  . 

These  terms  would  therefore  not  exist  in  the  absence 
of  the  fluid.  The  primary  influence  of  the  fluid  is 
clearly  on  the  transverse  waves  and  not  on  the 
longitudinal  waves. 
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APPENDIX  D.  MATRIX  OF  COEFFICIENTS  FOR 
THE  NUMERICAL  PREDICTION  RELATIONS 

The  coefficients  Mn  m  and  Pn  of  relations  (154) 
are  given  in  the  fo], lowing  pages.  These  are  arranged  in 
order  according  to  n. 
n  =  1 

Mn  =  {2  -  [  a2(o,k+l)  +  a3  (Ax  > 

*iz  =  -[aUk+l)  e(j,k+l)  +  a(A1)  P(AX  jj/B,  , 

^i3  =  +  aCAj  )  y(A1  )]a 

M ! 4  =  'I, 

*xs  -  *ie  =  M,7  =  Mle  =  0  , 

8 

P>  -  (2  !«.,«  VA>)  +  [Gnx<J-k+1)  +  GnxK )  . 

m=l 


+  w(a(j,k+l)  Y(j,k+1)  +  afA*)  Y(A1))J  , 

Bi  =  {[(1  +  e(j,k+l))  T(j,k+l)/m’]i/ 2 
+  [(1  +  e(A,  ))  T(A1)/m«]l/a  } 
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M*,  =  {2  -  [a2  (J,k+1)  +  a*  (A*)]}  A,  » 

Mss  =  -  [a(.i.k+l)  +a(Aj)  , 

Mb3  =  “  j_a( j>k+l)  y(«j*k+l)  +  a(Ag)  yC^) > 

^4  -  1  » 

=  H,.  =  =  M2e  =  0  , 

8 

P*  “  t2  I  *V>,»  Sn  (*.)  +  tGnx  («J>fc+1)  +  Gnx  (A.) 
m=l 

+  w(a(j,k+l)  y(j,k+l)  +  a(Ag)  y(4  )]  6t2/m‘j-/B2  , 


Ba  =  {[(1  +  e (tJjk+1) )  T(o, k-Uj/m’]1/  2 


+  C(  1  +  e(A.)  T^/m']1'3}  , 


CJ 


1 

2 


I  r  r  T(j  k+D 

1  Lm*(l  +  eU,k+l)) 


+ 


n  =  3 

^3i  ~  »  (&3  ~  A1  )  * 

ttja  *  {2  “  j_38(J*k+-0  +  ^(A,)]  ]A  , 

**33  =  ^  L£(«b*+1)  y(j>k+l)  +  g{Aj)  y(Ai)JA  j 

^4=0, 

^5  *  “I  3 

^33  =  i-J3  /  ~  8  —  0  , 

8 

P3  =  {  2  j  ^  (A  )  +  [Gny(j,k+1'  +  Gny(Al) 

m=X 

+  w  (&(3f**l)  y(j,fcfl)  +  p(Ai)  y(A  )J  Stj/m^A  , 

and 

rj »  —  ft 

°3  ~  W  • 


$ 


3 

I 

I 


I  • 

f 

! 


I  1  '],t;)i  iMWtaW*  «*w »v*v* 


202 


M41  =  M22  (A4  =  A.) 

m42  •=  {2  -  [e2  (3»k+1)  +  e2  (MTj-A 

M43  =  •  v(J»k+1)  +  P(Aa)  Y^jj/^a 

m44  =  0 

M4S  =  1 

M46  =  «47  =  m48  =  0 

8 

P4  =  {a  Z  M4m  ^  (A,)  +  [g  (5,w)  +  Gny  (A,) 
m=I 

+  w  (p(j,k+X)  y(«J,k+l)  +  p(A,)  yfA,))]  6ta/m*}/Ba 


C4' 


=  <y 


and 
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n  =  5 

Mgj  =  [<x(j,k+l)  +  <x(A5)]/B3 

Msa  =  [p(J,k+l)  +  p(As)]/B3 

M53  =  Lv(0>k+1)  +  y(a5)]/b3 

M54  «  =  0 

Mse  =  "I 

Ms7  =  =  0 

8 

Ps  =  {2  l  “sm  \  +  [GtUK+1)  +  Gt(*s) 

m=l 

-  »  (Y(j,k+1)  +  y^))]  fitE/m 
+  C"  [  P(j,k+1)  +  P'{AS )]  fits}  /% 

B3  =  1/m  C”  ,  %  (c  -  e„  -  ?)  +  S{T,  Eo ) 

C"  =  \  [  (Y(j,k+l)/m)i/  3  +  (Y^J/m)*'  3  j 
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n  =  7 


M  71  =  M,2  =  Hyg  =  M,4  =  Mys  =  M76  =  0 
M77  -  1,  M,a  =  o 

P7  =  {  [(l/b)4nLY(j,k+l)  Y(j,k)l  -  ?(J,k)]  At/2Tx 


+  /  (1  +  At/2Tl) 


n  =  8 

M81  =  M82  =  M0  3  —  M84  =  Mg  5  =  M86  =  0 

^7  =  0,  M88  =  1 

Pe  =  Uo  +  F(T(o>k))  At/To}/(l+  At/To)  , 

If  e0(j,k)  <  F(T(j,k)) 

Fe  =  e0(d,k).  If  e0(o,k)  >F(T(j,k))  . 


a 


LIST  OP  SYMBOLS 


Cross-sectional  area  of  mooring  line  at 
s,  t;  also  used  as  an  arbitrary  constant 

End  point  of  characteristic  line  (Pig.  l6) 

Value  of  s/as  for  point  A  ; 

Also  used  to  denote  norm?l  acceleration 

Archemedian  force  on  line 

End  point  of  characteristic  line  (Fig.  17) 

Coefficients  in  characteristic  equations 

Race  of  change  of  the  dynamic  spring 
coefficient  with  tension,  see  Eq  (55) 

Interpolated  relative  time  for  point  Bn 

Signal  speed  (C  *  for  transverse  modes, 

C"  for  longitudinal  modes) 

Transverse  and  longitudinal  signal  speeds 
corresponding  to  conditions  of  the  mean 
state  of  the  line 

See  Eq  (144) 

Added  mass  coefficient 

Drag  coefficient  (form  drag) 

Effective  diameter  of  the  line  at  s,  t 

Nominal  diameter  of  the  line  in  its  original 
unstrained  state 

Mean  rate  of  energy  dissipation  by  form  drag 
per  unit  material  length  of  line 
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dFnU) 

dHnU) 

dQn(ou) 


Stochastic,  complex  differentials  employed 
in  Fourier,  Stieltjes  representations  of 
stochastic  variables 


Total  mechanical  energy  per  unit  material 
length  of  line 

Elastic  energy  per  unit  material  length  of 
line 

Effective  kinetic  energy  per  At  material 
length  of  line 

Potential  energy  per  unit  material  length 
of  line 


F 


Total  fluid  force  on  line 


F(T) 


'V  fd»  fl’  f£ 


f(a) 

f 

A 

G 


Gv,  Gt 


G  ,G 
nx’  ny 

G’,G’ 
v'  U 


G(X) 


Saturation  function  giving  e  vs  T  for 
slow  unidirectional  straining 

Fluid  force  per  unit  material  length  of  line 

A 

Contribution  to  f  due  respectively  to 
static  buoyancy,  form  drag,  acceleration  of 
line  relative  to  fluid  and  lateral  thrust 
due  to  vortex  shedding 


A  scalar  function  of  a 

a 

Coefficient  in  relation  for  G  (Chap.  IV). 
see  Eq  (173) 

A 

That  part  of  residual  fluid  force  R  not 
dependent  upon  d*/dt,  see  Eq  (105)or  (106) 

Natural  components  of  G 

a 

x  and  y  components  of  the  normal  part  of  G 

Anomalies  of  G^  and  G  ,  see  Eq  (175) 

Kernal  function  in  an  integral  representation 
for  T,  see  Eq  (73) 


h(tu) 


Acceleration  due  to  gravity 

Hysteresis  for  a  closed  cycle  (energy  loss 
per  cycle  per  unit  material  length  of  line) 

Hysteresis  for  an  idealized  process 

Maximum  value  of  H 

Transfer  function  related  to  G(\) 

Moments  of  inertia  of  vessel  about  its 
principal  axes  for  roll,  pitch,  and  yaw 
(j  =  1,2,3  respectively) 

Unit  vector  in  the  x-direction  a  fixed 
Cartesian  frame  of  reference 

/-l;  also  used  to  denote  an  integer 

The  negative  of  the  torque  exerted  on  the 
surface  vessel  by  mooring  line  k 

A  A 

Anomaly  of  from  its  mean  value  ( J^0  ) 

Cartesian  components  of  (j  =  1,2,3) 

The  value  0  if  <  0,  N  if  an>  N 

The  integer  truncation  of  the  number 

Unit  vector  in  the  y-direction  for  a  fixed 
Cartesian  frame  of  reference 

An  integer 

Dynamic  spring  coefficient  at  zero  tension, 
see  Eq  (55) 

Constant';  with  units  of  load  per  unit  strain 
characterizing  the  elastic  properties  of  the 
line 

Apparent  dynamic  spring  coefficient  for 
cyclic  loading  at  frequency  m 


t 


Unit  vector  in  the  2-direction  (upwards) 
for  a  fixed  Cartesian  frame  of  reference 

An  integer;  also  denotes  a  complex  wave 
number  (Chapter  IV). 

Constants  in  the  representation  of  the 
damping  characteristics  for  transverse 
oscillations  of  the  line 

Original,  unstrained  length  of  li; 

(s  =  L  at  point  of  attachment  to  vr  t-?1) 

Total  mass  of  line  (directly  a  function  of 

s) 

The  effective  inertial  mass  of  the  surface 
vessel  for  surge,  sway,  and  heave 
(j  =  1,2,3  respectively) 

Matrix  of  coefficients  in  the  equations 
for  see  Eq  (154)  and  Appendix  D. 

Mass  of  line  per  unit  material  coordinate 

Effective  added  mass  of  water  per  unit 
material  coordinate 

Mass  of  water  displaced  by  line  per  unit 
material  coordinate 

m  +  ma 

Frictional  coefficient  (  =  Kt  )  for  line; 

also  an  integer  representing  the  maximum 
number  of  intervals  of  s  in  the  numerical 
grid 

Notation  for  anchor  point  of  mooring  line 

Denotes  a  material  point  on  the  line;  also 
denotes  period  (=2rr/w). 

Forcing  terms  in  the  characteristic  equations 
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A 

Pk 

Pkj 


Forcing  term  in  Eq  (154),  see  Appendix  D 

Separation  vector  between  center  of  mass 
of  vessel  and  point  of  attachment  of  mooring 
line  k  (See  Fig.  2) 


Mean  value  of  vector  p^ 

A 

Cartesian  components  of  p  ( j  =  1,2,3) 

Denotes  any  dependent  variable  of  the 
mooring  line  at  s  v.  t  ■  but 


Q, 


'n 


A 

V 

A 

R 


R* 


The  set  of  dependent  variables  defined  by 
Eq  (152),  n  =  1,  2  .  .  .8 

Cartesian  components  of  the  velocity  of  the 
center  of  mass  of  the  surface  vessel 

Residual  fluid  force  per  unit  material 
coordinate,  see  Eq  (lo) 

Anomaly  of  the  translational  displacement 
of  the  surface  vessel 

A 

Components  of  R*  in  surge,  sway  and  heave 
(j  =  1,2,3  respectively) 


r 


S(T,eo) 

s 


T 


Position  vector  a  material  point  on  the 
line  with  refererce  to  the  anchor  point 

Position  vector  at  point  of  attachment  of 
line  k  on  the  surface  vessel 

A 

Mean  value  of  r. 

k 

A  function  defined  by  Eq  (119c) 

Material  coordinate,  representing  the  arc 
length  between  0  and  P  in  the  original 
relaxed  state  of  the  line  (a  function  of  M 
alone) 

A  A 

Vector  tension  in  line  at  s,  t  (T  =  Tt) 
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Value  of  T  at  anchor  point 

a 

Value  of  T  at  point  of  attachment  of  line 
k.  to  surface  vessel 

a 

Mean  value  of 

A 

Cartesian  components  of  ( j  -  1,2,3) 

A  A 

Anomaly  of  from  its  mean  value  (T^0 ) 

Scalar  tension 

Dynamic  tension  for  given  e  arid  e0 

Saturation  value  of  T  for  given  e0 

Uluzmate  strength  of  line 

Mean  value  of  T  at  s 

Perturbation  of  T  from  its  mean  value 
at  s 

time 

Velocity  of  fluid  (in  absence  of  line)  at  s,t 

A 

That  part  of  U  normal  to  the  line  at  s,t 

A 

Natural  components  of  U 

A 

Cartesian  components  of  U 

A 

Cartesian  components  of  dU/dt 

A 

Velocity  of  line  at  s,t  (dr /St) 

A 

v  -  -  V 

'  T 

A 

Natural  components  of  V  at  s,t 

A 

Cartesian  components  of  V  at  s,t 

A 

Components  of  V  in  a  natural  coordinate 
system  defined  by  the  mean  configuration  of 
the  line  ™ 


Effective  mean  value  of  the  relative  normal 
velocity,  defined  by  Eq  (193) 

Net  weight  of  line  in  water  per  unit  material 
coordinate 

HorizontalCartesian  coordinate  in  the  plane  of 
the  equilibrium  configuration  of  the  line  (for 
conlanar  conditions) 

The  general  spring  coefficient  dT/de,  see 
Eq  (117c)  or  (121) 

Value  of  Y  corresponding  to  T  and  If 

Horizontal  Cartesian  coordinate  normal  to  the 
Coplanar  equilibrium  configuration  of  the  line 

.  See  Eq  (126) 

Vertical  Cartesian  coordinate 


A  vector  denoting  the  potential  rotational 
displacement  of  the  surface  vessel 

Cartesian  components  of  a  in  roll,  pitch, 
and  yaw  (j  =  1,2,3  respectively) 

A  dummy  variable  to  denote  e  e0  ;  also  used 
to  denote  the  x-component  of  t 

Damping  factors  pertinent  to  the  six  components 
of  motion  of  the  surface  vessel  (p^  for  translation, 

r.  ’  for  rotation,  j  =  1,2,3)  J 
J 

A 

The  y-conponent  of  t 

The  anomaly  of  p  from  its  mean  value 

See  Eq  (107);  also  used  in  a  different  sense  in 
Table  3 

a 

The  z -component  of  t 
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As 

AT 

At 

6 

6  '(a) 

6± 

6t 


e0 


ex 


e  ' 

C 

U 

0 

0* 

0> 

X 


Increment  of  s 

Amplitude  of  T  under  cyclic  loading 
Uniform  time  step  for  numerical  computation 
An  arbitrary  phase  angle 
Dirac  delta  function 
See  page  196 

Partial  time  step  given  by  Eq  (153c) 

Longitudinal  strain  in  the  line  at  s,t 
(defined  by  da/Ss  -1) 

Permanent  strain  (or  passive  analas+.ic  strain) 

A  partial  strain  referred  to  in  the  discussion 
of  the  Maxwell  model 

Percurbation  of  e  from  its  mean  value  at  s 
See;  Eq  (188) 

Ratio  of  lateral  thrust  on  line  to  drag  force 

Zenith  angle  cf  t,  see  Fig  15 

Mean  value  of  0  at  s 

Perturbation  of  G  from  its  mean  at  s 

Real  part  of  the  complex  wave  number  k  for 
transverse  wave  modes 

Potential  translational  displacement  of  the 
surface  vessel 

A 

Cartesian  component  of  X  in  surge,  sway,  and 
heave  (j=l, 2, 3  respectively) 

See  Eq  (107) 


X 
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a  .,a 

tJ  t J 


<*is>082 


Unit  vector  in  the  natural  coordinate  system 
normal  to  r  and  £  at  s,t 

A 

Mean  orientation  of  (j  at  s 

A  A 

Departure  of  ja  from  !Jo 

Imaginary  part  of  the  complex  wave  number  x 
for  transverse  wave  modes 

Unit  vector  in  the  natural  coordinate  system 
normal  to  r  and  at  s,t 

A 

Mean  orientation  of  v 

A  A 

Departure  of  v  from  vQ 

Imaginary  part  of  the  complex  wave  number  for 
longitudinal  waves  (in  the  absence  of  curvature) 

Modified  value  of  v  for  longitudinal  -waves  in 
the  presence  of  curvature  of  the  line 

Active  anelastic  strain  at  s,t  (its  meaning 
is  clarified  by  reference  to  Chap  III,  Sec.  2) 

Mean  value  of  £  at  s 

Perturbation  of  §  from  its  mean  at  s 

3.14159  •  .  • 

Average  density  of  sea  water  (about  2  slugs 
per  ft3 ) 

Actual  arc  length  of  line  between  0  and  P 
in  the  general  strained  state  at  time  t. 

Nominal  natural  frequencies  for  an  unmoored 
vessel  (ax ,  cs,  a3  *  being  zero) 

Variance  of  V^ 1  and  V  *  respectively  at  s 

Unit  vector  tangent  to  line  at  s,  t 

A 

Mean  oreintation  of  T  at  s 
Departure  of  t  from  t0  at  s 
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To , Characteristic  relaxation  times  for  the  line 

A 

$  Azimuth  of  t,  see  Pig.  15 

$  1  Anomaly  of  the  rotational  displacement  of  the . 

surface  vessel 

A 

\|(  . 1  Cartesian  components  of  $  *  in  roll,  pitch, 

J  and  yaw  (j  =  1,2,3  respectively) 

1)1  (uu)  Variance  spectrum  for  V  '  at  s  =  L 

(uj)  Variance  spectrum  for  V  *  at  s  =  L 

♦3(uu)  Variance  spectrum  for  VT '  at  s  =  L 

q. ’  Cartesian  components  of  the  anomaly  of  angular 

J  velocity  of  the  surface  vessel  (i  =  1,2,3 

respectively  for  roll,  pitch,  and  yaw) 

id  Radian  frequency 

uj^  Spectral  mean  value  of  w,  see  Eq  (248) 

u)  Strouhal  frequency 
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